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Statement  of  Problem  Studied 


This  contract  investigated  the  detection  of  electron  and  nuclear  spins  using  the 
technique  of  magnetic  resonance  force  microscopy  (MRFM).  MRFM  is  of  great 
interest  because  of  its  potential  for  atomic-resolution  magnetic  resonance  imaging  of 
molecules  and  materials.  The  primary  goal  of  the  contract  was  to  extend  the 
sensitivity  of  MRFM  to  detect  individual  electron  spins.  The  goal  was  accomplished. 
An  additional  goal  of  the  contract  was  the  development  of  the  technical  infrastructure 
(cantilevers,  magnetic  tips,  signal  processing,  theory,  etc.)  to  support  single  electron 
spin  detection  and  extension  of  MRFM  to  nuclear  spin  detection. 


Scientific  Accomplishments 

1.  Single-spin  magnetic  resonance  force  microscopy  (IBM) 


The  primary  goal  of  this  contract  was  successfully  achieved:  the  demonstration  of 
single  electron  spin  detection  by  magnetic  resonance  force  microscopy  (MRFM).  In 
the  key  experiment,  a  low  temperature  MRFM  apparatus  detected  the  weak 
(attonewton)  magnetic  force  from  an  individual  electron  spin  associated  with  a 
dangling  bond  on  a  silicon  atom  buried  beneath  the  surface  of  a  silicon  dioxide 
sample.  This  demonstration  represents  a  watershed  event  in  the  development  of 
MRFM  since  it  conclusively  demonstrated  that  MRFM  is  indeed  capable  of  single  spin 
detection,  as  previously  proposed  by  John  Sidles  (University  of  Washington). 

For  a  full  description  of  the  single  spin  experiment,  please  see  the  paper  published  in 
the  journal  Nature  (paper  #1  in  the  appendix). 

The  single  spin  demonstration  was  built  upon  a  base  of  science  and  technology 
developed  during  the  course  of  the  contract.  Some  of  the  key  elements  were: 

1 .  Ultrasensitive  cantilevers  that  were  engineered  to  minimize  spin  relaxation  by 
reducing  the  thermal  motion  in  high  order  modes.  See  Appendiix  paper  #2. 

2.  Methods  to  manipulate  and  detect  small,  randomly  polarized  spin  ensembles. 
See  Appendix  papers  #3  and  #4. 

3.  Detailed  understanding  of  spin  relaxation  in  the  context  of  MRFM.  See 
Appendix  paper  #5. 

4.  Signal  processing  techniques  optimized  for  detecting  stochastic  spin  signals. 
See  Appendix  paper  #6. 
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Figure  1.1  -  Diagram  of  the  MRFM  single  spin  experiment.  See  paper  #1  in  Appendix  for  details. 

The  single  spin  accomplishment  was  met  with  considerable  acclaim  in  both  the  main¬ 
stream  and  the  scientific  press.  The  American  Institute  of  Physics  called  the  single 
spin  detection  work  the  top  physics  story  of  2004  (see  the  AIP  web  site: 
http://www.aip.Org/pnu/2004/split/7 1 1  -1  .html).  The  single  spin  paper  was  ranked  #7 
of  the  most  frequently  downloaded  articles  on  the  Nature  web  site  during  July  2004 
and  was  the  only  physical  science  article  to  make  the  top  ten  list. 

News  of  the  single  spin  result  was  widely  disseminated.  Articles  appeared  in  the 
following  media:  New  York  Times,  Dallas  Morning  News  (full  length  feature),  San 
Jose  Mercury  News,  Kansas  City  Star,  Daily  Telegraph  (London),  II  Sole  24  Ore 
(Italy),  Science,  Physics  Today,  New  Scientist,  Discover  magazine,  Science  News, 
Chemical  and  Engineering  News,  and  Analytical  Chemistry  (feature  article).  This  is 
only  a  partial  list. 

2.  Detection  of  statistical  polarization  of  nuclear  spins  in  a  protein  (IBM) 

In  addition  to  our  extensive  work  on  electron  spin  detection,  we  also  investigated 
nuclear  spin  detection.  In  particular,  a  key  experiment  was  performed  that 
demonstrates  that  the  MRFM  measurement  protocols  perfected  for  electron  spin 
detection  also  work  well  with  nuclear  spins. 

In  an  initial  experiment,  the  statistical  polarization  of  fluorine  nuclear  spins  in  calcium 
fluoride  was  measured  using  the  interrupted  OSCAR  protocol,  where  OSCAR  stands 
for  oscillating  cantilever-driven  adiabatic  reversal.  In  a  second  experiment,  signals 
were  obtained  from  a  biomolecule:  the  protein  collagen.  This  result  is  significant  since 
a  key  long  term  goal  of  MRFM  research  is  to  image  the  structure  of  biomolecules. 

This  work  has  been  published  in  Physical  Review  (see  paper  #7  in  the  appendix  for 
technical  details). 
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3.  Image  reconstruction  algorithms  for  atomic  resolution  MRFM  (U. 
Michigan) 

A  major  motivation  for  magnetic  resonance  force  microscopy  (MRFM)  is  the 
possibility  that  three-dimensional  atomic  resolution  of  molecules  can  be  obtained.  To 
recover  atomic  positions  from  two-  and  three-dimensional  MRFM  data,  an  image 
reconstruction  procedure  must  be  used  since  the  MRFM  data  is  a  convolution 
between  the  atomic  positions  and  the  paraboloidal  point  spread  function  given  by  the 
shape  of  the  resonant  slice. 

The  image  reconstruction  problem  has  been  investigated  during  the  course  of  this 
contract  by  the  University  of  Michigan  (UM)  team,  led  by  Prof.  Al  Hero.  The  UM  team 
has  developed  a  sophisticated  algorithm  that  is  effective  at  atomic-scale 
reconstructions  even  when  the  signal-to-noise  ratio  is  low. 

Some  examples  of  the  reconstructions  are  shown  below  in  Figs.  3.1  -  3.3. 


Proton  Imaging  of  Molecules  -  Simulation 

A  simple  example:  Benzene  - 12  atoms,  6  protons 


N — N 
.5 


Proton  positions 


MRFM  x-v  scan 

cimi  ilatinn 


x[nmj 


r,ip=  3  nm 

tip  spacing  ~  3  nm 
tip  vibration  ~  0.5  A 


Figure  3.1  -  Simulation  of  MRFM  x-y  scan  of  protons  in  a  benzene  molecule.  The  six  rings 
originate  from  the  six  protons  in  the  molecule. 
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Reconstruction  of  proton  positions 


Advanced  reconstruction  algorithm  based  on 
expectation  maximization  with  empirical 
Baysian  denoising  (EM-EB) 

M.  Ting,  R.  Raich  and  A.  Hero,  Univ.  of  Michigan 

Figure  3.2  -  Reconstruction  of  proton  positions  from  simulated  MRFM  data  for  a  benzene 
molecule. 

Molecular  Imaging  -  Simulation  &  Reconstruction  m.  Ting,  r.  Raich  &  a.  Hero 

DNA  example:  772  atoms,  one  C13  per  nucleotide  univ.  of  Michigan 


*[nm]  Advanced  reconstruction  algorithm  based  on 


expectation  maximization  with  empirical  Baysian 
denoising  (EM-EB) 


Figure  3.3  -  Reconstruction  of  C-13  position  in  simulation  of  DNA  molecule.  The  molecule  is 
assumed  to  have  been  tagged  with  one  C-13  atom  per  nucleotide. 
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4.  MRFM  Signal  Detection  Studies  (University  of  Michigan) 

The  detection  of  the  single  electron  spin  signal  was  studied  under  the  oscillating 
cantilever-driven  adiabatic  reversal  (OSCAR)  and  interrupted  OSCAR  protocol.  The 
main  results  are  contained  in  paper  #6  attached  to  the  Appendix. 

Initially,  two  continuous-time  (CT)  models  were  studied:  firstly,  the  CT  model  obtained 
by  using  a  classical  dynamic  treatment  of  the  cantilever-single  spin  interaction; 
secondly,  the  CT  model  obtained  assuming  that  the  cyclic  adiabatic  inversion 
condition  holds,  and  modeling  the  electron  spin  decoherence  by  a  random  telegraph 
process.  Throughout  the  work  on  detection,  an  additive  white  Gaussian  noise  (AWGN) 
source  was  assumed.  The  first  CT  model  resulted  in  two  approaches,  which  will  be 
briefly  summarized  here.  Firstly,  a  parametric  approach  to  the  detection  was  explored 
[1 1].  A  variant  of  the  extended  Kalman  filter  (EKF)  was  derived;  together  with  a 
normal  Kalman  filter  (KF),  these  two  filters  formed  the  spin  detector.  This  detector 
was  based  on  comparing  the  residual  prediction  errors  of  the  Kalman  non-linear 
state  predictors  generated  by  the  classical  non-linear  dynamical  equations  of  motion 
of  the  cantilever  under  spin  and  no  spin  conditions,  respectively.  This  detector 
depended  on  knowing  the  exact  values  of  all  the  parameters  in  the  model.  The 
second  approach  was  a  non-parametric  “pattern  recognition”  method  applied  directly 
to  the  phase  space  generated  by  the  dynamical  equations.  This  detector  relied  on 
the  fact  that  under  spin  coupling  the  cantilever  position  oscillates  in  a  different  part  of 
phase  space  than  under  no  spin  coupling  [12].  This  phase  space  was  estimated 
using  correlation  dimension  and  Taken’s  embedding  technique  to  perform  pattern 
recognition  on  the  phase  trajectories.  None  of  the  above  techniques  were  tested  on 
experimental  data.  The  parametric  approach  was  computationally  expensive. 
Moreover,  not  all  of  the  model  parameters  were  exactly  known  in  the  experiment.  The 
non-parametric  approach  seemed  to  work  well  for  only  high  signal  to  noise  ratios 
(SNRs).  In  contrast,  the  experiment  was  conducted  in  a  relatively  low  SNR 
environment. 

We  next  considered  optimal  detection  for  the  random  telegraph  model.  The  first 
result  obtained  was  the  design  of  a  hybrid  Bayes/generalized  likelihood  ratio  (GLR) 
detector  for  the  standard  CT  random  telegraph  model  [13].  A  Gibbs  sampler  was 
used  to  approximate  the  integration  and  maximization  dictated  by  the  hybrid 
Bayes/GLR  principle.  The  major  difficulty  with  this  implementation  was  high 
computational  complexity.  This  led  us  to  study  an  analogous  random  telegraph 
detection  problem  in  discrete-time  (DT),  i.e.,  a  two  state  Markov  chain  derived  from 
discretization  of  the  random  telegraph  signal  in  AWGN.  A  representation  for  the 
optimal  Bayes  detector  was  obtained  and  it  was  shown  in  [6,14]  that  this  detector 
reduced  to  simple  low  complexity  detectors  under  various  relevant  limiting  conditions. 
In  particular,  we  proved  that  the  filtered  energy  (FE)  detector  is  approximately 
optimal  under  the  following  four  conditions:  symmetric  transition  probabilities,  low 
SNR,  long  observation  time,  and  a  small  probability  of  transition  between  two 
consecutive  time  samples.  When  only  the  last  three  conditions  hold,  the  hybrid  FE 
detector,  an  extension  of  the  FE  detector,  was  shown  to  be  approximately  optimal.  An 
alternative  model  that  was  proposed  for  the  single  electron  spin-cantilever  interaction 
is  the  DT  random  walk  model.  Simulations  were  presented  in  [6]  showing  that,  for 
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certain  random  walks,  the  FE  detector  is  approximately  optimal.  Theoretical  analysis 
is  presented  in  a  paper  under  preparation. 

Under  the  conditions  of  low  SNR  and  assuming  an  AWGN  noise  source,  it  was 
shown  that  the  optimal  detector  test  statistic  for  a  finite-state  Markov  signal  is 
approximately  the  matched  filter  statistic  with  the  minimum  mean  squared  error 
(MMSE)  predictor  values  used  in  place  of  the  known  signal  values  [6].  This  result 
provides  a  general  means  of  obtaining  a  suboptimal  detector  under  low  SNR  for  any 
finite-state  Markov  signal.  When  applied  to  the  DT  random  telegraph  model,  one 
obtained  an  estimator-correlator  detector.  By  using  a  suboptimal  linear  predictor 
instead  of  the  MMSE  predictor,  the  hybrid  FE  detector  was  recovered.  The  FE 
detector  was  successfully  used  in  the  detection  of  the  single  electron  spin. 


5.  Digital  signal  processing  and  control  (University  of  Washington) 


Work  at  the  UW  focused  on  system  control  and  end-to-end  simulation  in  MRFM.  This 
work  was  completed  as  planned,  culminating  in  a  PhD  defense  by  UW  graduate 
student  Tom  Kriewal  titled  “Heterodyne  Digital  Control  and  Frequency  Estimation  in 
Magnetic  Resonance  Force  Microscopy”.  This  thesis  documents  the  main  UW 
deliverables,  which  are  two  Pentek  VME  crates,  each  capable  of  simultaneous 
control  of  a  Larmor  cantilever,  and  emulation  of  that  cantilever. 

Both  program  tasks  -  -  emulation  and  control  -  -  are  running  in  closed-loop  mode. 

The  performance  of  the  system  is  in  good  accord  with  the  original  performance  goals. 
This  design  allowed  the  investigators  to  test  a  novel  method  of  heterodyne  control  for 
25,  50,  and  100  MHz  signals.  Heterodyne  control  is  a  method  for  reducing 
computational  load  and  noise  outside  the  passband  and  generating  lock-in 
(synchronous  detection)  signals  for  online  diagnostics,  system  identification,  and 
adaptive  control.  Using  one  crate  to  generate  the  signal,  and  the  other  to  control  it,  is 
the  final  deliverable  of  a  control  system  that  will  work  with  cantilever  signals  at  these 
frequencies. 


6.  Quantum  measurement  theory  and  simulation  (Univ.  of  Washington) 

Work  performed  at  Univ.  of  Washington  has  established  a  reasonably  satisfactory 
theoretical  understanding  of  the  IBM  single-spin  iOSCAR  experiments.  Both  quantum 
measurement  theory  and  IBM's  experimental  data  agree  that  these  experiments  can 
be  accurately  described  as  random  telegraph  signals  embedded  in  white  noise 
(RTS+WN).  This  description  is  valid  for  all  signal-to-noise  ratios,  and  there  are 
reasonable  theoretical  grounds  to  expect  that  this  RTS+WN  property  will  characterize 
future  single-spin  Larmor  experiments. 

Note  also  that  RTS+WN  processes  have  mathematical  properties  that  are  very 
desirable  from  an  engineering  and  signal  processing  point  of  view.  In  particular, 
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numerical  models  of  such  process  are  highly  efficient  and  straightforward  to 
implement,  and  theoretical  understanding  of  the  informatic  content  of  such 
processes  is  well  advanced.  The  significance  of  this  theoretical  finding  is  that  real¬ 
time  quantum  simulation  has  now  been  shown  to  be  feasible  for  MRFM  experiments. 

The  results  show  that  using  a  semi-classical  formalism  to  analyze  iOSCAR 
experiments  can  be  rigorously  justified  via  orthodox  POVM  theory,  and  that  Carlton 
Cave’s  amplifier  noise  bounds  can  be  applied  to  MRFM  experiments.  We  have  also 
shown  that  spin-spin  interactions  are  possible  for  Larmor  experiments  in  a  biological 
spin  environment,  provided  the  spin-lock  field  exceeds  20-30  Gauss  or  so,  possibly 
less. 


7.  Quantum  simulations  of  dipolar  coupled  spins  in  a  field  gradient  (GE 
Medical) 

Magnetic  Resonance  Force  Microscopy  (MRFM),  as  first  conceived  by  Sidles, 
involves  the  coupling  of  a  spin  and  a  cantilever,  both  resonant  at  the  Larmor 
frequency  of  the  spin.  In  subsequent  work,  it  has  proven  practical  to  resonate  the 
cantilever  at  a  frequency  small  compared  to  that  of  the  spin,  and  to  effect  a  coupling 
between  the  two  by  means  of  a  modulation  scheme,  such  as  adiabatic  passage, 
which  causes  the  spin  to  cyclically  reverse  at  some  low  frequency,  matching  that  of 
the  cantilever.  Such  a  strategy  is  also  convenient  for  calculations,  since  it  obviates 
(to  some  approximation)  the  necessity  of  transforming  the  cantilever  spatial 
coordinates  to  an  interaction  frame  at  the  Larmor  frequency. 

That  is,  the  oscillatory  field  of  the  cantilever,  in  any  direction  transverse  to  the 
quantization  axis,  will  produce  no  transitions  of  the  spin  system,  while  that  along  the 
axis  will  produce  at  most  a  frequency  or  phase  modulation.  We  have  performed 
calculations  of  the  spin  dynamics  -  focusing  particularly  on  dipole-mediated  spin 
exchange  in  the  presence  of  the  static  gradient  from  the  cantilever  tip--  for  simple 
some  models  in  this  regime,  and  with  an  eye  towards  the  following  question:  What  is 
the  likelihood  of  quenching  spin  diffusion  with  a  realizable  gradient? 


The  Two-Spin  Problem 

For  a  pair  of  dipolar-coupled  protons  in  a  static  B  field,  consisting  of  a  uniform 
component,  plus  a  strong  linear  gradient,  the  Hamiltonian  (in  units  of  radian/sec)  is: 

Hf  <%  +  k»z\  >z2  ~  >W+]  >-2  +  42  4H  • 
where  Jis  gradient  strength  times  the  internuclear  distance,  k  is  the  dipolar  constant, 

o  o  o 

Pq/(4ti)  y  h  r^2  (3cos‘  0-  1 ),  with  <9  the  angle  between  the  inter-proton  vector  and 

the  quantization  axis,  and  we  have  chosen  an  interaction  representation  in  which  the 
static  field  vanishes  at  one  of  the  protons. 

The  Multi-Spin  Problem 
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The  Hamiltonian  for  a  homonuclear  spin  cluster  is  given  by: 


H  =  2 :«*<r  i)Izi  +  2  kij  Uz  i  Izj  ~  1/4 (/+  f  /y  +  /+  y  /  /)] 


with  <*y(r  y)  as  the  field  of  the  cantilever  tip  at  the  /'th  spin,  and  k  jj  equal  to 

K/(47t:)  Y r  jj~^(3cos^0  jj  -  1).  This  admits  of  no  simple  solution  in  closed  form, 

although  it  is  straightforwardly  simulated.  We  have  employed  the  following  models: 
four  spins  in  a  straight  line,  of  arbitrary  direction,  with  internuclear  separation  of  2.5  A, 
and  seven  spins  in  a  cubic  lattice,  comprising  a  central  spin  surrounded  by  all  six  of 
its  nearest  neighbors,  also  with  internuclear  separation  of  2.5  A.  The  cantilever  tip  is 
modelled  as  a  sphere  of  50  A  radius,  with  juqM  =  1 .7  T. 

The  computation  consists  of  numerical  diagonalization  of  the  Hamiltonian,  followed 
by  evolution  of  the  lz  operators,  for  given  initial  conditions,  which  are  chosen,  given 

the  small  number  of  spins,  for  illustrative  purposes,  rather  than  on  physical  grounds. 

Note  that  we  follow  the  usual  practice  for  high  field  NMR  in  solids,  and  use  only  the  z 
component  of  the  dipole  field  produced  by  the  magnetic  tip;  this  is  easily  justifiable  for 
gradient  fields  (in  frequency  units)  <  1  %  of  the  Larmor  frequency. 

The  results  of  the  calculations  are  given  in  the  figures  below.  The  main  conclusion  of 
the  work  is  that  the  strong  field  gradients  needed  for  atomic  resolution  imaging  will 
suppress  the  flip-flop  interactions  between  dipolar-coupled  spin,  except  along  the  iso¬ 
field  lines  of  the  magnetic  tip.  In  general,  the  suppression  of  dipolar  coupling  effects  is 
believed  to  be  beneficial  for  atomic  resolution  imaging.  Further  work  will  be  required 
to  gauge  the  impact  of  the  iso-field  regions  where  there  is  inefficient  suppression. 

This  may  not  be  a  serious  issue  since,  during  imaging,  the  tip  is  mechanically 
scanned  with  respect  to  the  sample,  which  implies  that  dipolar  non-suppression  is 
only  a  temporary,  localized  condition. 


10 


Exemplary  calculation  of  flux  lines  for 
the  magnetic  tip;  axes  in  A 


Iso-contours  of  z-directed  dipole  field  of 
the  magnetic  tip,  in  units  of  Hz  (  proton 
larmor  frequency/27i).  Note  the 
characteristic  “d  orbital”  shape,  and  the 
nodal  lines  along  the  magic  angle 
(direction  cosines  [1  1  1]).  The  vertical  o 
and  horizontal  axes  are  z  and  y;  axes  in  A 


Figure  7.1:  Characteristics  of  the  model 
magnetic  tip. 


Dipole  z-field  strength  along  z  axis, 
in  Hz  (as  above)  axes  in  A 


Figure  7.2 :  Dynamics  of  linear  spin 
arrays  parallel  and  perpendicular  to 
maximum  gradient  direction 


■400  ‘300  ‘200  -100  0  100  200  3W  400  5ffl» 


Position  of  four-  spin  array  on  the 
dipople  field  iso-contour  map. 
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0 

MM 


Time  course  of  spin  evolution  (over  0.5 
ms)  for  four  spins  aligned  in  direction  of 
maximum  gradient  (~35  G/A).  The  initial 
condition  is  color-coded  in  the  transpierced 
circles  at  right;  note  total  suppression  of 
spin  exchange/diffusion. 


Time  course  of  spin  evolution  (over  0.5  ms)  for  four 
spins  aligned  approximately  on  tip-field  iso-contour 
(minimum  gradient).  The  initial  condition  is  per  color- 
coded  circles  below;  spin  exchange/diffusion  proceeds 
virtually  unhindered. 
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Single  spin  detection  by  magnetic 
resonance  force  microscopy 
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Magnetic  resonance  imaging  (MRI)  is  well  known  as  a  powerful 
technique  for  visualizing  subsurface  structures  with  three- 
dimensional  spatial  resolution.  Pushing  the  resolution  below 
1  pirn  remains  a  major  challenge,  however,  owing  to  the  sensitivity 
limitations  of  conventional  inductive  detection  techniques.  Cur¬ 
rently,  the  smallest  volume  elements  in  an  image  must  contain  at 
least  10 12  nuclear  spins  for  MRI-based  microscopy1,  or  107 
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electron  spins  for  electron  spin  resonance  microscopy2.  Magnetic 
resonance  force  microscopy  (MRFM)  was  proposed  as  a  means  to 
improve  detection  sensitivity  to  the  single-spin  level,  and  thus 
enable  three-dimensional  imaging  of  macromolecules  (for  ex¬ 
ample,  proteins)  with  atomic  resolution3,4.  MRFM  has  also  been 
proposed  as  a  qubit  readout  device  for  spin-based  quantum 
computers5,6.  Here  we  report  the  detection  of  an  individual 
electron  spin  by  MRFM.  A  spatial  resolution  of  25  nm  in  one 
dimension  was  obtained  for  an  unpaired  spin  in  silicon  dioxide. 
The  measured  signal  is  consistent  with  a  model  in  which  the  spin 
is  aligned  parallel  or  anti-parallel  to  the  effective  field,  with  a 
rotating-frame  relaxation  time  of  760  ms.  The  long  relaxation 
time  suggests  that  the  state  of  an  individual  spin  can  be  mon¬ 
itored  for  extended  periods  of  time,  even  while  subjected  to  a 
complex  set  of  manipulations  that  are  part  of  the  MRFM 
measurement  protocol. 

MRFM  is  based  on  the  detection  of  the  magnetic  force  between  a 
ferromagnetic  tip  and  spins  in  a  sample.  The  fundamental  challenge 
in  achieving  single-spin  sensitivity  is  that  the  force  from  a  single 
spin  is  exceedingly  small.  Even  with  tip  field  gradients  in  the  gauss 
per  nanometre  range,  the  force  from  an  electron  spin  is  only  a  few 
attonewtons.  This  force  is  roughly  a  million  times  smaller  than  is 
typically  detected  by  atomic  force  microscopy  (AFM).  Recently, 
major  strides  towards  single-spin  detection  have  been  made  with 
the  development  of  ultrasensitive  cantilever-based  force  sensors7,8, 
better  understanding  of  relevant  spin  relaxation  processes9-12  and 
the  successful  detection  of  statistical  polarization  in  small  spin 
ensembles13. 

The  basic  elements  of  our  MRFM  apparatus  are  shown  in  Fig.  1 
and  have  been  described  in  detail  in  ref.  13.  Briefly,  a  custom- 
fabricated  mass-loaded  silicon  cantilever8,13  with  an  attached  150- 
nm-wide  SmCo  magnetic  tip  is  used  to  sense  the  force  from  the 
electron  spin.  The  sample  consists  of  vitreous  silica  (Suprasil  W2) 
that  was  irradiated  with  a  2 -Gy  dose  of  Co60  gamma  rays.  The 
gamma  irradiation  produces  a  low  concentration  of  silicon  dangling 
bonds  containing  unpaired  electron  spins  known  as  E7  centres14. 


Cantilever 


Microwave 

Interferometer  ■  ^  coil 


Figure  1  Configuration  of  the  single-spin  MRFM  experiment.  The  magnetic  tip  at  the  end 
of  an  ultrasensitive  silicon  cantilever  is  positioned  approximately  1 25  nm  above  a  polished 
Si02  sample  containing  a  low  density  of  unpaired  electron  spins.  The  resonant  slice 
represents  those  points  in  the  sample  where  the  field  from  the  magnetic  tip  (plus  an 
external  field)  matches  the  condition  for  magnetic  resonance.  As  the  cantilever  vibrates, 
the  resonant  slice  swings  back  and  forth  through  the  sample  causing  cyclic  adiabatic 
inversion  of  the  spin.  The  cyclic  spin  inversion  causes  a  slight  shift  of  the  cantilever 
frequency  owing  to  the  magnetic  force  exerted  by  the  spin  on  the  tip.  Spins  as  deep  as 
1 00  nm  below  the  sample  surface  can  be  probed. 
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Estimated  spin  concentration  was  between  1013  and  1014  cm-3.  The 
experiments  were  performed  at  1.6  K  in  a  small  vacuum  chamber 
that  fits  within  the  bore  of  a  superconducting  magnet.  The  low 
operating  temperature  minimizes  the  force  noise  and  reduces  the 
relaxation  rate  of  the  spins. 

A  microwave  magnetic  field  ( B  j  ~  0.3  mT),  in  combination  with 
the  inhomogeneous  field  from  the  magnetic  tip,  sets  up  a  ‘resonant 
slice’  within  the  sample.  The  resonant  slice  corresponds  to  those 
points  in  the  sample  where  the  tip  field  Btip(x,  y,  z),  plus  a  static 
external  field  z  Bext,  satisfies  the  condition  for  electron  spin 
resonance:  B0(x,/,z)  =  |Btip(x,y,z)  +  z  Bextl  =  <*>rf/Y-  Here,  corf 
is  the  frequency  of  the  microwave  field  and  y  is  the  gyromagnetic 
ratio  (y/l'K  =  2.8  X  10loHzT-1).  In  the  present  experiment,  where 
corf/2 tt  =2.96  GHz,  the  resonant  slice  corresponds  to  B0(x,  y, 
z)  =  106  mT.  For  typical  conditions  (for  example,  Bext  =  30  mT), 
the  slice  is  a  bowl-shaped  surface  that  extends  roughly  250  nm  below 
the  tip.  Note  that  with  the  perpendicular  cantilever  orientation 
shown  in  Fig.  1 ,  the  cantilever  is  responsive  only  to  the  x  component 
of  force.  Thus  the  spin  must  be  located  either  slightly  in  front  of  or 
slightly  behind  the  cantilever  in  the  x  direction  to  produce  a 
substantial  response. 

To  generate  a  force  signal  that  can  be  distinguished  from  the 
much  larger  background  force  fluctuations,  we  use  the  recently 
developed  spin  manipulation  protocol  known  as  ‘interrupted 
OSCAR’  or  iOSCAR,  where  OSCAR  stands  for  oscillating  canti¬ 
lever-driven  adiabatic  reversals13.  In  the  iOSCAR  protocol,  the 
cantilever  is  part  of  a  gain-controlled  positive-feedback  loop  that 
drives  the  cantilever  to  oscillate  at  a  set  amplitude  (for  example, 
16  nm)  at  the  fundamental  frequency  of  the  cantilever 
( fc  =  5.5  kHz).  Because  the  cantilever  is  the  frequency- determining 
element  in  the  feedback  loop,  the  vibration  frequency  will  auto¬ 
matically  vary  in  response  to  tip-sample  interactions15. 

The  vibration  of  the  cantilever  tip  causes  the  resonant  slice  to 
sweep  back  and  forth  rapidly  through  the  sample.  If  the  slice 
happens  to  sweep  through  the  location  of  a  spin,  the  spin  will  be 
cyclically  inverted  in  synchrony  with  the  cantilever  motion  owing  to 
the  phenomenon  of  adiabatic  rapid  passage13,16.  This  synchronous 
inversion  of  the  spin  creates  an  alternating  magnetic  force  on  the 
cantilever  that  mimics  a  change  in  cantilever  stiffness.  The  resulting 
shift  in  cantilever  frequency  is  given  by13,17: 


^  =  +  2/cG/*b 
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Figure  2  Timing  diagram  for  the  iOSCAR  spin  manipulation  protocol.  The  zcomponent  of 
spin  follows  the  motion  of  the  cantilever  except  during  the  interruption  of  the  microwave 
field.  The  interruptions  last  for  one-half  of  a  cantilever  period  and  are  precisely  timed  to 
start  and  stop  at  vibration  extrema.  The  interruptions,  which  occur  every  64  cantilever 
cycles,  reverse  the  phase  of  the  spin  with  respect  to  the  cantilever,  resulting  in  a  square- 
wave  oscillation  of  the  cantilever  frequency. 


where  k  is  the  cantilever  spring  constant,  xpeak  is  the  peak  vibration 
amplitude  of  the  cantilever,  G  =  dB0/dx  is  the  lateral  field  gradient, 
and  /jlb  is  the  magnetic  moment  of  the  electron  (9.3X  10-24JT-1). 
The  sign  of  the  frequency  shift  depends  on  the  relative  phase  of  the 
spin  inversions  with  respect  to  the  cantilever  motion.  Using  the 
terminology  of  magnetic  resonance,  the  two  polarities  correspond 
to  adiabatic  rapid  passages  with  the  spin  either  aligned  or  anti¬ 
aligned  with  respect  to  the  effective  field  in  the  rotating  frame13,16. 
This  parallel  or  anti-parallel  alignment  is  expected  to  be  enforced 
by  the  quantum-mechanical  measurement  process,  resulting  in  the 
collapse  of  the  initial  spin  wavefunction  onto  an  eigenstate  of  the 
effective  field18,19.  For  the  parameters  of  the  current  experiment 
(G  =  2  X  105  Tm"1;  k  =  0.1 1  mN m-1;  xpeak=  16nm),  the 
expected  frequency  shift  is  \hfc\  =  3.7  ±  1.3 mHz.  The  estimated 
uncertainty  in  \5fc\  reflects  20%  uncertainties  in  the  calibration  of 
G,  k  and  xpeak. 

In  the  ‘interrupted’  version  of  the  protocol  (iOSCAR),  the 
microwave  field  B  x  is  turned  off  for  one-half  of  a  cantilever  cycle 
every  64  cycles  ( /int  =  fc/6 4  ~  86  Hz).  As  shown  in  Fig.  2,  each  time 
B  i  is  interrupted,  the  relative  phase  of  the  spin  and  cantilever  is 
reversed,  causing  the  frequency  shift  to  reverse  polarity.  (See 
Supplementary  Information  for  an  animation.)  The  net  result  of 
the  interruptions  is  that  the  frequency  shift  will  alternate  between 
positive  and  negative  values  in  a  square-wave-like  fashion  with  a 
frequency  given  by  /sig  =  /int/2,  or  approximately  43  Hz.  The 
resulting  frequency  shift  signal  can  be  written  as  the  Fourier  series: 

4 

A/(f)  =  —  |6/c |A(f)sin(2ir/sigf)  +  higher  harmonics  (2) 

TT 

where  the  4/tt  comes  from  the  first  harmonic  Fourier  amplitude  of  a 
square  wave. 

We  have  included  the  function  A(t)  to  take  into  account  the  fact 
that  the  signal  will  not  be  perfectly  periodic  owing  to  extra  random 
spin  flips  induced  by  the  environment.  Assuming  that  this  ‘quan¬ 
tum  jump’  model  is  correct,  A(t)  is  a  random  telegraph  function 
that  takes  on  values  of  ±1.  For  Poisson  distributed  jumps,  A(t)  is 
expected  to  have  a  lorentzian  power  spectrum20  and  the  statistical 
properties  (A(t))  =  0  and  ([A(f)]2)  =  1,  where  (•••)  denotes  time 
average.  We  detect  the  first  harmonic  of  the  signal  only,  so  we  will 
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Figure  3  Plots  showing  the  spin  signal  as  the  sample  was  scanned  laterally  in  the  x 
direction  for  two  values  of  external  field:  a,  Fext  =  34  mT,  and  b,  BeyX  =  30  mT.  The 
smooth  curves  are  gaussian  fits  that  serve  as  guides  to  the  eye.  The  1 9-nm  shift  in  peak 
position  reflects  the  movement  of  the  resonant  slice  induced  by  the  4-mT  change  in 
external  field.  The  difference  in  absolute  peak  height  is  primarily  due  to  different  lock-in 
amplifier  detection  bandwidths:  0.1 8  Hz  and  0.59  Hz  for  a  and  b,  respectively.  Power 
spectra  for  the  points  marked  1  and  2  are  shown  in  Fig.  4. 
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refer  to  the  quantity  A/i(t)  =  (4/Tr)|5/c|A(f)  as  the  spin  signal 
amplitude. 

Because  the  frequency  modulation  due  to  the  spin  is  only  a  few 
millihertz  and  the  frequency  noise  of  the  cantilever  due  to  thermal 
motion  and  tip-sample  interactions  is  relatively  large  (—25  mHz  in 
a  1-Hz  bandwidth),  we  must  use  signal  averaging  to  detect  the  spin 
signal.  However,  because  (Afi(t))  =  0,  we  average  the  square  of  the 
signal  (the  signal  ‘energy’),  rather  than  the  signal  amplitude. 
Specifically,  the  frequency  modulation  of  the  cantilever  is  detected 
using  an  analogue  frequency  discriminator15  followed  by  a  digital 
lock-in  amplifier  that  has  been  implemented  in  software.  The  lock- 
in  amplifier  uses  a  bank  of  low-pass  filters  and,  as  a  function  of 
detection  bandwidth,  determines  the  energy  (that  is,  variance)  of 
the  in-phase  and  quadrature  components  of  the  frequency- shift 
signal  Af(t).  The  spin  signal  and  the  measurement  noise  are 
uncorrelated,  so  we  can  write  the  in-phase  variance  as  a]  = 
°spin  +  anoise?  where  a2pin  is  the  variance  due  to  the  spin  signal 
and  a2oise  represents  the  measurement  noise.  This  in-phase  variance 
is  then  compared  to  the  quadrature  variance  Oq,  which  contains 
only  the  measurement  noise.  The  signal  energy  from  the  spin  can 
then  be  estimated  as  a2pin  =  o\  —  Oq.  This  energy  detection  scheme 


a 


b 
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Figure  4  By  measuring  the  spin  signal  energy  as  a  function  of  detection  bandwidth,  the 
power  spectral  density  of  the  spin  signal  amplitude  AG  (0  can  be  determined,  a,  Spectral 
density  for  the  two  positions  indicated  in  Fig.  3.  The  strong  spin  signal  at  position  1  has  a 
narrow  spectral  width  (0.21  Hz),  reflecting  the  long  (760  ms)  correlation  time  obtained 
using  iOSCAR  spin  manipulation.  The  data  are  well  fitted  by  a  lorentzian  function  (solid 
curve).  At  position  2,  approximately  70  nm  away  from  the  position  of  the  spin  signal,  the 
spectral  density  is  negligible,  b,  False-colour  plot  showing  power  spectral  density  as  a 
function  of  position.  The  spin  signal  is  localized  both  spatially  and  spectrally.  The  data 
were  interpolated  between  discrete  measurement  positions  to  create  a  smooth  image. 
The  apparent  fine  structure  in  the  signal  region  is  probably  an  artefact  of  the  limited 
signal-to-noise  ratio  of  the  original  data. 
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is  believed  to  have  performance  close  to  the  theoretically  optimal 
detector21. 

Figure  3a  shows  a  lateral  scan  where  we  plot  a2  •  as  a  function  of 
sample  position.  The  scan  shows  a  prominent  pealc  that  we  believe  is 
due  to  a  single  spin.  The  peak  width  is  25  nm,  or  roughly  1.6  times 
the  cantilever  oscillation  amplitude.  Note  that  the  baseline  of 
the  data  on  either  side  of  the  peak  is  essentially  zero  (within  the 
uncertainty  of  the  measurement),  thus  demonstrating  that  the 
iOSCAR  protocol  generates  no  systematic  feedthrough  artefacts. 
Because  the  signal-to-noise  ratio  was  so  low  (cTgpin/(7j  —  0.06), 
considerable  averaging  was  required.  In  Fig.  3a,  the  averaging 
time  was  13  h  per  point,  yielding  a  signal  peak  that  is  five  standard 
deviations  above  the  baseline  noise. 

To  confirm  that  the  observed  signal  is  truly  due  to  magnetic 
resonance,  a  number  of  basic  tests  were  performed.  As  expected  for 
an  iOSCAR  magnetic  resonance  signal,  the  signal  disappeared  if  the 
microwaves  were  absent  or  turned  on  continuously.  The  timing  of 
the  microwave  interruptions  was  also  varied.  The  signal  disap¬ 
peared,  as  expected,  when  the  starting  time  of  the  interruption  was 
shifted  from  the  vibration  peak  to  the  zero-crossing  of  the  vibration 
(that  is,  shifted  by  one-quarter  of  a  cantilever  cycle)  and  when  the 
interruption  duration  was  a  full  cantilever  cycle  instead  of  one-half 
cycle. 

One  key  test  for  magnetic  resonance  is  to  observe  the  field 
dependence  of  the  spin  signal.  If  the  external  field  is  reduced,  the 
resonant  slice  will  shrink  in  radius,  thus  shifting  the  scan  position  of 
the  signal  peak.  When  Bext  was  reduced  from  34  to  30  mT,  the 
expected  peak  shift  was  indeed  observed  and  found  to  be 
Ax  =  19  nm,  as  can  be  seen  by  comparing  the  scans  in  Figs  3a  and 
b.  The  ratio  ABext/Ax  suggests  that  the  field  gradient  G  is  approxi¬ 
mately  2  X  105Tm-1  (2Gnm_1). 

The  conclusion  that  the  signal  is  due  to  only  one  spin  is  based 
primarily  on  the  spatial  isolation  of  the  spin  signal.  By  design,  the 
spin  density  of  the  sample  was  very  low,  in  the  range  of  1013  to 
1014  cm-3,  giving  a  mean  spacing  between  spins  in  the  range  of  200 
to  500  nm.  The  sparseness  of  spins  implies  that,  for  most  sample 
locations,  there  is  no  spin  interacting  with  the  resonant  slice.  This  is 
why  the  data  in  Fig.  3  has  a  zero  baseline.  To  locate  a  spin  signal,  the 
sample  was  scanned  through  many  independent  locations,  of  the 
order  of  30,  before  a  strong  signal  from  a  well-positioned  spin  was 
found. 

By  measuring  a2pin  as  a  function  of  detection  bandwidth,  the 
power  spectral  density  of  the  spin  signal  amplitude  A/X(t)  can  be 
determined  (Fig.  4a).  Consistent  with  a  random  telegraph  model  of 
A/i(f),  the  spectrum  is  well  fitted  by  the  lorentzian  function  S(f)  = 
4rm([A/i(f)]2)/  [l  +  4tt2t^/2]  .  The  spectral  width  at  half-maxi¬ 
mum  is  0.21  Hz,  corresponding  to  an  impressively  long  rm  corre¬ 
lation  time  of  760  ms.  This  is  essentially  the  rotating  frame 
relaxation  time,  but  also  includes  possible  effects  of  the  iOSCAR 
spin  manipulation.  The  long  correlation  time  implies  that  the 
cantilever- driven  spin  inversions  are  coherent  for  thousands  of 
cycles.  The  false-colour  image  in  Fig.  4b  shows  that  the  spin  signal 
is  highly  localized  both  spatially  and  spectrally. 

The  total  magnitude  of  the  spin  signal,  obtained  by  integrating 
the  spectrum  in  Fig.  4a,  was  found  to  be  ([Afi(t)]2)  =  28mHz2. 
Using  this  experimental  value  and  the  relationship  ([A/i(f)]2)  = 
(4/7r)2|6/c|2([A(f)]2),  we  solve  for  \5fc\  using  the  assumption  that 
([A(t)]2)  =  1.  We  find  \5fc\  =  4.2  mHz,  in  excellent  agreement  with 
the  value  of  3.7  mHz  expected  from  equation  (1). 

In  conclusion,  we  have  presented  evidence  that  MRFM  is  now 
capable  of  detecting  individual  electron  spins.  Although  several 
other  single-spin  detection  methods  have  been  previously  demon¬ 
strated22-28,  MRFM  has  some  attributes  that  set  it  apart.  Perhaps  the 
most  important  of  these  is  the  ability  to  image  spins  below  the 
surface  with  nanometre  spatial  resolution.  Spins  as  deep  as  100  nm 
should  be  accessible  under  present  operating  conditions.  Although 
extensive  signal  averaging  is  currently  required,  even  a  modest 
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increase  in  field  gradient  (for  example,  five  times  larger)  will 
dramatically  speed  up  the  acquisition  time  and  thereby  enable 
two-  and  three-dimensional  imaging  applications.  (Because  we 
average  signal  energy,  rather  than  amplitude,  the  averaging  time 
will  decrease  as  the  inverse  fourth  power  of  the  gradient  in  the  limit 
of  low  signal-to-noise  ratio.)  If  the  measurement  time  can  be 
reduced  below  rm,  real-time  readout  of  the  spin  quantum  state 
will  become  possible,  enabling  a  wide  variety  of  quantum  measure¬ 
ment  experiments.  For  molecular  imaging  applications,  extension 
to  single  nuclear  spin  detection  is  necessary,  but  this  will  require  a 
roughly  1,000-fold  improvement  in  magnetic  moment  sensitivity. 
Considering  that  the  present  experiment  represents  a  sensitivity 
improvement  of  107  times  over  the  original  MRFM  experiment29, 
the  remaining  required  improvement  does  not  seem  out  of  the 
question,  especially  since  there  is  still  considerable  leeway  for 
increasing  the  field  gradient  and  lowering  the  operating 
temperature.  □ 
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Bulk  metallic  glasses  are  commonly  produced  by  the  rapid  cool¬ 
ing  of  liquid  alloys1.  They  have  emerged  over  the  past  decade  as  a 
novel  class  of  materials,  with  attractive  properties  and  techno¬ 
logical  promise1,2.  The  bulk  metallic  glasses  so  far  produced 
contain  three  or  more  component  elements3,4.  These  complex 
compositions  are  necessary  to  frustrate  the  crystallization  of  the 
liquid  melt  on  cooling,  but  can  also  lead  to  phase  separation, 
which  is  detrimental  to  the  thermal  and  mechanical  properties  of 
metallic  glasses5-8.  Here  we  report,  using  X-ray  diffraction 
measurements,  the  formation  of  a  bulk  metallic  glass  from 
elemental  zirconium  at  high  static  pressures  and  low  tempera¬ 
tures  (relative  to  its  melting  temperature  at  atmospheric  press¬ 
ure).  Amorphous  zirconium  can  be  recovered  at  ambient 
conditions  and  demonstrates  a  superior  thermal  stability  com¬ 
pared  to  amorphous  alloys3,9,  which  could  lead  to  new  high- 
temperature  applications  of  amorphous  metals. 

In  multi-component  systems,  glass-forming  ability  (GFA)  is 
viewed  as  the  resistance  to  precipitation  of  crystalline  phases  from 
supercooled  liquid  metals10,  and  alloys  with  high  GFA  all  have  three 
common  features3,4:  (1)  they  consist  of  at  least  three  components; 
(2)  there  is  significant  mismatch  of  the  atomic  size  of  the  constitu¬ 
ent  elements;  and  (3)  there  are  negative  heats  of  mixing  among  the 
major  alloying  elements.  Addition  of  elements  that  are  chemically 
and  topologically  different  from  the  other  species  not  only  creates 
an  energy  barrier  for  nuclei  to  form  but  also  effectively  increases 
melt  viscosity  or  fragility,  which  results  in  a  reduced  rate  of  both 
nucleation  and  growth  and  an  increase  in  GFA.  The  production  of 
bulk  glassy  materials  in  pure  metals,  however,  remains  a  long¬ 
standing  scientific  curiosity  and  technological  interest.  The  diffi¬ 
culties  arise  from  the  facts  that  the  equilibrium  melt  viscosity  of 
pure  metals  is  three  orders  of  magnitude  smaller  than  that  of 
amorphous  alloys11  and  that  current  technology  has  yet  to  reach  a 
cooling  rate  in  excess  of  the  10 10  °C  s-1  that  is  needed  to  make  pure 
metals  amorphous12. 

We  studied  zirconium  metal  at  pressures  and  temperatures  up  to 
17GPa  and  1,000  °C,  using  energy-dispersive  synchrotron  X-ray 
diffraction  and  time-of-flight  neutron  scattering.  In  X-ray  diffrac¬ 
tion  experiments,  we  used  both  DIA-type13  and  T-cup14  large- 
volume  high-pressure  apparatus  installed  at  Brookhaven  and 
Argonne  National  Laboratories.  Neutron  scattering  experiments 
were  performed  using  a  high-pressure/high-temperature  (high 
P-T )  cell  assembly15  in  a  TAP-98  toroidal -anvil  press  at  Los  Alamos 
Neutron  Scattering  Center.  The  starting  sample  of  zirconium  has  a 
close-packed  hexagonal  structure  (a -phase)  and  is  of  extremely 
high  purity,  with  35  p.p.m.  Hf,  less  than  25  p.p.m.  of  C,  N  and  Al, 
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ABSTRACT 


INTRODUCTION 


We  describe  the  design,  fabrication  and  testing  of  mass- 
loaded  cantilevers  for  electron-spin  magnetic  resonance 
force  microscopy.  These  single-crystal  silicon  cantilevers 
are  designed  to  have  large  gaps  in  their  thermal  mode 
spectra  so  as  to  reduce  thermal  noise  near  the  electron-spin 
Rabi  frequency.  Each  cantilever  typically  consists  of  a  2  pm 
thick  mass  suspended  at  the  end  of  a  0.1  pm  thick  hinge. 
The  fabrication  process  starts  with  an  SOI  wafer  followed  by 
selective  silicon  epitaxy,  cantilever  patterning,  and  backside 
release.  The  focus  of  the  process  is  on  precise  thickness 
control  and  material  homogeneity. 

We  will  present  characterization  results  for  these  canti¬ 
levers  at  room  and  low  temperatures  and  discuss  the  impact 
of  these  devices  on  magnetic  resonance  force  microscopy 
measurements. 


Magnetic  resonance  force  microscopy  (MRFM)  [1]  uses 
an  ultra-sensitive  cantilever  to  detect  magnetic  resonance  in 
small  ensembles  of  spins  (either  nuclear  spins  or  electron 
spins).  The  technique  requires  the  detection  of  a  very  small 
(attonewton)  oscillating  magnetic  force  exerted  on  a  mag¬ 
netic  tip  by  the  spins  in  the  sample  [2,  3].  An  RF  magnetic 
field  is  used  to  manipulate  the  spin  orientation  at  a  rate  that 
matches  the  kHz  resonance  frequency  of  the  lever  (Fig.  1). 

To  extend  this  technique  to  single-spin  sensitivity,  the 
magnetic  tip  must  generate  a  very  large  magnetic  gradient 
(e.g.,  several  gauss  per  angstrom)  in  order  to  create  a 
measurable  force.  Previous  experimental  and  theoretical 
studies  have  shown,  however,  that  sub -angstrom  thermal 
vibration  noise  in  the  upper  modes  of  the  cantilever,  in 
combination  with  the  large  gradient,  creates  enough 
magnetic  noise  to  destabilize  the  spin  [4-6]. 


RF  Field 
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Figure  1:  One  possible  configuration  for  a  magnetic  reso¬ 
nance  force  microscopy  measurement  (MRFM)  setup.  The 
magnetized  tip  of  an  ultra-sensitive  cantilever  is  brought  in 
close  proximity  to  the  surface  of  a  spin  sample  in  the 
presence  of  an  RF  field.  Vibration  of  the  cantilever  in 
combination  with  the  RF  field  causes  the  spins  within  the 
resonant  slice  to  undergo  cyclic  adiabatic  reversals.  The 
motion  of  the  spins  causes  a  slight  shift  in  the  cantilever’s 
natural  frequency.  The  coherence  of  the  derived  spin  signal 
is  affected  by,  amongst  other  factors,  thermal  noise  in  the 
cantilever  itself. 


To  overcome  this  problem,  the  MHz-frequency  thermo¬ 
mechanical  noise  of  the  cantilever  must  be  suppressed.  For 
detection  of  electron  spins,  the  frequency  range  that  is  most 
important  is  in  the  range  of  1 0  to  20  MHz,  corresponding  to 
the  Rabi  frequency  of  the  spin.  In  this  range,  vibration  noise 
of  only  a  few  milli-angstroms  can  be  detrimental. 

We  have  explored  several  approaches  for  suppressing 
the  high-order  mode  noise  [7].  This  work  focuses  on  mass- 
loaded  cantilevers  in  particular  (Fig.  2).  Finite-element 
analysis  shows  these  devices  to  have  much  lower  noise 
spectral  densities  at  MHz  frequencies.  This  design  can  be 
engineered  to  have  large  gaps  in  the  mode  spectrum  (Fig.  3). 


base 


Figure  2:  Schematic  of  a  mass-loaded  cantilever  (not  to 
scale).  In  this  design,  a  large  mass  is  suspended  at  the  end 
of  thin  hinge  which  is  connected  to  a  stiff,  thick  base. 
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Figure  3:  Predicted  thermal  noise  spectra  based  on  finite 
element  analysis  for  (a)  conventional  “flat”  cantilever,  and 
(b)  mass-loaded  cantilever.  Note  the  large  frequency  gaps 
in  graph  (b). 


DEVICE  FABRICATION 

Fabricating  cantilevers  for  the  present  application  poses 
special  challenges.  Unlike  previous  “flat”  cantilevers  [2], 
the  present  design  calls  for  a  dual-thickness  structure,  with  a 
very  thin  (0.1  pm)  “hinge”  section  coupled  to  a  considerably 
thicker  (2  pm)  “mass”.  To  minimize  dissipation,  the  canti¬ 
lever  has  to  be  made  of  single-crystal  silicon,  with  surfaces 


that  are  as  clean  and  defect-free  as  possible.  The  thickness 
of  each  section  needs  to  be  highly  uniform  as  well.  To  mini¬ 
mize  RF-induced  and  laser-induced  self-heating,  the  silicon 
material  must  have  very  low  electrical  conductivity,  i.e.  it 
must  be  undoped  or  very  lightly  doped.  Finally,  to  reduce 
clamping  losses,  the  overhang  at  the  base  of  the  cantilever 
(due  to  backside  misalignment)  must  be  minimized,  while 
the  base  itself  should  preferably  be  thickened  and  stiffened. 

To  address  these  requirements,  we  have  chosen  to  make 
the  cantilever  hinge  out  of  the  undoped  top  silicon  layer  of  a 
custom-made  silicon-on-insulator  (SOI)  wafer,  followed  by 
selective  undoped  silicon  epitaxy  to  form  the  mass. 

The  starting  material  is  an  SOI  wafer  with  a  0.1  pm 
thick  top  silicon  layer.  A  low-temperature  oxide  (LTO) 
layer  is  deposited  and  patterned  to  form  a  mask.  This  dark- 
field  mask  defines  the  thin  part  of  the  cantilever — the  hinge 
(Fig.  4a).  A  2  pm  thick  undoped  single-crystal  silicon 
epitaxy  is  grown.  The  epitaxy  is  selective,  i.e.  it  does  not 
grow  over  oxide.  The  oxide  mask  is  then  removed  with  HF 
(Fig.  4b). 

Another  LTO  layer  is  subsequently  deposited  and  pat¬ 
terned  to  form  a  second  mask.  This  mask  exposes  mainly 
the  base  of  the  cantilever  (Fig.  4c).  A  3  pm  thick  single¬ 
crystal  selective  silicon  epitaxy  is  then  grown,  and  the  LTO 
is  removed  with  HF  (Fig.  4d).  At  this  point  the  thickness  of 
the  base  is  effectively  5  pm,  providing  the  structural  rigidity 
required  for  reducing  clamping  losses  in  the  cantilever.  A 
lithographic  step  is  used  to  pattern  the  cantilever  and  its 
base,  and  a  silicon  plasma  etch  is  used  to  define  the  lever 
(Fig.  4e).  Note  that  this  etch  has  to  be  very  selective  against 
oxide,  since  the  thin  regions  of  the  top  silicon  layer  (where 
there  is  no  epitaxy)  etches  through  early  on  during  the  etch 
and  exposes  the  buried  oxide  to  the  plasma. 

Finally,  backside  lithography  followed  by  deep  reactive 
ion  etch  (DRIE)  is  used  to  create  through  holes  in  the  wafer. 
An  HF  etch  is  used  to  remove  the  buried  oxide,  and  the 
cantilevers  are  released  (Fig.  4f).  In  order  to  improve  yield, 
a  temporary  nitride-LTO  protective  layer  may  be  deposited 
on  the  front  side  of  the  wafer  before  the  backside  etch,  and 
the  nitride  removed  by  plasma  etching  before  the  HF  release. 
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Figure  3:  Fabrication  process  for  mass-loaded  cantilevers 
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Figure  5:  Mass-loaded  cantilevers:  (a)  long-hinge  design, 
(b)  short-hinge  design,  and  (c)  alternate  design  fabricated 
using  a  LOCOS-like  process.  Devices  similar  to  design  (a) 
were  used  in  the  majority  of  measurements  here. 


Table  1:  Predicted  and  measured  modes  for  long-hinge 
cantilever  (hinge  dimensions  72  x  3  x  0.1  pm;  mass  dimen¬ 
sions  15  x  3  x  2.0  pm).  Calculated  stiffness:  0.3  mN/m. 
*  torsional  modes  not  readily  detectable] 


Measured 

Predicted 

Type 

Fund,  mode  (kHz) 

6.3 

7.0 

Flexural 

2nd  mode  (kHz) 

94 

103 

Flexural 

3rd  mode  (kHz) 

* 

186 

Torsional 

4th  mode  (kHz) 

265 

323 

Flexural 

5th  mode  (kHz) 

* 

343 

Torsional 

Figure  6:  Thermal  mode  spectrum  of  long-hinge  cantilever 
measured  at  77  K  in  vacuum.  Numbers  in  parentheses  are 
theoretical  values  obtained  from  finite  element  analysis. 

Two  types  of  cantilevers  were  fabricated — “long-hinge” 
ones  (Fig.  5a)  and  “short-hinge”  ones  (Fig.  5b).  Both  types 
of  cantilevers  have  large  gaps  in  their  mode  spectra,  but  in 
practice  the  “long-hinge”  design  was  used  for  all  subsequent 
experiments  because  the  “short-hinge”  design  proved  to  be 
extremely  fragile  even  during  normal  handling.  While  both 
types  of  devices  are  sensitive  to  environmental  shock,  the 
high  stress  concentration  in  the  hinge  region  of  the  “short- 
hinge”  cantilever  makes  it  much  more  likely  to  break. 

In  an  alternate  fabrication  approach  (not  illustrated 
here),  an  SOI  wafer  with  a  thick  top  silicon  layer  is  used. 
The  hinge  regions  are  thinned  down  using  nitride-masked 
locdX  oxidation  of  silicon  much  like  the  VLSI  LOCOS 
process.  This  is  a  viable  option  if  silicon  epitaxy  capability 
is  not  available.  Possible  drawbacks  include  potentially 
large  variations  in  the  final  hinge  thickness  (corresponding 
to  thickness  variations  in  the  initial  top  silicon  layer),  and  the 
need  for  a  long,  precisely  timed  thermal  oxidation  step. 
Note  the  gradual  transition  in  thickness  between  the  hinge 
and  the  mass  (Fig.  5c),  characteristic  of  the  LOCOS  process. 


Figure  7:  Quality  factor  Q  of  a  long-hinge  mass-loaded 
cantilever  versus  temperature,  measured  in  vacuum.  The 
cantilever  has  a  magnetic  particle  glued  to  its  tip.  Q  at  room 
temperature  is  approximately  15000. 
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DEVICE  CHARACTERIZATION 

For  characterization  purposes,  a  long-hinge  cantilever 
was  mounted  in  a  setup  that  includes  a  fiber-based  laser 
interferometer.  The  thermal  mode  spectrum  of  the  cantilever 
was  measured  and  plotted  in  Fig.  6  (see  also  Table  1).  It  can 
be  seen  that,  as  expected,  there  is  a  large  gap  between  the 
fundamental  mode  and  the  higher  modes.  Theoretically 
obtained  modal  frequencies  (calculated  through  finite  ele¬ 
ment  analysis)  are  seen  to  be  in  good  agreement  with  the 
measured  values. 

The  quality  factor  Q  of  a  long-hinge  cantilever  was 
measured  in  vacuum  as  a  function  of  temperature.  Results 
are  shown  in  Figure  7.  It  can  be  seen  that  Q  can  reach 
80000,  and  that  it  varies  with  temperature.  In  general,  at  low 
temperatures  the  effect  of  surface  contaminants  are  “frozen 
out”  and  their  contribution  to  surface  losses  decreases, 
leading  to  higher  Q. 

In  the  experimental  configuration  used  for  our  MRFM 
measurements  (see  Fig.  1),  a  small  magnetic  particle  is 
attached  to  the  tip  of  the  cantilever  which  is  placed  in  close 
proximity  (approximately  0.1  pm)  to  an  electron-spin 
sample.  A  modulated  RF  field  is  applied  to  manipulate  the 


RF  power  (dBm) 


Figure  8:  Spin  signal  vs.  RF  power.  Varying  the  RF  power 
changes  the  associated  magnetic  field  strength,  which  in 
turn  changes  the  spin’s  Rabi  frequency.  The  suppression  of 
the  signal  at  3  dBm  may  be  due  to  thermal  noise  in  a 
cantilever  mode  that  coincides  with  the  Rabi  frequency  at 
that  point  (around  6  MHz). 


electron  spins,  and  the  motion  of  the  cantilever  resulting 
from  the  tip-sample  interaction  is  detected.  Because  the 
cantilever  is  very  soft  (0.3  mN/m),  attonewton- scale  forces 
can  be  detected  using  this  setup  with  the  help  of  a  high- 
performance  laser  interferometer. 

As  an  illustration  of  how  cantilever  modes  can  influence 
spin-detection  sensitivity,  MRFM  measurements  were  taken 
at  various  RF  power  levels  (Fig.  8).  Varying  the  RF  power 
changes  the  associated  magnetic  field  strength,  which  in  turn 
changes  the  spin’s  Rabi  frequency.  The  suppression  of  the 
spin  signal  at  3  dBm  may  be  due  to  thermal  noise  in  a 
cantilever  mode  that  coincides  with  the  Rabi  frequency  at 
that  point. 

Very  promising  results  have  been  obtained  so  far  with 
the  mass-loaded  cantilevers  in  our  MRFM  measurements. 
At  the  time  of  writing,  spin  signals  with  a  noise  floor  of  10 
spins  have  been  obtained.  Latest  results  will  be  reported  in 
detail  separately  [8]. 

This  work  was  supported  by  the  DARPA  Mosaic 
program.  The  cantilevers  were  fabricated  at  the  NSF 
National  Nanofabrication  User  Network  facility  at  Stanford 
University. 
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We  report  the  detection  of  the  v/V  statistical  polarization  in  a  small  ensemble  of  electron  spin  centers 
in  Si02  by  magnetic  resonance  force  microscopy.  A  novel  detection  technique  was  employed  that 
captures  the  statistical  polarization  and  cycles  it  between  states  that  are  either  locked  or  antilocked  to 
the  effective  field  in  the  rotating  frame.  Using  field  gradients  as  high  as  5  G/nm,  we  achieved  a 
detection  sensitivity  equivalent  to  roughly  two  electron  spins,  and  observed  ultralong  spin-lock 
lifetimes,  as  long  as  20  s.  Given  a  sufficient  signal-to-noise  ratio,  this  scheme  should  be  extendable 
to  single  electron  spin  detection. 
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In  his  classic  paper  on  nuclear  induction,  Bloch 
pointed  out  that  a  system  of  N  magnetic  moments  /. i 
will  give  rise  to  a  statistical  polarization  of  order  y[N / 1 
[1].  For  the  large  spin  ensembles  typically  used  in  con¬ 
ventional  magnetic  resonance  experiments,  this  statisti¬ 
cal  polarization  is  negligible  compared  to  the  thermal 
equilibrium  (Boltzmann)  polarization.  However,  for 
sufficiently  small  spin  ensembles  (especially  as  N  ap¬ 
proaches  unity),  the  statistical  polarization  can  exceed 
the  Boltzmann  polarization  and  even  dominate.  Given  a 
sufficiently  sensitive  means  to  detect  magnetic  resonance, 
it  should  be  possible  to  make  use  of  the  self-polarizing 
nature  of  small  spin  ensembles  to  perform  magnetic 
resonance  experiments,  just  as  one  manipulates  and  mea¬ 
sures  conventionally  polarized  ensembles  [2,3]-  This  ap¬ 
proach  offers  one  route  for  performing  magnetic 
resonance  experiments  on  small  numbers  of  spins,  per¬ 
haps  even  a  single  spin. 

In  the  present  work,  we  exploit  the  exquisite  sensitivity 
of  magnetic  resonance  force  microscopy  (MRFM)  [4-8] 
and  demonstrate  a  detection  sensitivity  equivalent  to  two 
spins  in  a  0.1  Hz  measurement  bandwidth.  This  represents 
nearly  2  orders  of  magnitude  improvement  over  previous 
MRFM  results  [8].  We  apply  the  technique  to  spin  en¬ 
sembles  consisting  of  fewer  than  102  electron  spins  and 
study  statistical  polarizations  of  order  10  spins.  Rather 
than  detecting  the  free  precession  of  the  transverse  mag¬ 
netization,  which  was  demonstrated  previously  for  very 
large  (N  ~  1023)  ensembles  of  nuclear  spins  [9-13],  we 
detect  the  longitudinal  component  in  the  “rotating 
frame”  using  the  technique  of  adiabatic  rapid  passage. 
By  applying  a  microwave  field  that  is  initially  off- 
resonance,  and  then  bringing  the  spins  into  resonance, 
the  spins  stay  “spin-locked”  to  the  effective  field  for 
times  of  order  Tlp,  the  spin-lattice  relaxation  time  in 
the  rotating  frame  [14,15].  Since  Tip  is  orders  of  magni¬ 
tude  longer  than  the  usual  spin  decoherence  time  T2,  the 
statistical  spin  signals  we  detect  exhibit  much  greater 
coherence  than  signals  from  the  free  precession  of  the 
transverse  spin  component. 


As  shown  in  Fig.  1,  our  MRFM  experiment  is  based  on 
an  ultrasensitive  cantilever  mounted  perpendicular  to  the 
sample.  At  the  end  of  the  cantilever  is  a  micron-size 
SmCo  magnetic  particle  that  generates  a  strong  magnetic 
field  gradient  (^2  G/nm).  A  microwave  field  Bi  ~  3  G 
from  a  superconducting  resonator  [16]  is  applied  to  excite 
electron  spin  resonance.  The  inhomogeneity  of  the  tip 
field  confines  the  magnetic  resonance  to  the  region  that 
satisfies  the  condition  B0(x,y,z )  =  (o/y,  where  co  is  the 
frequency  of  the  microwave  field,  y  is  the  gyromagnetic 
ratio  (y/27r  =  2.8  X  106  Hz/G),  and  B0  is  the  tip  field 


(b) 


magnetic  tip 

0I0-^ 

silicon  mass 

i - 1 

10  pm 

FIG.  1.  (a)  Schematic  of  the  MRFM  experiment.  As  the 

cantilever  oscillates,  the  spins  near  the  resonant  slice  are 
cyclically  inverted,  resulting  in  a  small  cantilever  frequency 
shift.  Cantilever  position  is  detected  using  an  optical  fiber 
interferometer  (not  shown),  (b)  Scanning  electron  micrograph 
of  the  single-crystal  silicon  cantilever  with  the  larger  (micron- 
size)  magnetic  tip.  The  cantilever  had  a  mass-loaded  geometry 
in  order  to  suppress  motion  of  the  tip  for  higher  order  flexural 
modes. 
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plus  an  optional  external  field.  In  our  case,  microwaves 
were  applied  at  2.96  GHz,  so  the  resonance  condition  was 
Bq(x,  y,  z)  =  1060  G.  The  entire  microscope  was  operated 
at  200  mK  in  a  small  vacuum  chamber  attached  to  the 
bottom  of  the  mixing  chamber  of  a  dilution  refrigerator. 

The  two  single-crystal  silicon  cantilevers  used  for  this 
study  had  a  mass-loaded  design,  consisting  of  a  roughly 
100  nm  thick,  70  jmm  long  shaft,  with  a  2  jmm  X  15  jmm 
silicon  mass  at  the  end  [Fig-  1(b)].  The  purpose  of  the 
mass  loading  was  to  suppress  the  thermal  motion  at  the 
end  for  the  higher  order  cantilever  modes  [17].  These 
modes  can  cause  unwanted  spin  relaxation  due  to  vibra- 
tionally  induced  magnetic  noise  near  the  Rabi  frequency 
(<z>Rabi  =  yB l)  [18-20].  One  cantilever  had  a  ~1  jmm 
magnetic  tip  attached  to  its  end  to  generate  the  field 
gradient.  In  order  to  increase  the  field  gradient,  the  second 
cantilever  had  a  magnetic  particle  that  was  shaped  using  a 
focused  ion  beam  so  that  the  tip  tapered  down  to  a  250  nm 
wide  apex.  The  cantilevers  had  spring  constants  k  of  about 
6  X  10-4  N/m,  frequencies  fc  of  6600  and  8600  Hz, 
respectively,  and  Q*  s  of  order  50  000  for  temperatures 
below  4  K. 

Two  samples  of  optically  polished  vitreous  silica  were 
studied.  Both  samples  were  irradiated  by  60Co  gamma 
rays  to  produce  silicon  dangling  bonds  known  as  E' 
centers  [15,21].  The  high  spin  density  sample  (Corning 
7943)  had  a  spin  density  of  ~1018  cm-3  and  was  studied 
with  the  larger  tip.  The  sample  with  lower  spin  density 
(Suprasil  W2)  had  a  spin  density  of  ~1015  cm-3  and  was 
studied  with  the  smaller  tip. 

To  detect  the  electron  spins,  we  look  for  a  shift  in 
cantilever  resonance  frequency  using  the  OSCAR  (oscil¬ 
lating  cantilever-driven  adiabatic  reversal)  protocol  [8]. 
The  cantilever  is  self-oscillated  at  its  resonance  fre¬ 
quency  through  the  use  of  a  gain-controlled  positive 
feedback  loop  [22].  As  the  cantilever  position  oscillates 
sinusoidally  according  to  xc(t)  =  vpkcos(<^),  the  field  B0 
at  a  given  sample  location  is  modulated  because  of  the 
field  gradient  from  the  tip  G  =  dB0/dx.  In  the  rotating 
frame,  the  effective  field  Beff  can  be  written  as  [14] 
BeffW  =  Bxx  +  [B0(t)  (o/y\ z  =  B{x  +  Gxc(t)  z, 

where  we  have  assumed  that  the  resonance  condition  is 
fulfilled  for  =  0.  If  Beff  changes  sufficiently  slowly,  the 
spins  will  be  “spin  locked”  (or  antilocked)  to  Beff .  Under 
these  circumstances,  the  z  component  of  magnetization 
will  therefore  oscillate  synchronously  with  the  cantilever 
position  as 


frequency  demodulator  [22].  In  the  limit  Gvpk  »  Bx,  the 
frequency  shift  Sfc  is  calculated  to  be  [23] 

<5 fc  =  — (2/cG/  7r&xpk)meff.  (2) 

Note  that  Sfc  is  directly  proportional  to  meff,  the  spin 
moment  aligned  with  the  effective  field  in  the  rotating 
frame.  In  the  present  experiment,  meff  is  the  result  of 
statistical  polarization  and,  over  time,  will  fluctuate  and 
even  reverse  sign  with  correlation  time  on  the  order  of 
Tip  [24]. 

The  above  description  applies  to  a  localized  region  of 
the  resonant  slice.  The  total  frequency  shift  will  be  a  sum 
of  contributions  from  the  entire  slice.  Since  the  spins  that 
are  in  different  regions  of  the  slice  experience  various 
values  of  G,  they  will  contribute  with  differing  weights  to 
the  overall  frequency  shift.  In  particular,  because  of  the 
perpendicular  cantilever  geometry,  spins  near  the  sur¬ 
face,  off  to  either  side  of  the  cantilever,  contribute  with 
the  strongest  weight,  while  spins  directly  below  the  tip, 
where  G  =  0,  contribute  minimally.  Because  of  symme¬ 
try,  the  measurement  responds  only  to  the  left-right  sta¬ 
tistical  imbalance  of  meff. 

To  improve  detectability  in  the  presence  of  low  fre¬ 
quency  noise,  and  to  endow  the  spin  signal  with  a  dis¬ 
tinctive  signature,  we  cycle  the  spins  between  locked  and 
antilocked  using  a  variation  of  the  protocol  we  refer  to 
as  “interrupted  OSCAR.”  The  microwave  field  B{  is 
normally  on,  but  it  is  then  interrupted  periodically  for 
one-half  of  a  cantilever  cycle,  starting  at  a  cantilever 
extremum,  as  shown  in  Fig.  2  (curve  B ).  During  the 
time  the  microwaves  are  off,  mz  remains  essentially 
static.  Since  the  cantilever  continues  to  oscillate,  when 


—  I — — - "*eii> 

V»i  +  (Gxc(t))2 

where  meff  is  the  component  of  the  magnetization  along 
Beff  and  the  sign  depends  on  whether  meff  is  aligned  or 
antialigned  with  Beff  (i.e.,  locked  or  antilocked). 

Because  of  the  back-action  force  on  the  magnetic  tip 
from  the  spins,  the  oscillating  mz(t )  results  in  a  frequency 
shift  of  the  cantilever,  which  is  detected  by  an  analog 


FIG.  2.  Timing  diagram  for  the  interrupted  OSCAR  protocol. 
The  cantilever  is  oscillated  continuously.  The  microwave  field 
(curve  B )  is  normally  on,  but  is  periodically  interrupted  for 
one-half  cantilever  cycle.  The  z  component  of  the  magnetiza¬ 
tion  (C)  oscillates  in  response  to  the  cantilever  motion  due  to 
adiabatic  rapid  passage  when  the  microwaves  are  on,  but  is  left 
static  when  they  are  off.  The  oscillating  magnetization  reverses 
phase  with  respect  to  the  cantilever  for  each  microwave  inter¬ 
ruption,  giving  a  cantilever  frequency  shift  ( D )  that  oscillates 
at  one-half  the  microwave  interrupt  frequency. 
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the  microwaves  are  turned  back  on  after  the  half-cycle 
gap,  Beff  will  have  reversed  orientation,  and  the  magne¬ 
tization  will  have  changed  from  locked  to  antilocked. 
This  change  from  locked  to  antilocked  and  vice  versa 
occurs  every  interrupt  cycle,  resulting  in  a  cantilever 
frequency  shift  (curve  D)  that  oscillates  at  one-half  the 
microwave  interrupt  frequency  /int.  The  fact  that  the 
signal  is  at  a  subharmonic  of  /int  gives  it  a  very  distinctive 
signature  that  is  free  of  spurious  feedthrough  artifacts. 

Although  somewhat  subtle  in  concept,  the  interrupted 
OSCAR  scheme  is  actually  quite  simple  in  practice.  Once 
the  cantilever  is  self-oscillating  and  pulses  with  correct 
timing  are  used  to  gate  the  microwaves,  one  can  simply 
look  for  a  peak  at  /int/2  in  the  power  spectrum  of  the 
frequency  demodulated  signal  D.  Alternatively,  phase 
sensitive  detection  via  a  lock-in  amplifier  can  be  em¬ 
ployed  to  detect  both  the  magnitude  and  the  sense  of 
the  net  polarization  meff.  In  either  case,  this  scheme  has 
the  appealing  feature  that  there  is  no  need  to  wait  a  spin- 
lattice  relaxation  time  Tx  between  measurements  for  the 
sample  to  repolarize,  since  we  do  not  rely  on  the 
Boltzmann  polarization.  This  aspect  is  potentially  quite 
useful  in  increasing  the  efficiency  of  taking  data  in  sys¬ 
tems  with  long  Tx  times  [3,10]. 

Figure  3(a)  shows  a  power  spectrum  of  the  frequency 
demodulated  cantilever  signal  obtained  for  the  case  of  the 
low  density  sample.  When  the  microwaves  were  inter- 


FIG.  3.  (a)  Power  spectral  density  of  the  frequency  demodu¬ 

lated  signal  for  the  low  spin  density  sample.  The  peak  at 
49.5  Hz  is  the  statistical  spin  signal.  The  integrated  signal  in 
the  peak  is  equivalent  to  that  from  about  6  spins,  and  the 
baseline  noise  is  1.8  spins  in  the  0.12  Hz  natural  bandwidth. 
The  cantilever  oscillation  amplitude  was  10  nm,  and  the 
resonant  slice  extended  roughly  100  nm  below  the  surface, 
(b)  Similar  spectra  from  the  high  spin  density  sample  at 
various  tip-sample  spacings.  As  the  resonant  slice  was  pulled 
out  of  the  sample,  the  spin  signal  decreased  monotonically. 


rupted  at  99  Hz,  a  prominent  peak  was  seen  at  49.5  Hz, 
as  expected.  The  peak  disappeared  when  the  microwave 
power  was  turned  off,  when  the  microwaves  were  applied 
cw,  when  the  sample  was  retracted,  or  when  the  external 
field  was  adjusted  so  that  the  resonant  slice  was  no  longer 
in  the  sample. 

The  sharpness  of  this  spectral  feature  (A/sig  =  0. 12  Hz 
FWHM)  implies  that  the  statistical  polarization  has  an 
extraordinarily  long  correlation  time  rm,  nearly  3  s,  as 
given  by  rm  =  l/7rA/sig.  The  time  rm  depended  on  the 
cantilever  oscillation  amplitude,  among  other  parame¬ 
ters,  and  was  observed  to  be  as  long  as  20  s  for  the  case 
of  60  nm  amplitude.  [The  larger  the  amplitude,  the  more 
complete  the  adiabatic  passages,  as  seen  by  Eq.  (1),  and 
the  more  time  the  spins  are  off-resonance.]  The  time  rm  is 
related  to  Tip  and  is  ultimately  limited  by  spin  relaxation 
caused  by  magnetic  noise  produced  by  the  tip  [8,18,20]. 

Figure  3(b)  shows  a  sequence  of  power  spectra  taken 
with  the  high  spin  density  sample  as  the  tip  was  retracted 
using  a  piezoelectric  actuator.  Each  curve  was  taken  with 
the  spacing  increased  by  25  nm.  The  signal  decreased, 
corresponding  to  the  smaller  volume  of  the  resonant  slice 
within  the  sample,  until  finally  the  slice  was  pulled  out  of 
the  sample  completely. 

By  using  synchronous  detection  to  detect  the  inter¬ 
rupted  OSCAR  signal,  we  could  monitor  fluctuations  of 
the  statistical  polarization  in  real  time.  Figure  4(a)  shows 
the  spin  fluctuations  for  the  low  density  sample.  Note  that 
significant  deviations  from  the  mean  can  persist  for  sev¬ 
eral  or  even  tens  of  seconds.  When  the  microwaves  were 
applied  cw,  only  the  background  noise  remained,  giving 
the  background  signal  shown  in  Fig.  4(b). 


Time  (s) 


FIG.  4.  (a)  Output  of  lock-in  amplifier  showing  real-time 

spin  fluctuations  for  the  low  spin  density  sample.  The  net 
polarization  wanders  between  being  in-phase  or  antiphase 
with  the  lock-antilock  cycling.  The  correlation  time  is  on  the 
order  10  s.  The  lock-in  filter  time  constant  was  3  s.  (b)  Lock-in 
amplifier  output  showing  the  baseline  noise  level  equivalent  to 
2.1  spins  rms. 
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If  we  know  the  field  gradient,  we  can  deduce  the 
number  of  spins  responsible  for  the  signal  using  Eq.  (2). 
We  do  not  have  the  ability  to  measure  G  =  3B0/dx  di¬ 
rectly,  but  we  have  measured  the  related  derivative 
dBz/ dz,  which  can  be  used  to  estimate  G  given  a  model 
for  the  tip.  To  determine  dBz/ dz,  we  measured  the  lock-in 
signal  as  a  function  of  external  field,  for  two  different 
values  of  tip-to-sample  spacing  [8,25].  For  the  smaller 
tip,  we  found  that  retracting  the  tip  by  55  nm  was  equiva¬ 
lent  to  offsetting  the  field  by  300  G.  Thus,  the  axial  field 
gradient  dBz/dz  in  the  center  of  that  field  range  is  given 
by  A5/Az  =  5.5  G/nm.  The  tip  was  then  modeled  sim¬ 
ply  as  a  uniformly  magnetized  spherical  tip  0.4  jmm  in 
diameter,  a  value  chosen  to  fit  the  measured  dBz/ dz .  The 
resulting  value  for  G  off  to  the  side  of  the  tip  at  the  sample 
surface  was  then  calculated  to  be  4.3  G/nm. 

Assuming  this  field  gradient,  we  can  use  Eq.  (2)  to 
assign  units  of  Bohr  magnetons  (jmB)  to  the  curves  in 
Fig.  3(a)  and  4.  The  integrated  peak  in  Fig.  3(a)  corre¬ 
sponds  to  5.5julb  rms.  The  baseline  noise  level  in  the 
signal  bandwidth  of  0.12  Hz  corresponds  to  1.8 julb  rms. 
In  Fig.  4(a),  the  spin  signal  fluctuations  correspond  to 
6.3 lx B  rms,  and  the  background  noise  level  in  Fig.  4(b)  is 
2.1  julb  rms  in  the  0.1  Hz  noise  bandwidth  of  the  lock-in 
amplifier. 

The  above  numbers  are  not  literally  the  number  of  net 
spins  responsible  for  the  signal.  Many  of  the  spins,  such 
as  those  directly  under  the  tip,  are  in  regions  where  G  is 
much  lower  than  the  assumed  value.  In  addition,  some  are 
not  optimally  placed  within  the  middle  of  the  resonant 
slice,  and  therefore  Eq.  (2)  overestimates  their  contribu¬ 
tion.  We  simply  point  out  that  the  signals  observed  are 
equivalent  in  strength  to  the  given  number  of  spins,  pro¬ 
vided  that  they  were  all  optimally  situated  at  the  stated 
field  gradient  and  undergoing  full  adiabatic  reversals. 

Unambiguous  detection  of  a  single,  isolated  spin  will 
require  a  sample  with  a  reduced  density  of  spin  centers 
(i.e.,  less  than  one  spin  on  average  within  the  resonant 
slice  volume).  Improved  signal-to-noise  ratio  (SNR)  will 
also  be  necessary  in  order  to  locate  the  spin  in  a  reason¬ 
able  amount  of  time.  Because  the  sign  of  the  spin  signal 
will  fluctuate,  signal  averaging  must  be  performed  on  a 
positive  definite  quantity  such  as  the  signal  power.  Such 
averaging  is  very  inefficient  when  the  single-shot  SNR  is 
low,  improving  only  as  n ^4,  where  n  is  the  number  of 
averages  [26,27].  To  reduce  the  rms  noise  from  2  spins  to 
0.2  spin,  for  example,  would  therefore  require  a  104  times 
increase  in  averaging  time.  Thus,  for  practical  single-spin 
detection,  more  fundamental  improvements  are  required, 
such  as  increasing  the  field  gradient  to  increase  the  mag¬ 
nitude  of  the  force  signal,  or  reducing  the  noise,  currently 
dominated  by  tip-surface  interactions  [28]. 
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Creating  Order  from  Random 
Fluctuations  in  Small 
Spin  Ensembles 

R.  Budakian,* *  H.  J.  Mamin,  B.  W.  Chui,  D.  Rugar 

We  demonstrate  the  ability  to  create  spin  order  by  using  a  magnetic  resonance 
force  microscope  to  harness  the  naturally  occurring  statistical  fluctuations  in 
small  ensembles  of  electron  spins.  In  one  method,  we  hyperpolarized  the  spin 
system  by  selectively  capturing  the  transient  spin  order  created  by  the 
statistical  fluctuations.  In  a  second  method,  we  took  a  more  active  approach 
and  rectified  the  spin  fluctuations  by  applying  real-time  feedback  to  the  entire 
spin  ensemble.  The  created  spin  order  can  be  stored  in  the  laboratory  frame  for 
a  period  on  the  order  of  the  longitudinal  relaxation  time  of  30  seconds  and  then 
read  out. 


Creating  order  from  random  thermal  fluctua¬ 
tions  has  been  of  interest  to  physicists  since  the 
development  of  statistical  mechanics  in  the 
19th  century  (1).  In  a  more  modern  context, 
creating  order  in  microscopic  physical 
systems  is  an  essential  part  of  quantum 
information  processing  and  quantum  compu¬ 
tation  (2-4),  where  the  ability  to  set  the  state 
of  a  collection  of  quantum  objects  to  a  desired 
configuration  is  required.  The  device  used  to 
perform  this  operation  must  be  capable  of 
controlling  the  microscopic  degrees  of  free¬ 
dom  of  the  system  while  being  subjected  to 
environmental  fluctuations. 

Here,  we  take  advantage  of  the  outstanding 
sensitivity  of  magnetic  resonance  force  mi¬ 
croscopy  (MRFM)  to  follow  statistical  y/N 
fluctuations  in  small  ensembles  of  electron 
spins  (5,  6)  with  a  real-time  sensitivity 
corresponding  to  1.3pB,  where  pB  is  the  Bohr 
magneton.  The  spin  manipulation  protocols 
we  have  developed  allow  us  to  monitor  and 
respond  to  the  instantaneous  spin  imbalance 
in  the  rotating  frame.  By  monitoring  the  spin 
system  and  selectively  capturing  the  large 
positive  fluctuations,  we  have  created  a  mean 
polarization  corresponding  to  ~6pB  in  an 
ensemble  of  A  ~  70  spins.  We  also  used  real¬ 
time  feedback  to  effectively  cool  the  spin 
system  and  create  a  mean  polarization 
corresponding  to  ~7pB.  The  spin  order  was 
then  transferred  to  the  laboratory  frame, 
stored,  and  later  read  out. 

In  MRFM  detection,  spins  are  manipu¬ 
lated  by  using  magnetic  resonance,  and  the 
longitudinal  component  of  the  magnetization 
is  detected  mechanically  by  measuring  the 
interaction  between  the  spins  and  a  small 
permanent  magnet  attached  to  the  end  of  a 
sensitive  silicon  cantilever.  Typically,  the 
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force  generated  by  a  spin  on  the  cantilever  is 
quite  small.  An  electron  spin  will  produce  a 
force  of  only  2  x  10-18  N  when  subjected  to  a 
magnetic  field  gradient  from  the  tip  as  large 
as  2  x  105  T  m_1.  Sensitivity  to  such  small 
forces  requires  the  ability  to  coherently 
manipulate  spins  for  many  cycles  of  the 
cantilever.  Recently,  through  the  use  of 
specially  engineered  cantilevers  that  reduce 
disturbance  to  the  spins  (7-9),  we  observed 
longitudinal  relaxation  times  in  the  rotating 
frame  (i.e.,  during  measurement)  of  up  to 
several  seconds.  This  has  allowed  us  to 
realize  a  single-shot  detection  sensitivity 
approaching  the  single  spin  level  (5)  and 
the  detection  of  an  isolated  electron  spin  by 
signal  averaging  (6). 

In  the  experimental  setup  of  the  MRFM 
apparatus  (Fig.  1A),  a  custom-fabricated 
mass-loaded  silicon  cantilever  with  a  sub¬ 
micrometer  SmCo  magnetic  particle  attached 
to  the  tip  is  used  as  the  force-sensing  element 
(5,  6).  The  sample  consists  of  vitreous  silica 
(Suprasil  W2,  Heraeus  Quarzglas  GmbH  and 
Company  KG,  Hanau,  Germany)  that  has 


been  irradiated  by  gamma  rays  from  a  60Co 
source  to  produce  spin- 1/2  paramagnetic 
defects  or  E'  centers  (unpaired  electron  spins 
on  Si)  (10).  Experiments  were  performed 
with  two  different  cantilever  and  sample 
combinations.  In  setup  1,  a  cantilever  having  a 
fundamental  resonance  frequency^  =  8.7  kHz 
and  stiffness  k  =  0.6  mN  m1  with  a  250-nm- 
wide  SmCo  tip  was  used  with  a  sample  that 
had  a  spin  concentration  of  ~1015  cm-3.  In 
setup  2,  a  cantilever  with  fc  =  5.5  kHz,  k  = 
0.11  mN  m-1,  and  a  1 50-nm-wide  SmCo  tip 
was  used  with  a  ~1014  cm-3  concentration 
sample.  The  MRFM  apparatus  was  operated 
in  vacuum  and  cooled  to  300  mK  to  reduce 
the  thermal  vibrations  of  the  cantilever. 

Electron  spin  resonance  was  excited  at 
corf/27r  =  2.96  GHz  with  the  use  of  a  micro- 
wave  field  with  amplitude  B]  ~  0.3  mT.  In  the 
presence  of  the  inhomogeneous  field  from  the 
tip,  only  those  spins  within  a  thin  resonant 
slice  satisfying  the  resonance  condition  given 
by  B0(x,y,z)  #  |B tip(x,y,z)  +  Bex tz|  =  corf/y 
will  interact  with  the  microwave  field.  Here, 
B0(x,  y,  z)  is  the  sum  of  the  tip  field,  B[  jp(x,  y,  z), 
and  a  uniform  external  field,  Bextz,  produced  by 
a  superconducting  magnet,  and  y  is  the  gy- 
romagnetic  ratio  (y/2 n  =  2.8  x  1010  Hz  T1). 
For  the  vertical  orientation  of  the  cantilever 
shown  in  Fig.  1A,  only  those  spins  that  are 
slightly  to  the  left  or  right  of  the  tip  contrib¬ 
ute  to  the  signal.  Furthermore,  because  of 
symmetry,  the  cantilever  will  respond  only  to 
the  left-right  imbalance  of  spin  polarization. 

To  detect  spins,  we  use  the  recently  de¬ 
veloped  spin  manipulation  protocol  OSCAR 
(oscillating  cantilever-driven  adiabatic  rever¬ 
sal),  which  measures  the  shift  in  the  funda¬ 
mental  frequency  of  the  cantilever  in  response 
to  tip-spin  interactions  (5,  6,  11).  The 
cantilever  is  self-oscillated  at  its  fundamental 
resonance  frequency  by  using  a  piezoelectric 
transducer  that  drives  the  cantilever  to  a 
fixed  amplitude  xpk  (11,  12).  As  the  cantile¬ 
ver  position  oscillates  according  to  xc(t)  = 
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Fig.  1.  (A)  Schematic  of  the  MRFM  apparatus.  The  cantilever  with  the  attached  magnetic  tip  is 
oriented  vertically  about  140  nm  away  from  the  surface  of  a  sample  containing  a  low  concentration 
of  f'centers.  A  small  coil  placed  near  the  sample  generates  microwaves  at  2.96  GHz  to  excite 
electron  spin  resonance.  At  2.96  GHz,  the  resonance  condition  for  electrons  is  met  for  B0  =  106  mT. 
For  Bext  =  30  mT,  the  resonant  slice  is  a  paraboloidal  shell  that  extends  about  250  nm  below  the  tip. 
(B)  Timing  diagram  for  the  interrupted  OSCAR  protocol.  For  setups  1  and  2,  the  microwave  (MW) 
power  is  interrupted  every  88  cycles  (/sig  ~  50  Hz)  and  64  cycles  (/sjg  ~  43  Hz),  respectively. 
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xpksin(27 ifct),  the  field  experienced  by  a  spin 
near  the  tip  is  modulated  sinusoidally  at  the 
cantilever  frequency  with  peak  amplitude 
A B  =  Gx  k,  where  G  =  dBJdx  is  the  lateral 
gradient  from  the  tip.  For  typical  experimen¬ 
tal  parameters  of  xpk  ~  20  nm  and  G  ~  2  x 
105  T  m_1,  the  external  field  experienced  by 
spins  within  the  resonant  slice  is  modulated 
by  AB  =  4  mT.  In  the  presence  of  the  micro- 
wave  field,  the  slow  variation  of  the  static 


field  at  the  cantilever  frequency  causes  each 
spin  within  the  resonant  slice  to  be  adia- 
batically  inverted  synchronously  with  the 
cantilever  motion.  The  back  action  of  the  spin 
on  the  cantilever  in  turn  shifts  the  fundamen¬ 
tal  frequency  of  the  cantilever  by  a  small 
amount 
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Fig.  2.  (A)  Trace  showing  a  time  record  of  the  statistical  fluctuations  recorded  in  a  83-mHz 
bandwidth  with  the  use  of  setup  1.  The  frequency  shift  is  converted  to  equivalent  number  of  spins 
(right-hand  axis)  by  dividing  Af  by  the  average  frequency  shift  per  spin  (\Sf  \  =  0.8  mHz/spin).  (B) 
Histogram  of  a  1-hour  continuous  record  of  the  in-phase  lock-in  output  signal.  The  Gaussian 
distributed  spin  fluctuations  have  a  zero  mean,  indicating  that  the  time-averaged  spin  imbalance 
does  not  point  in  a  preferred  direction  with  respect  to  the  effective  field.  This  observation  is 
consistent  with  the  fact  that  the  spin  temperature  in  the  rotating  frame  approaches  Ts  =  oo  for 
measurement  times  tmeas  »  xm.  (C)  Histogram  of  the  quadrature  channel  showing  the  detection 
noise.  (D)  Power  spectrum  of  a  1-hour  record  of  the  fluctuations. 


Fig.  3.  (A)  Trace  showing 
the  average  of  2800  indi¬ 
vidual  capture-store-readout 
sequences  with  a  storage  time 
of  1.0  s,  taken  with  use  of 
setups  The  average  frequency 
shift,  A f{t),  and  stored  magne¬ 
tization,  mstored,  are  shown  on 
the  left-  and  right-hand  axes, 
respectively.  Equivalent  spins 
are  calculated  by  using  \Sf  \  = 
1.7  mHz/spin.  (B)  Data  taken 
under  the  same  conditions  as 
(A),  except  the  storage  time 
was  increased  to  5.0  s.  The  de¬ 
caying  magnetization  observed 
during  readout  in  (A)  and  (B) 
fits  well  to  an  exponential  indi¬ 
cating  aim  =  5.0  s  spin  relaxa¬ 
tion  time  during  measurement. 


where  G.  is  the  gradient  at  the  position  of  the 
zth  spin  and  a.  is  a  random  variable  that  takes 
values  +1  or  -1,  depending  on  whether  the 
spin  moment  is  aligned  or  anti-aligned  with 
respect  to  the  effective  field  in  the  rotating 
frame  (5,  13 ,  14).  The  two-level  nature  of  ai 
reflects  the  quantum  measurement  character¬ 
istic  of  MRFM  spin  detection  (i.e.,  a  rotating 
frame  Stern-Gerlach  behavior)  (15-17).  For 
an  ensemble  of  N  spins,  we  define  the 

instantaneous  statistical  spin  imbalance  in 

N 

the  rotating  frame  to  be  AN  =  a*  and  the 

i=  1 

corresponding  magnetic  moment  to  be  Am  = 
A/VpB.  The  cantilever  frequency  shift  for  the 

N 

ensemble  is  given  by  Af  =  $fi- 

i=  1 

Near  the  surface  of  the  sample,  the 
cantilever  frequency  is  affected  not  only 
by  the  presence  of  the  spins  but  also  by 
the  more  dominant  electrostatic  and  van  der 
Waals  forces.  In  order  to  make  the  spin 
component  of  the  signal  distinctive,  we  pe¬ 
riodically  reversed  the  sign  of  the  frequency 
shift  by  interrupting  the  microwave  power 
for  1/2  cycle  of  the  cantilever  vibration  every 
Nint  cycles  (Fig.  IB).  The  periodic  inter¬ 
ruptions  encode  the  spin  signal  with  a  char¬ 
acteristic  frequency  fsi  =  fJ2Nint.  A  lock-in 
amplifier  referenced  to  /sig  is  used  to  de¬ 
modulate  the  frequency  shift  and  determine 

4 m  6). 

A  200-s  record  of  the  cantilever  frequen¬ 
cy  shift,  recorded  from  the  in-phase  channel 
of  the  lock-in  amplifier,  is  shown  (Fig.  2A). 
A  histogram  of  the  frequency  fluctuations 
(Fig.  2B)  reveals  a  Gaussian  shape  with  var¬ 
iance  o2  =  69  mHz2.  We  note  that  the  dis¬ 
tribution  has  a  zero  mean,  indicating  that 
the  spin  ensemble  does  not  have  a  preferred 
direction  with  respect  to  the  effective  field. 
To  show  that  the  measured  frequency  shift 
is  dominated  by  the  spin  signal  rather  than 
measurement  noise,  we  show  the  histogram 
of  the  measurement  noise  recorded  from 
the  quadrature  channel  of  the  lock-in  am¬ 
plifier  (oQ2  =  5  mHz2)  (Fig.  2C).  The  vari¬ 
ance  of  the  spin  signal  without  measurement 
noise  is  simply  related  to  the  in-phase  and 
quadrature  variances  by  ospin2  =  o^  -  oQ2 
(6).  From  the  measured  variance  ospi2  and 
a  model  of  the  field  produced  by  the  tip 
(18),  we  estimate  the  total  number  of  spins 
contained  within  the  resonant  slice  to  be 
N  ~  70. 

The  power  spectrum  of  the  frequency 
fluctuations  (Fig.  2D)  fits  well  to  a  Lorentzian 
with  spectral  width  Avs  =  21  mHz,  cor¬ 
responding  to  a  correlation  time  xm  = 
l/27iAvs  =  7.6  s.  This  long  correlation  time, 
closely  related  to  the  rotating  frame  longi¬ 
tudinal  relaxation  time,  allows  us  to  mea¬ 
sure  real-time  fluctuations  of  the  spin 
ensemble  with  a  lock-in  time  constant  of 
up  to  3  s. 
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Fig.  4.  (A)  A  200-s  record  of  spin  fluctuations,  taken  with  use  of  setup  1,  that  has 
been  rectified  by  using  active  feedback  to  maintain  positive  spin  polarization. 
The  7i  inversions  (indicated  by  the  vertical  bars)  are  applied  by  omitting  a  single 
1/2  cycle  interrupt  from  the  standard  iOSCAR  protocol.  (B)  Histogram  of  a  1- 
hour  record  of  the  rectified  signal  showing  a  positive  mean  frequency  shift  of  5.6 
mHz.  An  analysis  of  the  histogram  indicates  an  80%  probability  that  the  spin 
polarization  will  be  positive.  For  comparison,  in  the  absence  of  feedback  (dashed 
curve),  the  spin  polarization  is  within  4%  of  having  no  preferred  direction.  (C)  An 
average  polarization  of  about  3.7pB  was  created  in  an  ensemble  of  N  ~  70  spins 
through  active  feedback,  stored  in  the  laboratory  frame  for  2.5  s,  and 
subsequently  read  out.  The  data  represents  2500  averages  taken  with  use  of 
setup  2.  (D)  Data  taken  under  the  same  conditions  as  (C)  except  for  a  10-s 
storage  time. 
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As  an  initial  demonstration  of  harnessing 
spin  fluctuations,  we  created  a  hyperpolar- 
ized  spin  state  by  selectively  capturing  and 
storing  the  especially  large  statistical  fluctu¬ 
ations.  To  do  this,  we  continuously  monitor 
the  spin  signal  and  wait  until  a  fluctuation 
exceeds  a  predetermined  threshold  value 
threshold  (7p)-  uP°n  registering  a  suitable 
fluctuation,  the  microwave  power  is  turned 
off  at  a  maximum  of  the  cantilever  motion, 
which  leaves  the  instantaneous  spin  polar¬ 
ization  pointing  along  B0.  In  the  absence  of 
the  microwave  field,  the  spins  no  longer 
respond  to  the  cantilever  motion,  and  the 
nonequilibrium  state  of  the  spin  ensemble 
can  be  stored  in  the  laboratory  frame  for  as 
long  as  a  spin-lattice  relaxation  time,  Tv 
The  stored  magnetization  can  then  be  read 
out  by  reapplying  the  microwave  field  at  a 
maximum  of  the  cantilever  motion  and  using 
the  standard  interrupted  OSCAR  (iOSCAR) 
protocol. 

The  average  of  2800  individual  capture- 
store-readout  sequences,  taken  with  use  of 
setup  2,  is  shown  (Fig.  3).  For  the  data  in  Fig. 
3 A,  the  captured  magnetization  was  stored 
for  1  s  in  the  laboratory  frame  (i.e.,  with  the 
microwave  field  off)  and  then  read  out  with 
use  of  the  iOSCAR  protocol.  The  readout 
signal  had  a  peak  amplitude  of  9.4  mHz, 
which  we  estimate  corresponds  to  an  average 
magnetization  of  mstored  =  5.5jib.  When  the 
storage  time  was  increased  to  5  s  (Fig.  3B), 
the  peak  stored  magnetization  decreased  by 
14%  to  mst0red  =  4.8jib.  This  drop  in  mstored 
is  the  result  of  depolarization  due  to  longi¬ 


tudinal  relaxation,  indicating  Tx  ~  30  s.  The 
observed  Tx  includes  the  contribution  from 
the  lattice  as  well  as  tip-induced  relaxation 
(77). 

In  addition  to  simply  selecting  and  cap¬ 
turing  desired  fluctuations,  we  can  also  take 
a  more  active  approach  by  applying  real¬ 
time  feedback  to  the  spin  system  in  order  to 
continuously  guide  its  evolution.  As  a  demon¬ 
stration  of  feedback  control,  we  have  rectified 
the  spin  fluctuations  by  monitoring  the  spin 
signal  and  applying  a  n  inversion  to  the  entire 
spin  ensemble  whenever  Af  <  0  (19).  The  n 
inversions,  accomplished  with  the  use  of 
adiabatic  inversion,  flip  the  sign  of  the  spin 
imbalance  so  as  to  always  keep  Am  positive 
in  the  iOSCAR  reference  frame.  Figure  4A 
shows  a  200-s  record  of  the  iOSCAR  sig¬ 
nal,  taken  with  the  use  of  setup  1,  along 
with  vertical  bars  indicating  times  when  n 
inversions  were  applied.  In  contrast  to  Fig. 
2B,  the  histogram  of  the  signal  with  feedback 
control  (Fig.  4B,  solid  curve)  now  shows  a 
nonzero  mean  value  of  5.6  mHz  corre¬ 
sponding  to  ~7.0jib.  Thus,  through  the  use 
of  feedback,  we  have  essentially  hyper- 
polarized  the  spins  in  the  rotating  frame  of 
the  iOSCAR  measurement.  This  spin  order 
can  once  again  be  transferred  to  the  labora¬ 
tory  frame,  stored,  and  then  read  out  (Fig.  4, 
C  and  D). 

We  have  demonstrated  real-time  con¬ 
trol  of  electron  spins  in  small  ensembles 
using  two  spin  manipulation  protocols:  fluctu¬ 
ation  capture  and  fluctuation  rectification. 
Because  the  present  single-shot  detection 


sensitivity  is  already  approaching  the  sin¬ 
gle  spin  level,  relatively  modest  improve¬ 
ments  in  detection  signal-to-noise  ratio 
should  allow  real-time  quantum  state  de¬ 
tection  and  control  of  individual  electron 
spins. 
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Slip-Rate  Measurements  on  the 
Karakorum  Fault  May  Imply 
Secular  Variations  in  Fault  Motion 

M.-L.  Chevalier,1,2  F.  J.  Ryerson,2*  P.  Tapponnier,1  R.  C.  Finkel,2 
J.  Van  Der  Woerd,3  Li  Haibing,4  Liu  Qing5 

Beryllium-10  surface  exposure  dating  of  offset  moraines  on  one  branch  of  the 
Karakorum  Fault  west  of  the  Gar  basin  yields  a  long-term  (140-  to  20- 
thousand-year)  right-lateral  slip  rate  of  ~10.7  ±  0.7  millimeters  per  year.  This 
rate  is  10  times  larger  than  that  inferred  from  recent  InSAR  analyses  (~1  ±  3 
millimeters  per  year)  that  span  -8  years  and  sample  all  branches  of  the  fault. 
The  difference  in  slip-rate  determinations  suggests  that  large  rate  fluctuations 
may  exist  over  centennial  or  millennial  time  scales.  Such  fluctuations  would  be 
consistent  with  mechanical  coupling  between  the  seismogenic,  brittle-creep, 
and  ductile  shear  sections  of  faults  that  reach  deep  into  the  crust. 


The  Karakorum  Fault  in  Tibet  is  the  main 
Quaternary  right-lateral  fault  north  of  the 
Himalayas.  Determining  its  past  and  present 
motion  is  critical  to  understanding  the  kine¬ 
matics  of  Asian  continental  deformation  and 
the  rheology  of  the  continental  lithosphere 
(7, 2).  The  fault  trends  roughly  parallel  to  the 
western  Himalayan  range,  extending  from  at 
least  Kailas  to  the  Pamirs,  a  length  of  >1200 
km  (Fig.  1).  Its  Quaternary  slip  rate  remains 
poorly  constrained,  compared  to  that  of  other 
large  faults  in  Asia  such  as  Kunlun,  Haiyuan, 
and  Altyn  Tagh  (3-5).  Previous  attempts  to 
determine  the  rate  have  produced  disparate  val¬ 
ues  ranging  from  1  to  30  mm/year  (2,  6-11). 
Such  disparities  may  result  from  the  different 
techniques  applied  and  time  periods  observed, 
the  part  of  the  fault  investigated,  or  its  com¬ 
plex  geometry,  which  displays  multiple  splays 
with  oblique  slip  (12).  We  present  measure¬ 
ments  of  the  Mid-  to  Late  Pleistocene  slip 
rate  on  the  southern  stretch  of  the  fault, 
based  on  10Be  surface  exposure  dating  of  two 
moraine  crests  displaced  by  the  fault  at  the 
Manikala  glacial  valley  terminus  (32°2.529'N, 
80°1.212'E,  4365  to  4760  m  above  sea  level) 
(Figs.  1  and  2). 
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Globe  de  Strasbourg,  UMR  7516,  CNRS,  Strasbourg, 
France,  laboratory  of  Continental  Dynamics,  Institute 
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Beijing  100037,  China.  5Total  Exploration  China,  Total- 
Fina-Elf,  Beijing,  100004,  China. 
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The  Manikala  moraine  complex  lies  at  the 
base  of  the  faulted  Ayilari  range  front,  which 
bounds  the  west  side  of  the  Gar  valley,  a  large 
pull-apart  basin  floored  by  marshland  that 
hides  other  strands  of  the  Karakorum  Fault 
system  (13)  (Fig.  1).  The  moraines,  Ml  and 
M2,  lie  southeast  of  the  U-shaped  Manikala 
Valley,  a  glacial  trough  deeply  entrenched 
into  the  range’s  igneous  basement  (Fig.  2  and 
fig.  SI).  The  range  front  is  marked  by  tri¬ 
angular  facets  as  high  as  800  m  that  testify  to 
a  normal  component  of  slip  on  the  fault.  The 
principal  strand  of  the  fault  shows  discrete 
right-lateral  offsets  (10  +  2  m,  35±5m,  and 
75  ±  5  m)  of  young  rills  of  different  depths, 
incised  into  postglacial  colluvium  (fig.  SI). 
Within  the  till  complex,  two  main  groups  of 
moraines  are  recognized  (Fig.  2  and  fig.  S2). 
All  were  emplaced  by  the  Manikala  Daer  gla¬ 
cier,  whose  terminus  is  today  ~7  km  upstream. 
The  morphology  of  the  moraines  indicates 
that  they  correspond  to  major  advances  of 
the  glacier  and  were  later  abandoned  when 
the  glacier  retreated  upstream  (fig.  S3). 

The  relative  ages  of  the  moraine  groups 
can  be  qualitatively  assessed  from  their 
surface  characteristics  (Fig.  2  and  fig.  S2). 
The  Ml  surface  is  rough  and  composed  of 
chaotically  distributed,  imbricate  blocks  (as 
large  as  3  m  in  diameter)  surrounded  by 
coarse  debris.  The  smoother  surface  of  M2 
appears  older,  with  blocks  (tens  of  centi¬ 
meters  to  a  meter  in  diameter)  protruding 
above  a  mantle  of  smaller  debris  (Fig.  2  and 
fig.  S2).  The  morainic  ridges  thus  appear  to 
become  younger  from  east  to  west,  consist¬ 
ent  with  right-lateral  motion  on  the  fault. 

The  M2  moraine  complex  is  divided  into 
eastern  and  western  sections  (M2E  and 
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M2W,  respectively)  by  a  deep,  beheaded, 
flat-floored  paleovalley  that  is  truncated  by 
the  fault  and  flanked  by  well-defined  lateral 
moraines  (Fig.  2).  The  crest  of  the  lateral 
moraine  east  of  the  paleovalley  is  well  pre¬ 
served,  and  its  eastern  edge  extends  to  the 
base  of  the  faceted  range  front.  There  is  no 
catchment  on  the  mountain  slope  facing 
this  valley,  indicating  that  it  must  corre¬ 
spond  to  a  former  channel  of  the  Manikala 
glacier  (Fig.  2  and  fig.  S3).  The  youngest 
moraine  group,  Ml  (Fig.  2),  is  the  only  one 
present  on  both  sides  of  the  Manikala 
outwash  valley  and  displays  well-preserved 
terminal  lobes  and  sharply  defined  ridge 
crests. 

Upstream  from  the  fault,  the  limits  of 
glacial  incision  reach  the  base  of  the  trian¬ 
gular  facets  that  border  the  Manikala  valley. 
Downstream,  the  Ml  and  M2E  moraine  ridge 
crests  extend  linearly  to  the  fault  and,  when 
realigned  with  the  sharp  edge  of  glacial 
bedrock  incision  south  of  the  fault,  provide 
the  only  piercing  points  accurate  enough  to 
obtain  offset  estimates  (14)  (Fig.  2  and  fig. 
S4).  Once  restored  from  satellite  images,  the 
Ml  and  M2E  offsets  are  220  ±  10  m  and 
1520  ±  50  m,  respectively.  Another  moraine 
complex  with  morphologically  similar 
surfaces  and  offsets  is  found  at  the  terminus 
of  the  Tajiang  Daer  glacial  valley  ~10  km 
to  the  west  (fig.  S3),  lending  additional 
support  to  this  reconstruction  (fig.  S4). 

10Be  model  ages  for  samples  collected 
along  the  moraine  ridge  crests  (Fig.  2B) 
define  consistent  age  clusters  that  can  be 
used  to  date  their  abandonment  (15).  Ages 
on  the  two  M2  ridges  (Fig.  3)  fall  mostly 
between  103  and  204  thousand  years  (ky) 
(15  samples;  mean  age,  152  +  28  ky),  with  a 
subset  of  seven  samples  between  132  and 
150  ky  old  (mean  age,  140  ±  5.5  ky). 
Samples  WG-3,  WG-4,  and  WG-7,  on  the 
eastern  moraine  (M2E),  are  more  than  55  ky 
older  than  the  main  M2E  population  (132  to 
150  ky  old).  We  consider  these  three  samples 
to  be  outliers,  probably  originating  from 
more  ancient  till  upstream.  Examination  of 
the  M2  population  with  the  outliers  excluded 
reveals  two  distinct  subgroups:  a  younger 
cluster  of  nine  samples  with  ages  of  103  to 
149  ky  (mean  age,  133  ±  15  ky)  and  an  older 
group  of  six  samples  of  1 60  to  204  ky  (mean 
age,  181  ±  14  ky).  The  nine  samples  from  the 
Ml  moraine  yield  a  younger  mean  age  of 
35  ±  9  ky.  Seven  samples  fall  between  36 
and  45  ky  (mean  age,  40  +  3  ky),  and  two 
distinctly  younger  samples  have  ages  of 
21  +  1.0  ky. 
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Suppression  of  Spin  Diffusion  near  a  Micron-Size  Ferromagnet 

R.  Budakian,  H.  J.  Mamin,  and  D.  Rugar 

IBM  Research  Division,  Almaden  Research  Center,  650  Harry  Road,  San  Jose,  California  95120-6099,  USA 

(Received  1  July  2003;  published  23  January  2004) 

We  have  used  the  large  gradients  generated  near  the  ferromagnetic  tip  of  a  magnetic  resonance  force 
microscope  to  locally  suppress  spin  diffusion  in  a  silica  sample  containing  paramagnetic  electron  spins. 

By  controlling  the  slice  location  with  respect  to  the  tip,  the  magnetic  field  gradient  was  varied  from  0.01 
to  36  mT/yam,  resulting  in  a  fourfold  decrease  in  Tfl  and  a  similar  decrease  in  .  The  observed 
dependence  of  the  relaxation  rates  on  field  gradient  is  consistent  with  the  quenching  of  flip-flop 
interactions  that  mediate  the  transport  of  magnetization  between  slow  and  fast  relaxing  spins. 

DOI:  10. 1 103/PhysRevLett.92.037205  PACS  numbers:  75.40.Gb,  07.55.-w,  76.30.-v,  85.70.-w 


In  recent  years,  numerous  proposals  have  been  made 
for  solid-state  devices  that  operate  using  the  spin  degree 
of  freedom  of  electrons  or  nuclei  [1-4].  The  successful 
function  of  these  spintronic  and  quantum  computing  de¬ 
vices  depends  critically  on  minimizing  sources  of  spin 
relaxation  and  decoherence  from  the  environment  [5,6]. 
Long  relaxation  times  are  also  essential  for  the  detection 
and  manipulation  of  single  spins  by  magnetic  resonance 
force  microscopy  (MRFM)  [7]. 

Spin  diffusion  has  long  been  recognized  as  a  highly 
efficient  mechanism  for  mediating  spin  relaxation.  The 
concept  was  originally  proposed  by  Bloembergen  to  ac¬ 
count  for  the  anomalously  fast  nuclear  spin-lattice  relaxa¬ 
tion  observed  in  CaF2  and  other  ionic  insulating  crystals 
[8].  According  to  Bloembergen,  fluctuating  fields  gener¬ 
ated  by  paramagnetic  impurities  rapidly  relax  nearby 
nuclear  spins.  Dipolar  coupling  between  adjacent  nuclear 
spins  gives  rise  to  flip-flop  interactions  that  mediate  the 
transport  of  magnetization  from  the  bulk  to  the  paramag¬ 
netic  sinks,  greatly  reducing  the  overall  Tx.  Waugh  and 
his  colleagues  showed  that  the  effect  of  spin  diffusion 
could  be  reversed  and  the  nonequilibrium  state  of  slow 
relaxing  spins  could  be  recovered  by  a  suitable  sequence 
of  strong  rf  pulses  which  reverse  the  sign  of  the  dipolar 
Hamiltonian  [9]. 

A  more  direct  route  to  controlling  spin  diffusion 
is  to  prevent  flip-flop  processes  from  taking  place. 
The  dipolar  Hamiltonian  contains  the  term  Hss  = 
—  jm0y2h2(\  —  3cos 20)(SxSf  +  Sx  S^/^irr3  which  per¬ 
mits  two  like  spins  separated  by  r  to  exchange  magneti¬ 
zation  while  conserving  the  total  Zeeman  energy.  A 
sufficiently  intense  magnetic  field  gradient  could  be 
used  to  prevent  simultaneous  flips  from  occurring  if  the 
difference  in  Zeeman  energies  of  two  adjacent  spins  could 
be  made  larger  than  their  dipolar  coupling.  If  we  consider 
a  homogeneous  solid  like  CaF2,  for  example,  we  find  that 
the  gradient  required  to  suppress  flip-flops  between  two  F 
nuclei  is  G  ~  ji^yh/Airr 4  ^  4  X  102  mT/ jmm,  which  is 
orders  of  magnitude  larger  than  the  largest  gradients  used 
in  conventional  magnetic  resonance  experiments  [10]. 


While  such  large  gradients  are  difficult  to  generate  in 
the  laboratory,  there  are  instances  where  they  occur  natu¬ 
rally  on  the  atomic  scale  [8,11].  For  instance,  the  same 
paramagnetic  impurities  that  drive  macroscopic  spin  dif¬ 
fusion  in  solids  are  surrounded  by  microscopic  regions  or 
“frozen  cores”  where  nuclear  spin  diffusion  is  quenched 
due  to  the  large  field  gradient  from  the  impurity  [12,13]. 

Spin  measurement  using  MRFM  relies  on  large  field 
gradients  generated  near  a  micron-size  ferromagnet  (the 
MRFM  tip)  in  order  to  detect  small  ensembles  of  spins 
with  high  sensitivity  and  nanometer-scale  spatial  resolu¬ 
tion  [14-16].  In  the  present  work,  we  exploit  these  intense 
gradients  to  investigate  the  spin-lattice  relaxation  of  a 
dilute  concentration  of  paramagnetic  electron  spins  as  a 
function  of  the  field  gradient  near  an  MRFM  probe  tip.  We 
show  that  by  controlling  the  relative  tip  position,  we  can 
locally  suppress  spin  diffusion  and  thereby  affect  spin- 
lattice  relaxation  rates. 

The  MRFM  apparatus,  shown  schematically  in  Fig.  1, 
was  placed  in  an  evacuated  vessel  and  cooled  to  4.2  K 
The  custom  designed  single-crystal  silicon  cantilever 
with  dimensions  210  jmm  X  3.3  jmm  X  0.86  jmm  had  a 
fundamental  flexural  frequency  fc  =  5.8  kHz,  a  spring 
constant  k  =  0.02  N/m,  and  a  Q  of  90  000.  Displacement 
of  the  cantilever  was  detected  using  an  optical  fiber 


FIG.  1.  Schematic  of  MRFM  apparatus.  Two  slice  loca¬ 
tions  are  indicated  for  different  values  of  the  external 
field.  Effective  slice  thickness  (£)  is  calculated  for  £1/2tt  = 
20  MHz.  (a)  Bexi  =  80  mT,  dslice  =4.3  jam,  £  =  240  nm. 
(b)  Bext  =  50  mT,  dslice  =  2.5  jmm,  £  =  90  nm. 
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interferometer  coupled  to  a  1550  nm  distributed  feedback 
diode  laser.  The  magnetic  tip  was  fabricated  by  gluing  a 
small,  high  anisotropy  SmCo  particle  near  the  end  of  the 
cantilever  in  the  presence  of  an  external  orienting  field 
After  gluing,  the  bottom  portion  of  the  particle  was 
milled  flat  using  a  focused  ion  beam.  Magnetic  charac¬ 
terization  of  the  particle  was  carried  out  using  cantilever 
magnetometry  [17]  which  showed  that  the  shaped  ferro- 
magnet  with  dimensions  6  jmm  X  7  jmm  X  11  jmm  had  a 
total  moment  of  jmtip  =  2.3  X  10-10  J/T  and  was  oriented 
to  within  20°  of  the  cantilever  normal. 

The  sample  consisted  of  optically  polished  vitreous 
silica  (Corning  7943)  that  had  been  y  irradiated  by  a 
60Co  source  (150  Mrad)  to  produce  paramagnetic  defects 
or  E'  centers  (unpaired  electrons  on  Si)  with  a  concen¬ 
tration  of  ns  ~  2  X  1018  cm-3.  The  spectrum  of  the  E' 
centers  was  obtained  at  room  temperature  using  a  cw  ESR 
spectrometer  operating  at  9.3  GHz.  The  spectrum,  shown 
in  Fig.  2,  is  inhomogeneously  broadened  and  has  the 
familiar  powder  line  shape  characteristic  of  a  vitreous 
sample  [18].  Prior  to  the  MRFM  measurements,  the 
sample  surface  was  coated  with  a  100  A  thick  Ti/Au 
layer  to  minimize  unwanted  electrostatic  tip-surface 
interactions. 

Electron  spins  were  excited  at  corf/27r  =  2.96  GHz 
using  a  superconducting  microstripline  resonator  coupled 
to  a  200  jmm  diameter  coil  which  generated  a  microwave 
field  Bx  =  0.16  mT  [19].  The  spins  contained  within  the 
resonant  slice  satisfy  the  magnetic  resonance  condition 
given  by  <urf  =  yB0(R),  where  y  is  the  gyromagnetic 
ratio  (y/27r  =  2.8  X  1010  Hz/T)  and  2?0(R)  is  the  sum 
of  the  inhomogeneous  field  from  the  tip  Z?tip(R)  and  a 
uniform  external  field  Z?ext  produced  by  a  superconducting 
magnet.  At  2.96  GHz,  the  resonance  condition  for  elec¬ 
trons  is  met  when  B0  =  106  mT.  By  changing  Z?ext,  both 
the  location  of  the  resonant  slice  with  respect  to  the  tip 
(dslice)  as  well  as  the  strength  of  the  gradient  within  the 
slice  can  be  varied.  This  makes  it  possible  to  study  the 
dependence  of  T ,  on  the  strength  of  the  field  gradient. 


FIG.  2.  ESR  spectrum  of  silica  sample  containing  E'  centers. 
The  plot  was  obtained  by  integrating  the  derivative  line  shape 
measured  using  a  cw  ESR  spectrometer. 
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We  used  a  frequency  modulated  MRFM  measurement 
protocol  similar  to  that  described  by  Wago  et  al.  [20]. 
Prior  to  measurement,  the  microwave  source  was  turned 
off  and  the  spins  were  allowed  to  polarize  in  the  static 
field  i?o(R)  f°r  a  time  long  compared  to  Tx.  To  generate 
the  MRFM  signal,  the  microwave  field  was  frequency 
modulated  sinusoidally  at  the  cantilever  frequency  with 
a  peak  frequency  deviation  of  0/27 t  =  20  MHz 
(AZ?mod  =  Cl/y  =  0.7  mT),  causing  the  spins  within  the 
slice  of  thickness  £  —  2 El/yGzz  to  undergo  cyclic  adia¬ 
batic  inversions,  where  Gzz  =  dBz/ dz-  In  the  presence  of 
the  axial  field  gradient  of  the  magnetic  particle  Gzz ,  the 
oscillating  magnetization  creates  a  time  dependent  force 
that  resonantly  drives  the  cantilever  from  its  equilibrium 
position.  Hence,  the  peak  vibration  amplitude  of  the 
cantilever  (the  MRFM  signal)  is  proportional  to  the 
number  of  polarized  spins  within  the  slice.  For  Z?ext  = 
50  mT,  where  Gzz  =  15  mT/  jmm  and  the  slice  is  approxi¬ 
mately  0.5  jmm  below  the  surface,  the  resonant  slice  con¬ 
tained  roughly  N  ~  6  X  103  net  polarized  spins  resulting 
in  a  peak  detected  force  Fpeak  =  NjmBGzz  ^  8  X  10“ 16  N. 

To  calibrate  the  gradient  from  the  tip,  the  location  of 
the  resonant  slice  was  measured  as  a  function  of  external 
field  [14].  For  a  given  tip-to-surface  distance  <7surface  there 
is  a  corresponding  external  field  Z?ext  =  Z?onset  below 
which  the  slice  does  not  penetrate  the  surface  and  no 
signal  is  observed.  The  gradient  is  obtained  by  measuring 
the  shift  in  the  onset  field  Z?onset  at  two  different  locations 
of  the  surface.  We  observed  that  when  the  surface  was 
moved  closer  to  the  tip  by  Az  =  340  nm  (from  2.50  to 
2. 16  jmm)  the  onset  field  changed  by  A Bz  =  4.9  mT  (from 
45.8  to  40.9  mT).  Thus,  the  average  value  of  the  axial 
gradient  was  Gzz  =  ABZ/Az  =  15  mT / jmm  for  Z?ext  = 
43.4  mT.  From  the  known  geometry  and  the  total  moment 
of  the  magnetic  particle,  a  model  for  Z?tip(R)  was  devel¬ 
oped  which  was  consistent  with  the  measured  gradient. 
This  model  was  then  used  to  calculate  the  gradient  from 
the  tip  for  an  arbitrary  value  of  5ext. 

The  slice  penetration  depth  d  =  dslice  —  dsurface  was 
chosen  so  as  to  obtain  a  relatively  uniform  axial  gradi¬ 
ent  within  the  slice  volume,  while  having  a  detectable 
number  of  spins.  If  we  assume  a  dipole  distribution  for 
the  tip  field,  then  the  axial  gradient  is  given  by  Gzz  = 
3jm0jmtip(3cosa  —  5cos3 a) / AttRa ,  where  a  is  the  angle 
made  between  the  moment  |mtip  and  the  z  axis  and  R  is 
the  distance  from  the  center  of  the  particle  to  the  center 
of  the  slice  (see  Fig.  1).  It  is  evident  from  this  expression 
that  the  gradient  within  the  slice  is  most  uniform  for 
small  opening  angles  amSiX  —  yfd/~R,  hence  small  penetra¬ 
tion  depths.  By  making  d  small,  however,  both  the  slice 
volume  as  well  as  the  signal  decrease.  For  these  measure¬ 
ments,  we  chose  d^  0.5  jmm  which  gave  a  sufficiently 
high  signal-to-noise  ratio  while  maintaining  reasonable 
gradient  uniformity.  For  the  largest  gradients  that  were 
used,  we  estimate  that  the  average  gradient  within  the 
slice  is  within  15%  of  the  peak  value. 
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To  obtain  T { ,  the  spins  were  first  allowed  to  polarize, 
then  adiabatically  inverted  by  briefly  turning  on  the 
microwave  field  20  MHz  above  resonance  and  sweeping 
the  frequency  to  20  MHz  below  resonance.  The  spins 
were  then  allowed  to  recover  for  a  variable  period  of 
time  and  the  polarization  state  was  measured  using  the 
frequency  modulated  MRFM  protocol.  Tx  was  deter¬ 
mined  by  numerically  fitting  the  data  to  the  exponen¬ 
tial  recovery  curves  shown  in  Fig.  3.  As  can  be  seen  in 
Fig.  4,  Txl  has  a  roughly  constant  value  of  Txl  = 
0.16  s-1  for  gradients  ranging  between  0.01  <  Gzz  < 
1  mT/^tm.  For  Gzz>  1  mT/yum,  T{  1  decreases  mono- 
tonically  with  increasing  gradient,  reaching  Txl  = 
0.04  s-1  for  Gzz  =  36  mT/ jmm. 

In  order  to  explain  the  observed  dependence  of  Tx  on 
field  gradient,  we  propose  the  following  model.  In  a 
disordered  sample,  there  can  be  a  distribution  of  relaxa¬ 
tion  times  resulting  from  local  inhomogeneities  in  spin 
concentration  (e.g.,  clusters  [21])  or  from  variations  in 
coupling  strength  between  a  paramagnetic  electron  and 
the  host  defect  site  [22].  If  the  spin  temperature  of  the 
system  is  perturbed,  as  is  the  case  when  we  adiabatically 
invert  the  spin  population,  faster  relaxing  spins  rapidly 
relax  to  the  lattice  while  the  slower  relaxing  spins  remain 
inverted.  Thermal  equilibrium  within  the  in  homogene  - 
ously  broadened  spectrum  is  established  via  spin  diffu¬ 
sion  which  carries  Zeeman  energy  from  slower  to  faster 
relaxing  spins,  thus  reducing  the  overall  T ,  of  the  entire 
ensemble  [23-25].  If  the  field  gradient  from  the  magnetic 
particle  is  of  sufficient  strength  to  quench  flip-flop  inter¬ 
actions,  spin  diffusion  will  be  suppressed  and  the  overall 
relaxation  rate  will  become  slower. 

If  spin  diffusion  is  the  operative  mechanism  by  which 
thermal  equilibrium  is  established  within  the  spin  en¬ 
semble,  then  as  the  gradient  is  increased  Txx  should  begin 
to  decrease  when  the  difference  in  field  between  adjacent 
spins  exceeds  the  homogeneous  linewidth,  i.e.,  zGzz  > 


A Bh,  where  z  =  r/ 2  is  the  average  spacing  between  spins 
in  the  z  direction.  In  our  case,  r  =  8  nm  and 

A Bh  —  S.&juiQyR/ATrT3  —  1.4  X  10-2  mT  is  an  estimate 
for  the  homogeneous  linewidth  for  a  dilute  concentration 
of  spins  [26].  For  our  sample,  the  minimum  gradient 
required  to  spectrally  separate  two  overlapping  spin  pack¬ 
ets  is  approximately  Gzz  —  2.8  mT/ jmm,  which  is  consis¬ 
tent  with  the  observed  onset  gradient  in  Fig.  4. 

For  the  case  of  an  in  homogeneously  broadened  spec¬ 
trum,  the  probability  for  two  spins  with  different  g  values 
to  exchange  magnetization  while  conserving  Zeeman 
energy  is  proportional  to  the  normalized  overlap  func¬ 
tion  d>(Gzz)  [26], 


®(GZZ) 


f  -  zG-z)dB' 

f  f2(B')dB ' 


(1) 


where  f(B )  represents  the  shape  of  the  inhomogeneous 
spectrum.  The  solid  line  in  Fig.  4  is  a  plot  of  Eq.  (1) 
calculated  using  the  f(B )  shown  in  Fig.  2,  which  at 
2.96  GHz  has  a  width  —0.1  mT.  Excellent  agreement 
with  the  measured  data  is  obtained  for  z  —  5  nm.  This 
value  for  z  is  in  reasonable  agreement  with  that  estimated 
from  the  spin  concentration  (z  —  4  nm). 

An  estimate  for  the  fraction  of  fast  relaxing  spins 
nj/ns  needed  to  account  for  the  observed  decrease  in  Tx 
in  the  presence  of  spin  diffusion  can  be  derived  from  the 
work  of  Vugmeister  [25].  Within  this  model,  the  slow 
relaxing  spins  are  coupled  to  the  lattice  only  through 
cross  relaxation  with  the  fast  relaxing  spins.  In  the  limit 
t1l  «  tis ,  where  til  is  the  spin-lattice  relaxation  time 
for  the  fast  relaxing  spins  and  rIS  is  the  cross-relaxation 
time  between  fast  and  slow  relaxing  spins, 


nj_  ~  28 

ns  9Txn2syAR2 


(2) 


Recovery  time  (s) 

FIG.  3.  Polarization  as  a  function  of  recovery  time  for  three 
different  field  gradients:  •,  Gzz  ^  0.9  mT/jmm,  Tx  = 
5.6  ±0.1  s;  ■,  Gzz  ^  9  mT / /xim,  Tx  =  11  ±  0.5  s;  A,  Gzz  ^ 
27  mT/yum,  Tx  =  17.8  ±  1.2  s.  The  lines  show  exponential 
recovery  curves  fit  to  the  data. 


External  field  (mT) 
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FIG.  4.  Spin-lattice  relaxation  rate  Tf 1  versus  gradient.  The 
measured  data  (•)  is  plotted  as  a  function  of  external  field  (top 
axis).  Values  for  the  tip  gradient  (bottom  axis)  for  a  given 
external  field  are  calculated  numerically  using  a  model  for  the 
tip.  The  overlap  function  O(G^)  (solid  curve)  is  calculated 
assuming  an  average  axial  separation  z  =  5  nm  between  spins. 
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FIG.  5.  Decay  of  magnetization  in  the  rotating  frame  as  a 
function  of  the  spin-lock  time  rlock.  The  solid  curves  are  double 
exponential  fits  to  the  data  that  serve  as  visual  guides, 
(a)  Gzz  —  2  mT / /xm,  T]p  =*  0.33  s.  (b)  Gzz  —  6  mT/jum,  Tlp 
0.53  s.  (c)  Gzz  ^  12  mT/yum,  Tlp  ^  0.80  s. 

If  we  take  8/y  ^  0. 1  mT  for  the  width  of  the  inhomo- 
geneously  broadened  spectrum  and  Tx  ^  6  s  for  the  re¬ 
laxation  time  in  the  presence  of  spin  diffusion,  we  find 
ni/ns  ~  10-6*  Thus,  our  results  are  consistent  with  a 
spin  population  where  only  a  very  small  fraction  of  the 
total  spins  have  fast  intrinsic  relaxation  rates.  When  the 
spin  diffusion  process  is  suppressed,  the  slower  in¬ 
trinsic  relaxation  rate  of  the  vast  majority  of  spins  is 
revealed 

Flip-flop  interactions  in  the  rotating  frame  can  also 
give  rise  to  spin  diffusion  and  enhance  the  rotating 
frame  spin-lattice  relaxation  rate  Tfp  [27].  In  the  pres¬ 
ence  of  the  field  gradient  from  the  tip,  one  might  expect 
such  interactions  to  be  suppressed  when  the  variation 
in  the  magnitude  of  the  effective  field  |Beff|  between 
adjacent  spins  becomes  comparable  to  the  dipolar  line- 
width,  where  Beff  =  Bxx  +  zGzzz  and  z  is  the  distance 
from  the  center  of  the  resonant  slice.  For  an  average 
nearest  neighbor  spacing  z  =  5  nm  and  a  dipolar  line- 
width  of  A Bh  —  0.014  mT,  we  would  expect  TXp  to  be 
affected  for  Gzz  ^  13  mT/yam. 

To  measure  TXp,  the  polarized  spins  within  the  resonant 
slice  were  first  spin  locked  to  the  effective  field  Beff  by 
turning  on  the  microwave  field  20  MHz  above  resonance 
and  then  adiabatically  sweeping  the  microwave  frequency 
onto  resonance.  After  waiting  a  variable  time  rlock  with 
the  microwave  field  turned  on,  the  resulting  magnetiza¬ 
tion  was  read  out. 

Figure  5  shows  the  magnetization  as  a  function  of 
spin-lock  time  rlock.  The  nonexponential  decay  curves 
can  be  characterized  with  two  time  constants.  For  rlock  < 
0.2  s,  the  decay  rate  was  found  to  be  —2.5 X  faster  than 
for  subsequent  times.  As  the  gradient  from  the  tip  was 
raised  from  2  to  12  mT/jum ,  the  time  for  the  magnetiza- 
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tion  to  decay  to  e  1  of  its  initial  value  increased  from 
0.33  to  0.80  s. 

In  summary,  we  have  used  MRFM  to  investigate  and 
control  the  spin-lattice  relaxation  rate  of  E'  centers  in 
Si02.  We  find  that  Tf1  and  TXp  are  significantly  reduced 
in  the  presence  of  large  gradients  that  suppress  the  dipolar 
coupling  between  adjacent  spins.  The  ability  to  quench 
spin-spin  interactions  by  using  large  gradients  is  not  only 
relevant  for  MRFM  detection,  but  may  also  be  of  interest 
for  the  design  of  solid-state  quantum  devices. 
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D.  Pearson  for  experimental  support.  We  gratefully  ac¬ 
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Abstract — Detection  of  a  finite  state  Markov  signal  in  additive 
white  Gaussian  noise  (AWGN)  can  be  done  in  an  intuitive  manner 
by  applying  an  appropriate  filter  and  using  an  energy  detector. 
One  might  not  expect  this  heuristic  approach  to  signal  detection 
to  be  optimal.  However,  in  this  paper,  we  show  that  for  a  certain 
type  of  finite  state  Markov  signal,  namely  the  discrete-time  (d-t) 
random  telegraph,  this  filtered  energy  detector  is  approximately 
optimal  under  the  following  conditions  of:  symmetric  transition 
probabilities,  low  signal-to-noise  ratio  (SNR),  long  observation 
time,  and  small  probability  of  transition  between  two  consecutive 
time  instances.  When  these  last  three  conditions  hold,  but  the 
transition  probabilities  are  not  symmetric,  we  show  that  a  variant 
of  the  filtered  energy  detector  is  approximately  optimal.  It  is  also 
shown,  under  low  SNR  conditions,  that  the  Likelihood  Ratio 
Test  for  a  d-t  finite  state  Markov  signal  in  AWGN  reduces  to 
the  matched  filter  statistic  with  the  Minimum  Mean- Squared 
Error  predictor  used  in  place  of  the  “known”  signal  values.  This 
estimator-correlator  receiver  structure  is  well  known  in  signal 
processing  literature.  This  work  is  relevant  to  magnetic  resonance 
force  microscopy,  an  emerging  technology  that  uses  ultrasensitive 
force  sensing  to  detect  magnetic  resonance.  Sensitivity  down  to 
the  single  spin  level  was  demonstrated  in  a  recent  experiment. 

I.  Introduction 

Detection  of  a  discrete-time  (d-t)  finite  state  Markov  signal 
in  additive  white  Gaussian  noise  (AWGN)  is  widespread  in 
many  different  fields.  Detection  of  a  random  telegraph  signal 
is  used  in  the  study  of  particle  tunnelling  [1],  and  in  the 
study  of  low-frequency  noise  characteristics  of  light-emitting 
diodes  [2].  In  [3],  statistical  network  anomaly  detection  based 
on  Markov  chains  is  presented,  and  in  [4],  land  mine  detection. 
The  focus  application  of  this  paper  is  magnetic  resonance 
force  microscopy  (MRFM),  which  is  a  promising  technique 
for  three-dimensional  imaging  on  the  nanometer  scale.  Recent 
experiments  at  IBM  have  shown  that  MRFM  is  capable  of  de¬ 
tecting  and  localizing  individual  electron  spins  associated  with 
sub-surface  atomic  defects  in  silicon  dioxide  [5].  This  single¬ 
spin  detection  milestone  represents  a  factor  of  107  improve¬ 
ment  over  conventional  electron  spin  resonance  detection  and 
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was  achieved  using  energy  detection  methods  similar  to  those 
described  in  this  paper.  Other  recent  MRFM  experiments  have 
demonstrated  the  ability  to  detect  and  manipulate  naturally 
occurring  statistical  fluctuations  in  small  spin  ensembles  [6]. 
With  further  development,  single- spin  MRFM  may  eventually 
lead  to  atomic-resolution  magnetic  resonance  imaging  and  find 
application  in  quantum  computing  experiments  [7]. 

There  are  two  main  results  of  this  paper.  Firstly,  when  used 
to  detect  the  d-t  random  telegraph  in  AWGN,  the  filtered 
energy  (f-e)  detector  is  approximately  optimal  under  the  fol¬ 
lowing  four  conditions:  symmetric  transition  probabilities,  low 
signal-to-noise  ratio  (SNR),  long  observation  time,  and  a  small 
probability  of  transition  between  two  consecutive  instances. 
The  f-e  detector  is  no  longer  approximately  optimal  when 
the  transition  probabilities  are  asymmetric.  We  extend  the  f-e 
detector  to  a  hybrid  second-order  detector  which  combines  the 
filtered  energy,  amplitude,  and  energy  statistics.  It  is  shown 
that  the  hybrid  detector  is  approximately  optimal  for  the  d-t 
random  telegraph  model  under  only  the  last  three  conditions. 
The  second  result  of  this  paper  is  a  new  interpretation  of 
the  optimal  Likelihood  Ratio  Test  (LRT)  for  a  d-t  finite 
state  Markov  signal  under  low  SNR  conditions.  It  is  shown 
that,  under  low  SNR,  the  LRT  reduces  to  the  matched  filter 
statistic  with  the  Minimum  Mean-Squared  Error  (MMSE) 
predictor  used  in  place  of  the  “known”  values.  Current  single 
spin  experiments  operate  under  conditions  of  very  low  SNR; 
consequently,  we  are  interested  in  the  performance  of  detectors 
in  the  regime  of  low  SNR  and  long  observation  time.  Our 
results  are  therefore  applicable  to  [1],  [2],  [3],  [4]  under 
conditions  of  low  SNR. 

Equations  for  the  recursive  implementation  of  the  LRT  can 
be  found  in  [8].  The  recursive  form  of  the  LRT  is  used  to 
derive  the  two  reported  results.  The  first  result  is  intuitively 
pleasing:  the  idea  of  performing  detection  of  a  finite  state 
Markov  signal  by  filtering  the  noisy  observations  and  applying 
an  energy  detector  is  one  that  comes  naturally.  Moreover,  we 
present  simulations  that  suggest  the  optimality  of  a  similar 
approximation  for  the  d-t  random  walk  process.  A  similar  sim¬ 
plification  effect  has  been  observed  in  the  linear  vs.  nonlinear 
estimation  of  a  noisy  continuous-time  (c-t)  random  telegraph 
process.  See  Remark  1  in  [9].  The  second  result  has  a  c-t 
analog:  in  c-t,  the  LRT  for  detecting  a  random  signal  in  AWGN 
has  the  form  of  the  matched  filter  statistic  with  the  MMSE 
predictor  used  in  place  of  the  “known”  signal.  It  is  exact 
under  all  SNR  conditions  [10],  [11].  In  d-t,  this  estimator- 
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correlator  receiver  structure  appears  in  the  LRT  of  problems 
whose  probability  density  functions  (pdfs)  have  an  exponential 
form  [12].  In  particular,  it  applies  when  wanting  to  detect  a 
Gaussian  signal  in  AWGN. 

The  outline  of  this  paper  is  as  follows.  In  Section  II,  we 
briefly  review  the  basic  principles  of  MRFM.  This  is  followed 
by  a  discussion  in  Section  III  of  two  d-t  signal  models:  the 
random  telegraph  and  random  walk  models.  In  Section  IV,  we 
describe  existing  detectors  that  are  commonly  used,  namely  the 
amplitude  and  filtered  energy  detectors,  and  compare  them  to 
the  optimal  detectors.  We  derive  a  new  interpretation  of  the 
LRT  under  low  SNR  conditions.  As  well,  the  f-e  detector  is 
extended  to  a  hybrid  version.  Simulation  results  are  presented 
in  Section  V. 

II.  Basic  Principles  of  MRFM  spin  detection 

MRFM  experiments,  in  general,  involve  the  measurement 
of  magnetic  force  between  a  submicron- size  magnetic  tip 
and  spins  in  a  sample.  The  details  of  spin  manipulation 
and  signal  detection  depend  on  the  exact  MRFM  protocol 
used.  One  particularly  successful  protocol  is  called  OSCAR, 
which  stands  for  OScillating  Cantilever-driven  Adiabatic  Re¬ 
versal  [13],  [14].  A  variation  of  this  protocol,  “interrupted 
OSCAR”,  was  used  in  recent  single  spin  experiments  [5]. 

A  schematic  diagram  of  an  OSCAR-type  MRFM  experi¬ 
ment  is  shown  in  Fig.  1.  As  shown  in  the  figure,  a  submicron 
ferromagnet  is  placed  on  the  tip  of  a  cantilever  and  positioned 
close  to  an  unpaired  electron  spin  contained  within  the  sample. 
An  applied  radio-frequency  (rf)  field  serves  to  induce  magnetic 
resonance  of  the  spin  when  the  condition  Bq  =  cjrf /7  is  met. 
Here,  Bq  is  the  magnitude  of  the  magnetic  field  from  the  tip, 
plus  any  externally  applied  static  field  that  may  be  present. 
Because  the  magnetic  field  emanating  from  the  tip  is  highly 
inhomogeneous,  magnetic  resonance  is  confined  spatially  to  a 
thin  bowl-shaped  region  called  the  “resonant  slice”. 


Fig.  1.  Schematic  of  the  iOSCAR  experiment. 

In  an  OSCAR  experiment,  a  gain-controlled  positive  feed¬ 
back  loop  is  used  to  oscillate  the  cantilever  with  a  preset 
amplitude  (typically  10-20  nm).  The  cantilever  oscillation 
frequency  is  determined  by  the  cantilever  itself  (specifically, 
by  the  fundamental  flexural  mode  eigenfrequency),  as  well  as 
by  tip- sample  interactions.  As  the  tip  of  the  cantilever  vibrates, 
the  resonant  slice  passes  back  and  forth  through  the  spin  and, 


as  a  result,  the  spin  direction  is  cyclically  inverted  due  to 
an  effect  called  adiabatic  rapid  passage  [13],  [15],  [16].  The 
cyclic  inversion  is  synchronous  with  the  cantilever  motion 
and  affects  the  cantilever  dynamics  by  slightly  shifting  the 
cantilever  resonance  frequency.  The  frequency  shift  depends 
on  the  angle  0  of  the  spin  with  respect  to  a  vector  called  the 
“effective  field  in  the  rotating  frame”.  See  [14],  [15]  for  further 
details.  The  frequency  shift  can  be  written  as 

Acj  Acjmax  cos  0,  (1) 

where  Ao;max  =  2ujqG p / 7rkxp^  [14].  Here  cjo  is  the  unper¬ 
turbed  cantilever  frequency,  G  is  the  field  gradient  from  the 
tip,  /i  is  the  magnetic  moment  of  the  spin,  k  is  the  cantilever 
spring  constant  and  xv k  is  the  peak  amplitude  of  the  cantilever 
vibration.  The  factor  cos  6  represents  the  normalized  projection 
of  the  spin  in  the  direction  of  the  effective  field. 

There  are  several  impediments  to  single  spin  detection. 
Firstly,  because  the  force  from  a  single  spin  is  so  tiny  (a 
few  attonewtons),  the  maximum  cantilever  frequency  shift 
is  only  about  one  part  per  million  for  typical  experimental 
parameters.  This  small  frequency  shift  must  be  detected  in 
the  presence  of  the  cantilever  phase  (or  frequency)  noise 
that  originates  from  cantilever  thermal  vibrations  and  sample- 
induced  force  fluctuations.  The  resulting  low  SNR  necessitates 
long  integration  times  for  signal  detection.  Secondly,  the 
detection  is  complicated  by  environmental  disturbance  to  the 
spin  (i.e.,  relaxation  effects)  that  can  randomly  flip  the  spin 
orientation  and  reverse  the  signal  polarity  during  the  signal 
integration  time.  A  low  operating  temperature,  on  the  order  of 
IK,  can  help  reduce  both  the  cantilever  thermal  excitations  and 
the  random  spin  flip  rate.  Nevertheless,  the  signal  processing 
methodology  must  take  these  effects  into  account. 

We  consider  two  models  for  the  spin  behaviour  in  the  pres¬ 
ence  of  environmental  disturbances.  In  a  quantum  mechanical 
measurement  model,  the  spin  is  always  found  to  be  either 
aligned  or  anti-aligned  with  the  effective  field,  so  that  cos  0  = 
±1.  Thus  the  frequency  shift  signal  has  only  two  levels: 
Auj  =  ±A cjmax,  and  the  time  sequence  of  the  frequency 
shift  is  a  random  telegraph  signal  with  a  transition  rate  that 
depends  on  the  spin  relaxation  rate.  In  a  classical  measurement 
model  (which  we  consider  for  the  sake  of  completeness),  cos  0 
can  take  arbitrary  values  between  +1  and  -1.  As  a  result  of 
environmental  disturbance,  Ac u  will  be  uniformly  distributed 
between  —  Acjmax  and  +Acjmax.  For  this  case,  a  bounded 
random  walk  model  is  appropriate. 

III.  MRFM  SIGNAL  MODELS 
A.  Model  1:  Discrete-time  Random  Telegraph  model 

In  the  quantum  measurement  model,  the  frequency  shift  is 
characterized  by  random  jumps  between  two  discrete  levels. 
The  jumps  are  taken  to  be  Poisson  distributed.  Denote  the 
d-t  random  telegraph  signal  by  Q,  where  ti  =  iTs  are  the 
sampling  times,  and  Ts  is  the  sampling  time  interval.  In  this 
paper,  a  Markovian  process  with  a  finite  number  of  states 
will  have  a  state  space  denoted  by  4/  =  {^i, . . . ,  Vv}>  where 
r  is  the  number  of  states.  Let  the  state  space  of  the  d-t 
random  telegraph  be  it  has  r  =  2  states  and  we  shall 
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take  ipi  =  —A,ip2  =  A,  where  A  is  the  amplitude  of  the 
random  telegraph  (A  corresponds  to  Atumax  for  the  case  of  a 
MRFM  signal).  As  an  initial  condition,  (0  is  equally  likely  to 
be  either  ±A.  Then,  a  probability  transition  matrix  Prt  can  be 
associated  with  Q  such  that  the  (j,  k)- th  value  of  Prt  equals 
P{Ci  —  'ipk I Ci— i  =  ^j)  f°r  1  <  j,  &  <  2  and  i  >  1.  Assume 
that  Prt  has  the  form: 


Prt  = 


q  1  -q 

1  —  p  p 


(2) 


where  0  <  p,  q  <  1.  If  p  =  q,  we  say  that  the  transition 
probabilities  are  symmetric ,  whereas  if  P  ^  q,  we  shall 
say  that  they  are  asymmetric.  Define  the  signal  vector  (  = 
[Co,  •  •  • ,  CiV-i]T>  the  noise  vector  w  =  [wo, . . . ,  rc/v- i]t, 
and  the  observation  vector  y  =  [yo,  •  •  • ,  Vn-i]T,  where  the 
superscript  (*)T  denotes  the  transpose  operator.  The  Wi  s  are 
modelled  as  independent  and  identically  distributed  (i.i.d.) 
Gaussian  random  variables  (r.v.s)  with  zero  mean  and  variance 
a2.  The  detection  problem  is  then  to  decide  between: 


Ho  (spin  absent) 
Hi  (spin  present) 


y  =  w 
y  =  (  +  w 


(3) 


Examples  of  noiseless  and  noisy  random  telegraph  signals 
are  given  in  Fig.  2.  For  the  random  telegraph  signal,  the  SNR 
is  defined  to  be  SNR  =  A2 /cr2.  The  SNR  in  dB  is  defined  in 
the  usual  way  as  SNRdB  =  10  log10  SNR. 


(c)  nP 


I ;  we  shall  refer  to  this  model  as  the  d-t  random  walk  model. 
The  state  space  4/rw  of  the  d-t  random  walk  will  then  have 
r  =  2 M  +  1  states,  where  i/jj  =  (j  —  M  —  l)s  for  j 

1 . (2 M  +  1).  Associate  with  Ci  the  probability  transition 

matrix  Prw,  so  that,  as  before,  the  (j,  k)- th  element  of  Prw  is 
P(Ci  =  V’/cICi-i  =  i>j)  for  1  <  j,  k  <  (2 M  +  1)  and  i  >  1. 
Prw  is  defined  such  that,  at  each  time  step,  Q  changes  by 
either  ±s.  We  assume  reflecting  boundary  conditions,  and  Co 
is  equally  likely  to  be  either  ±s.  These  conditions  imply  that 
Prw  is  a  tridiagonal  matrix. 

The  detection  problem  is  now  to  test  (3)  when  (  is  modelled 
by  a  random  walk.  Note  that  the  d-t  random  walk  model  can 
be  regarded  as  a  multi-state  generalization  of  the  d-t  random 
telegraph  model.  In  the  limit  as  s  —>  0,  M  —>  oo,  the  random 
walk  converges  to  Brownian  Motion  over  the  interval  I  [17]. 

Examples  of  noiseless  and  noisy  random  walk  signals  are 
given  in  Figs.  3  and  4,  where,  at  each  state,  a  change  of 
±8  is  equally  likely.  Here,  we  shall  define  the  SNR  as 
SNR  =  (lim^oo  E[(2])/a2,  where  Q  is  the  d-t  random 
walk  process.  In  other  words,  the  SNR  is  the  ratio  of  the 
steady-state  expected  energy  of  Q  to  the  noise  variance.  Note 
that  this  definition  is  consistent  with  that  provided  for  the 
random  telegraph  signal.  If  Q  represented  the  d-t  random 
telegraph  signal,  then  Q  =  A2  at  all  time  instances,  leading 
to  SNR  =  A2 /a2. 


0  50  100  150  200  250  300  350  400  450  500 

time  index 


Fig.  2.  a:  Noiseless  random  telegraph  signal  with  symmetric  transition 
probabilities  p  =  q  =  0.98;  b:  Noisy  version  of  (a)  at  SNR  =  -3  dB;  c: 
Noiseless  random  telegraph  signal  with  asymmetric  transition  probabilities 
p  =  0.98,  q  =  0.6;  d:  Noisy  version  of  (c)  at  SNR  =  -3  dB 


Fig.  3.  a:  Noiseless  random  walk  signal  with  5  levels;  b:  Noisy  version  of 
(a)  at  SNR  =  -7.3  dB 


IV.  Detection  strategies 


B.  Model  2:  Discrete-time  Random  Walk  model 

In  the  classical  spin  detection  model,  the  frequency  shift 
signal  is  well  approximated  by  a  one  dimensional  random  walk 
confined  to  the  interval  /  =  [—A,  A],  where  A  =  Acjmax  for 
the  case  of  a  MRFM  signal.  We  discretize  I  into  (2 M  +  1) 
states  using  a  step  size  of  s ,  where  MgZ  and  M,  s  >  0  and 
define  Q  to  be  the  random  walk  restricted  to  the  discretized 


The  detectors  considered  here  can  be  placed  into  three 
categories:  versions  of  existing  detectors  that  are  currently  in 
use  for  MRFM;  LRTs;  and  approximations  to  the  LRT  for 
the  d-t  random  telegraph  model.  The  LRT  is  a  most  powerful 
test  that  satisfies  the  Neyman-Pearson  criterion:  it  maximizes 
the  probability  of  detection  (Pd)  subject  to  a  constraint  on 
the  probability  of  false  alarm  (Pp)  [18],  which  is  set  by  the 
user.  Consequently,  it  can  be  used  as  a  benchmark  with  which 
to  compare  the  other  detectors.  When  the  random  transition 
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Fig.  4.  a:  Noiseless  random  walk  signal  with  21  levels;  b:  Noisy  version  of 
(a)  at  SNR  =  -7.8  dB 


d-t  random  telegraph  model.  When  the  transition  probabilities 
are  symmetric,  this  figure  corresponds  to  the  mean  number  of 
transitions  per  second  A  of  the  c-t  random  telegraph  model. 
Given  a  value  p  =  q  =  po,  one  can  equate  the  expected  number 
of  transitions  in  both  d-t  and  c-t  models  to  obtain  A  =  (1  — 
Po)/Ts.  In  practice,  since  the  mean  number  of  transitions  is 
only  approximately  known  to  the  experimenter,  a  bank  of  LPFs 
with  different  a’s  are  used  to  perform  detection  [5].  The  energy 
and  filtered  energy  detector  can  then  be  expressed  as 

N- 1  Hi 

77  (Vi  *  hi)2  ^  y  (6) 

i= 0  Hq 

where  represents  the  convolution  operator.  For  the  energy 
detector,  hi  is  taken  to  be  the  unit  impulse  function  S[i\, 
while  for  the  filtered  energy  detector,  hi  =  hp?  [i] ,  the  impulse 
response  of  HLF(z)  in  (5). 

Note  that  the  computational  complexity  for  the  amplitude, 
filtered  energy,  and  energy  detectors  is  O(N). 


times  are  known,  the  optimal  LRT  is  the  matched  filter,  called 
the  omniscient  matched  filter  (MF)  in  this  paper.  Although 
unimplementable,  the  MF  detector  provides  an  absolute  upper 
bound  when  comparing  the  various  detectors’  Receiver  Oper¬ 
ating  Characteristic  (ROC)  curves. 


A.  Amplitude,  energy,  filtered  energy  detectors 
The  d-t  amplitude  detector  is 


1 

N 


N-l 


52  y< 


Hx 

it  V 

H0 


(4) 


where  p  is  set  to  satisfy  the  constraint  on  Pp.  This  is  the 
optimal  test  under  the  assumption  that  yi  is  the  sum  of  an 
unknown  constant  and  AWGN.  This  assumption  would  be  true 
if  there  were  no  random  spin  flips.  In  this  case,  the  amplitude 
detector  is  simply  a  MF  detector.  However,  as  the  number 
of  random  transitions  in  yL  increases,  the  performance  of  the 
amplitude  detector  degrades.  An  alternative  test  statistic  is 
the  d-t  signal  energy,  i.e.  the  sum  of  the  squares  of  the  yi 
instead  of  the  magnitude  of  the  sum  in  (4).  As  the  signal  and 
noise  are  assumed  to  be  independent,  under  hypothesis  Hi , 
one  would  expect  y  to  have  a  higher  energy  on  average  than 
under  hypothesis  H0.  This  can  be  reliably  detected  under  a 
sufficiently  high  SNR.  A  natural  improvement  to  the  energy 
detector  is  to  reject  out-of-band  noise  by  prefiltering  y  over  the 
signal  passband.  As  the  signal  Q  is  baseband,  a  lowpass  filter 
(LPF)  is  appropriate.  In  particular,  one  might  use  a  simple 
first-order,  single-pole  filter  given  by 


H^{z) 


1  —  a  1  +  z  1 

2  1  —  OLZ~X 


(5) 


where  we  require  |a|  <  1  for  stability  [19].  The  time  constant 
a  should  be  chosen  based  on  the  bandwidth  of  the  signal; 
if  ujc  is  the  desired  -3dB  bandwidth  of  the  filter,  one  should 
set  a  =  (1  —  sino;c)/ coso;c.  The  -3dB  bandwidth  depends 
on  the  mean  number  of  transitions,  i.e.  (1  —  p)/Ts  for  the 


B.  Optimal  LRT  detectors  and  their  approximations 

A  recursive  form  of  the  LRT  for  a  d-t  finite  state  Markov 
signal  in  AWGN  is  derived  in  [8].  As  well,  it  is  possible  to 
obtain  a  closed- form  expression  for  the  LRT:  the  LRT  can  be 
expressed  as  a  product  of  matrices  pre-multiplied  by  a  row 
vector  of  the  initial  state  probabilities  and  post-multiplied  by 
a  column  vector  of  all  ones  [20].  Here,  we  focus  on  recursive 
implementations  of  the  LRT  for  the  d-t  random  telegraph  and 
random  walk  process.  See  [8],  [20]  for  more  details. 

Define  yl  =  [yo, . . .  ,^]T  for  i  >  0.  Let  Eq[-]  denote 
the  expectation  under  the  H0  (spin  absent)  hypothesis,  and 
Ei[-]  the  expectation  under  the  Hi  (spin  present)  hypothesis. 
Consider  first  the  d-t  random  telegraph  model:  define  Rfiu)  = 
P(Q  =  uly1-1),  where  u  G  47 rt  and  i  >  1.  There  exists  a 
recursive  formula  for  RfiA)  and  Rfi—A): 


(  Ri(-A) 
V  Ri(A) 


= p7  ( 1  ^  ^ )  where 

= _ efiyyp-ijA) _ 

efr™-1  Ri-i(A)  +  e“  R^-A) 

(7) 


for  i  >  1  and  with  initial  conditions  Ro(A)  =  Ro(—A)  = 
1/2.  Let  fi(y)  be  the  pdf  of  y  under  hypothesis  Hi  and 
fo{y)  be  the  pdf  of  y  under  hypothesis  Ho.  The  LRT  is 

Art (y)  =  fi(y)/fo(yj-  Define  ip(x;  fi,  a2)  =  exp[-(x  - 
H)2/2a2]/VHr<J.  Applying  (7), 

A^-1)  =  fi(yN-i\yN-2)  ■  ■  ■  fi(yi\yo)fi(yo) 

N- 1 

=  H  [RiiAjtpiVuXv2)  +  Ri(-A)ip(yi;-A,(T2)] 

i= 0 

Using  this,  the  LRT  for  the  d-t  random  telegraph  can  be 
expressed  as: 

N- 1  Hi 

K(y) II  k'(A)e^?"  •  />’,(  ,1)C  ;'M  Sty  (8) 

*=o  Ho 
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Note  that  the  running  time  of  (8)  is  0(N ),  where  N  is  the 
number  of  observations.  Under  low  SNR  conditions  (|A/cr|  <C 
1),  the  log  LRT  becomes 


N-i 

log  Art  (y)  ~  Y  log 
2=0 
N—l 

=  Ylo& 

2= 0 


Ri(A)(  1  H 2  Hi)  +  —  A){1 - 2^) 

crz  <jz 

1  +  \yi(ARi(A)  -  ARi(-A)) 

az 


where  we  use  the  approximation  e 6  ~  1  +  5  for  small  S. 
Now,  AR0(A)  -  ARo(-A)  =  [Co]  =  Co  and  ARi(A)  - 
ARi(-A)  =  E^Qy^1]  =  &  for  i  >  1,  i.e.  it  is  the  MMSE 
predictor  of  Q  given  the  past  observations  yl~x . 


N-l 

log  An (y)  ~  ^  log 

i= 0 


N-l 


-^Yy&  (9) 


2= o 


where  the  approximation  log(l  +  5)  ~  S  for  small  S  was 
used.  Under  low  SNR,  the  LRT  is  effectively  correlating  the 
observations  yi  with  the  MMSE  predictor  of  Q.  This  is  the 
matched  filter  statistic  with  the  MMSE  predictor  substituted 
in  place  of  the  “known”  values  of  Q,  and  is  known  as 
the  estimator-correlator  detector.  In  particular,  the  estimator- 
correlator  structure  is  known  to  be  optimal  for  Gaussian 
signals  in  AWGN  [12]. 

Under  the  regime  of  low  SNR,  long  observation  times  (N 
1),  and  p  -\-  q  —  l  —  l  (the  probability  of  transition  between 
consecutive  samples  is  small),  the  second-order  expansion  of 
log  Art  (y)  is  approximately  equal  to  the  hybrid  detector  with 
the  test  statistic 

Y^Vi * M*])2  +  1  2qQ  Ca  YVi  +  1  2a  Ce  ^  (10) 

i  i  i 

where  the  constants  Ca ,  Ce  are  given  in  the  appendix,  and 
a  =  p  +  q  —  1.  Therefore,  in  the  aforementioned  regime,  one 
expects  the  hybrid  detector  to  have  performance  similar  to  the 
optimal  LRT.  When  p  =  q,  the  second-order  expansion  of  the 
LRT  is  approximately  equal  to  the  filtered  energy  detector.  See 
the  appendix  for  more  details.  In  light  of  the  computational 
complexities  for  the  filtered  energy,  energy,  and  amplitude 
detectors,  the  complexity  of  (10)  is  also  O(N). 

Next,  consider  the  d-t  random  walk.  Define  the  vectors: 


(  Ri(-Ms)  \ 

/  fw(Vi  +  Ms)  \ 

R,  = 

Ri(  0) 

,m  = 

fwiVi) 

\  Ri(Ms) 

V  fw(Vi  ~  Ms)  j 

for  i  >  0 

.  Here,  Ri(u) 

=  P(Ci 

=  G 

and  fw(-) 

=  <£>(•;  0,  <72). 

If  a  = 

(ai,...,a„)T  and 

(&i, . . . ,  bn)T,  define  the  operation  a© 6  =  (ai&i, . . . ,  anbn)T . 
There  exists  a  recursive  formula  for  R{  that  is  similar  to  that 
for  the  d-t  random  telegraph  process,  cf.  (7): 


R:  =  PrJ 


Wi^QRi- 


(id 


for  i  >  1  and  with  the  initial  condition  R0  =  |(eM_1  + 
— m+i)’  where  {e^}^  are  the  standard  basis  vectors  for 
M2M+i.  The  LRT  for  the  d-t  random  walk  can  then  be 
expressed  as 


N-l 


KAv)  -  n 


rTws 


r=o 


H i 

i: 

Hn 


(12) 


The  running  time  of  the  LRT  for  the  d-t  random  walk 
is  0(MN ),  as  Prw  is  a  tridiagonal  matrix.  As  with  the 
d-t  random  telegraph  model,  an  insightful  understanding  of 
Arw(y)  can  be  gained  by  considering  the  regime  of  low  SNR. 
Lollowing  the  same  steps  that  were  used  to  obtain  (9): 

N-1  I”  r  1  1  - 

log Arw(y)  ~  Y  loe  YRi<^^1  ~ 

i= 0  ^  j=l 

-j^Y  ViEiiCily^1]  -  ^2  Y  ^i[C2|yi_1] 

i=0  i=0 

(13) 

We  see  that  (13)  is  the  matched  filter  statistic,  but  with  the 
MMSE  predictor  of  Q  and  Q  used  instead  of  the  “known” 
values.  Since  (11),  (12),  and  (13)  are  independent  of  the 
structure  of  Prw,  this  result  applies  to  all  finite  state  Markov 
processes  under  low  SNR  conditions.  It  is  interesting  to 
note  that  the  c-t  analog  is  given  in  [10],  [11].  In  the  c-t 
case,  however,  the  result  is  exact,  while  in  the  d-t  case,  we 
have  shown  the  result  approximately  holds  under  low  SNR 
conditions. 

This  result  is  consistent  with  the  approximation  that  was 
earlier  derived  for  the  d-t  random  telegraph  LRT  in  (9).  Since 
(—A)2  =  A2,  EiKfly1-1]  =  A2  for  the  d-t  random  telegraph. 
So  the  second  term  in  (13)  will  become  a  constant  independent 
of  the  observations  y  and  can  be  omitted. 


V.  Simulation  results 

The  objective  in  this  section  is  to  compare  the  detection 
methods  discussed  in  the  previous  section.  The  class  of  LRT 
detectors  is  optimal  for  their  respective  signal  models,  and 
provides  a  good  comparison  benchmark.  Comparison  of  the 
various  detectors  is  done  using:  (1)  ROC  curves,  each  of 
which  is  a  plot  of  probability  of  detection  (Pd)  vs.  probability 
of  false  alarm  (Pf),  and  (2)  power  curves,  each  of  which 
is  a  plot  of  Pd  vs.  SNR  at  a  fixed  Pp.  Some  of  the 
parameters  used  in  the  simulation  of  the  d-t  random  telegraph 
and  random  walk  models  are  as  follows:  k  =  10-3  N  m-1, 
w0  =  2tt  •  104  rad  s“\  By  =  0.2  mT,  G  =  2  •  106  T  m  \ 
The  sampling  period  was  Ts  =  1  ms,  and  signal  durations  of 
T  =  60  s  and  T  =  150  s  were  used.  The  performance  of 
the  detectors  varies  as  a  function  of  T;  in  general,  a  larger  T 
results  in  better  performance.  Values  of  T  used  in  iOSCAR 
MRLM  experiments  are  on  the  order  of  tens  of  hours  [5]. 
Nevertheless,  the  comparative  results  obtained  from  using  the 
two  values  of  T  above  are  representative  of  larger  values. 
Indeed,  our  approximations  to  the  optimal  detectors  improve 
with  increased  T. 


Appendix  paper  #6 


SUBMISSION  TO  IEEE  TRANSACTIONS  ON  SIGNAL  PROCESSING  (EDICS  CATEGORIES:  1-MRFM,  2-DETC,  2-NDIR)  6 


A.  Discrete-time  random  telegraph  model 

First,  consider  the  d-t  random  telegraph.  Fig.  5  depicts 
the  simulated  ROC  curves  at  SNR  =  -35  dB,  A  =  0.5  s_1, 
and  with  symmetric  transition  probabilities  (p  =  q).  With 
Ts  =  1  ms,  this  results  in  p  =  q  =  0.9995.  We  examine  the 
matched  filter,  d-t  random  telegraph  LRT  (RT-LRT),  filtered 
energy,  hybrid,  amplitude,  and  unfiltered  energy  detectors. 
The  RT-LRT,  filtered  energy,  and  hybrid  detector  curves  are 
virtually  identical,  which  is  consistent  with  our  analysis.  The 
unfiltered  energy  and  amplitude  detectors  have  performance 
that  is  poorer  than  the  RT-LRT,  as  it  should  be  since  the  RT- 
LRT  is  the  optimal  detector.  The  unfiltered  energy  detector 
has  the  worst  performance  out  of  the  five  detector  methods 
considered.  Lastly,  the  omniscient  MF  detector  has  the  best 
performance. 

A  power  curve  was  generated  over  a  range  of  SNRs  under 
the  same  conditions  as  before  with  a  fixed  PF  =  0.1;  it  is 
illustrated  in  Fig.  6.  For  spin  detection,  an  acceptable  range  for 
PF  is  on  the  order  of  0.05  to  0.1.  The  RT-LRT,  filtered  energy, 
and  hybrid  detector  have  similar  performance  from  -30  dB 
to  -45  dB.  With  this  particular  value  of  PF  and  A,  the  RT- 
LRT,  filtered  energy,  and  hybrid  detector  perform  from  5  dB 
to  10  dB  worse  than  the  MF  detector.  Although  the  amplitude 
detector  has  worse  performance  than  the  RT-LRT  and  filtered 
energy  detector,  all  three  have  comparable  performance  at 
-45  dB. 


Fig.  5.  Simulated  ROC  curves  for  the  d-t  random  telegraph  model  with 
symmetric  transition  probabilities  at  SNR  =  -35  dB,  T  =  60  s,  and  A  = 
0.5  s-1  for  the  omniscient  matched  filter,  d-t  random  telegraph  LRT  (RT- 
LRT),  filtered  energy,  hybrid,  amplitude,  and  unfiltered  energy  detectors.  The 
RT-LRT  is  theoretically  optimal. 

Fig.  7  shows  the  power  curves  generated  using  the  bigger 
value  of  T  =  150  s.  Again,  the  RT-LRT,  filtered  energy,  and 
hybrid  detectors  have  the  same  performance  from  -30  dB 
to  -45  dB.  Note  that  the  values  of  PF  have  increased  as 
compared  to  Fig.  6. 

The  ROC  and  power  curve  simulations  were  repeated  with 
different  values  of  A,  and  the  same  relative  performance  was 
observed.  In  the  interest  of  space,  however,  they  will  not  be 
shown.  Note  that  performance  degrades  as  A  increases  while 
Ts  is  held  constant. 


Fig.  6.  Simulated  power  curves  (Pd  vs.  SNR)  for  the  d-t  random  telegraph 
model  with  Pp  fixed  at  0.1  and  A  =  0.5  s-1,  T  =  60  s.  The  RT-LRT  is 
theoretically  optimal. 


Fig.  7.  Simulated  power  curves  (Pd  vs.  SNR)  for  the  d-t  random  telegraph 
model  with  Pp  fixed  at  0.1  and  A  =  0.5  s-1,  T  =  150  s.  The  RT-LRT  is 
theoretically  optimal. 


In  the  second  set  of  simulations,  we  investigate  the  case  in 
which  the  transition  probabilities  are  asymmetric,  i.e.  p  ^  q. 
Consider  the  scenario  where  p  =  0.9998,  q  =  0.9992,  and  all 
of  the  other  parameters  values  are  unchanged.  The  ROC  curves 
for  these  parameter  values  are  presented  in  Fig.  8.  There  are 
noticeable  differences  between  the  curves  of  the  RT-LRT  and 
filtered  energy  detectors.  The  hybrid  detector’s  curve  is  slightly 
below  that  of  the  LRT,  and  it  is  better  than  that  of  the  filtered 
energy  detector.  In  fact,  the  filtered  energy  detector  has  worse 
performance  than  the  amplitude  detector.  An  asymmetry  in  p,  q 
leads  to  a  non-zero  mean  signal,  which  is  probably  why  the 
amplitude  detector’s  performance  improves.  Indeed,  for  the  d-t 
random  telegraph  model,  lim^oo  E[Q]  =  A2p_~(^q  =  0.6 A 
for  the  values  of  p  and  q  used  here.  Asymmetric  transition 
probabilities  can  arise  in  some  situations,  e.g.  the  feedback- 
cooling-of-spins  MRFM  protocol  proposed  by  Budakian  [6]. 

Power  curves  from  SNR  =  -55  dB  to  -35  dB  were  generated 
for  the  asymmetric  case  in  Fig.  9.  A  larger  value  of  T  is 
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Fig.  8.  Simulated  ROC  curves  for  the  d-t  random  telegraph  model  with 
asymmetric  transition  probabilities  (p  =  0.9998,  q  =  0.9992)  at  SNR  = 
-45  dB,  T  =  150  s.  The  RT-LRT  is  theoretically  optimal. 


required  when  p  ^  q  for  the  hybrid  detector  to  be  a  good 
approximation  to  the  optimal  LRT;  hence,  T  =  150  s  was 
used  for  simulations  of  the  asymmetric  random  telegraph 
model.  The  hybrid  detector  has  better  performance  than  the 
amplitude  and  filtered  energy  detectors.  It  has  performance 
that  is  comparable  to  the  RT-LRT  for  lower  SNR  values. 


Fig.  9.  Simulated  power  curves  (Pd  vs.  SNR)  for  the  d-t  random  telegraph 
model  with  Pp  fixed  at  0.1,  p  =  0.9998,  q  =  0.9992,  and  T  =  150  s.  The 
RT-LRT  is  theoretically  optimal. 


B.  Discrete-time  random  walk  model 

Recall  that  for  the  d-t  random  walk  model,  Prw  is  tridi¬ 
agonal.  For  the  simulations,  a  particular  subset  of  tridiagonal 
matrices  was  studied.  Suppose  for  the  moment  that  M  is  even. 
Recall  that  the  random  walk  Q  is  confined  to  the  interval 
{—Ms,Ms\.  Define  the  lower-quartile  transition  probabilities 
as  pi,hPi,2  and  the  upper-quartile  transition  probabilities  as 
Pu,hPu, 2-  Let  Prw’^  be  the  (j,  k)-th  element  of  Prw.  Here, 
we  examine  the  performance  of  the  detectors  assuming  the 
following  reflecting  boundary  conditions:  Prw’2^  =  1,P^  = 


0  for  i  ±  2  and  P^M+1’2M)  =  1,  Pr(2M+1’i}  =  0  for  i  ±  2 M. 
The  rest  of  Prw  is 


pCfi*0 


Pi,  i 
Pi,  2 

0.5 

Pu,  1 
Pu,  2 


2<  j  <  M/2  +  1,  fe  =.7-1 
2  <  j  <  M/2  +  1,  k  =  j  +  1 
M/2  +  1  <  j  <  3M/2  +  1,  k  =  j  ±  1 
3M/2  +  1  <  j  <  2M,  k  =  j  -  1 
3M/2  +  1  <  j  <  2M,  k  =  j  +  1 

(14) 


Let  An(pi,p2)  be  a  n  x  (n  +  2)  matrix  that  looks  like: 


A„(pi,p2) 


t  Pi  0  p2 

Pi  o  P2 


\ 


V 


Pi  0  P2  j 


where  the  unspecified  parts  of  the  matrix  are  taken  to  be 
all  zeros.  In  this  section,  the  following  subset  of  transition 
matrices  for  the  d-t  random  walk  was  studied: 


(  0  1 

A|4(Piti,Pi)2) 

F 

Am_1  {j)u,  !•>  Pu,2  ) 

V  2  i 


\ 


0  / 


where  F  =  Am+i(0.5,  0.5).  Note  that  since  each  row  of  a 
probability  transition  matrix  must  sum  to  1,  one  has  pip  + 

Pi, 2  =  1  and  pUi i  +pu,2  =  I- 

In  the  case  of  M  odd,  the  ranges  for  the  indices  j,  k  would 
change  in  an  obvious  way.  When  pip  =  pu?2  p^2  = 

Pu,  i,  we  say  that  the  transition  probabilities  are  symmetric , 
and  if  not,  that  they  are  asymmetric.  The  matched  filter,  d-t 
random  walk  LRT  (RW-LRT),  d-t  random  telegraph  LRT 
(RT-LRT),  filtered  energy,  amplitude,  and  unfiltered  energy 
detectors  are  compared.  In  order  to  run  the  RT-LRT  in  the  case 
of  the  symmetric  d-t  random  walk,  an  average  autocorrelation 
function  of  the  random  walk  was  empirically  generated;  then  p 
was  selected  (and  one  used  q  =  p)  so  that  the  autocorrelation 
function  of  the  symmetric  d-t  random  telegraph  matched  the 
empirical  result.  From  this,  the  optimal  a  for  the  LPF  of  the 
filtered  energy  detector  was  also  obtained. 

The  ROC  curves  for  two  symmetric  cases  are  illustrated 
in  Figs.  10  and  11.  In  the  former,  p^  1  =  p^2  —  Pu, l  = 
pu, 2  =  0.5,  while  in  the  latter,  pp%  =  pu,2  =  0.52  and 
Pi, 2  =  Pu, l  =  0.48.  In  both  cases,  the  performance  of 
the  RW-LRT,  RT-LRT,  and  filtered  energy  detector  are  all 
approximately  the  same,  i.e.  the  latter  two  detectors  are  nearly 
optimal.  When  the  transition  probabilities  of  the  d-t  random 
walk  are  asymmetric  however,  as  in  the  case  of  Fig.  12,  the  d-t 
random  walk  LRT  is  noticeably  better  than  the  filtered  energy 
detector. 


VI.  Conclusion  and  Discussion 

We  have  developed  and  compared  optimal  and  non-optimal 
detectors  under  two  single  spin  MRFM  signal  models.  Recent 
experiments  using  the  approximately  optimal  filtered  energy 
detector  have  resulted  in  the  successful  detection  of  a  single 
electron  spin.  This  is  strong  evidence  that  the  random  telegraph 
signal  model  accurately  describes  the  cantilever-single  spin 
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Fig.  10.  Simulated  ROC  curves  for  the  d-t  symmetric  random  walk  pi:±  = 
Pi, 2  =  Pu, l  =  Pu, 2  =  0.5  at  SNR  =  -39.9  dB,  T  =  60  s  for  the  matched  filter, 
RW-LRT,  RT-LRT,  filtered  energy,  amplitude,  and  unfiltered  energy  detector. 
The  RW-LRT  is  theoretically  optimal. 
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Fig.  11.  Simulated  ROC  curves  for  the  d-t  symmetric  random  walk  pii  = 
pu, 2  =  0.52,  pij2  =  Pu, l  =  0.48  at  SNR  =  -37.4  dB,  T  =  60  s.  The 
RW-LRT  is  theoretically  optimal. 


interaction.  Indeed,  recent  theoretical  results  have  surfaced  that 
rigorously  justifies  this  model  [21]. 

The  results  of  this  paper  lend  strong  theoretical  and  practical 
support  to  the  use  of  the  simple  filtered  energy  detector  for 
the  current  MRFM  single  spin  research  community.  It  has  been 
shown  that  the  existing  baseband  filtered  energy  detector  that 
is  in  current  use  is  approximately  optimal  in  the  case  of  the 
symmetric  d-t  random  telegraph  model  under  the  regime  of 
low  SNR,  long  observation  times,  and  p  close  to  1.  The  last 
condition  can  be  achieved  by  sampling  at  a  sufficiently  fast 
rate  as  compared  to  the  rate  of  random  transitions.  This  result 
has  been  extended  to  the  case  of  the  asymmetric  d-t  random 
telegraph  by  using  a  hybrid  filtered  energy/amplitude/energy 
detector.  Simulations  were  presented  showing  that  the  near 
optimality  of  the  baseband  filtered  energy  detector  extends  to 
the  case  of  the  symmetric  d-t  random  walk  model.  In  the  case 
of  the  asymmetric  d-t  random  walk,  the  filtered  energy  detector 
does  not  perform  as  well  as  the  optimal  LRT.  We  expect  that 
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Fig.  12.  Simulated  ROC  curves  for  the  d-t  asymmetric  random  walk  pp±  = 
pu  l  =  0.45,  piy 2  =  pu,2  =  0.55  at  SNR  =  -41.0  dB,  T  =  60  s.’The 
RW-LRT  is  theoretically  optimal. 


a  hybrid  detector  along  the  lines  of  that  formulated  for  the 
d-t  random  telegraph  will  perform  close  to  the  optimal  for  the 
asymmetric  d-t  random  walk.  Mathematical  analysis  of  the  d-t 
random  walk  model  will  be  presented  in  a  future  paper. 

A  new  interpretation  of  the  LRT  for  a  d-t  finite  state 
Markov  signal  in  AWGN  under  low  SNR  conditions  was 
presented.  Specifically,  the  LRT  is  approximately  the  matched 
filter  statistic  with  the  MMSE  predictor  used  in  place  of  the 
“known”  signal  values. 
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Appendix 

Let  f  (//o  •>  •  •  •  9  2/ tv — l )  denote  the  log  LRT  function  of  the  d-t 
random  telegraph;  this  is  obtained  by  taking  the  log  of  the  left- 
hand  side  of  (8).  Let  g(yo, . . . ,  2/jv-i)  be  the  filtered  energy 
detector  function  in  (6).  Let  us  analyze  the  two  functions  / 
and  g  under  the  regime  of  low  SNR  (|A/cr|  <C  1)  and  long 
observation  time  (N  >>  1). 

We  want  to  obtain  the  approximate  Taylor  series  expansion 
of  /  about  y  =  0  and  compare  that  with  g.  Define: 


the  approximate  Taylor  series  expansion  of  /  as 


/  ~  L\  +  L2a  +  L2b  +  h.o.t.,  with 
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In  (17),  “h.o.t.”  denotes  the  higher-order  terms;  specifically, 
terms  of  degree  three  of  higher.  Cm  is  a  parameter  that 
indicates  the  mismatch  between  the  transition  probabilities  p 
and  q.  In  the  symmetric  case,  p  =  q  =>  Cm  =  0,  and  one 
obtains  a  simpler  expression  for  /.  Let  /sym  be  the  function  / 
under  symmetric  transition  probabilities,  i.e.  p  =  q.  Then, 
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For  sufficiently  large  N,  it  can  be  shown  that 
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for  i  >  0.  From  (7),  a  recursive  equation  for  Oi  can  be  derived. 
Its  approximate  solution  is 


M  =  0i  +  4  E  UmM  >  o  where 
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and  r  =  p  +  q  —  1.  Note  that  p,  q  G  (0, 1)  =>  |r|  <  1.  Define 
Si  =  -fyi-  Then’ 
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By  solving  for  Ri(A)  in  terms  of  Oi  and  using  (15),  one  obtains 


i  —  2 

where  D  =  is  a  constant;  note  that  D  plays  no  role  in 
the  performance  of  the  test  statistic.  Comparing  (21)  and  (22), 
we  see  that  they  are  nearly  identical  in  form  if  a  =  2p  —  1.  If 
a  =  2p  —  1  =>  |/sym  -  g\  ~  f  Ei  Vi-  Now,  £1  [EE*  %2]  “ 
^oEjlo1  Vi]  =  A2N.  On  the  other  hand,  for  large  N, 


Ei 


N-l  i-l 


EE 

L  i=  1  7=0 

~  ga2n 


a 1  JyiVj 


—  En 


N-l  i-l 

EE 

L  2=1  7=o 


oi1  JyiVj 


(23) 


where  G  =  When  «  =  2p  -  1,  G  =  = 

4(11_p)  +  A  -  1-  For  p  close  to  1,  G  >  and  > 

^d2JV.  So  to  the  first  moment,  the  difference  of  A  E(  vl 
between  /sym  and  p  does  not  represent  a  significant  difference 
when  p  cx  1.  Under  these  conditions,  we  expect  that  the 
performance  of  the  filtered  energy  detector  and  the  d-t  random 
telegraph  LRT  to  be  similar. 

It  is  possible  to  obtain  an  approximation  to  the  d-t  random 
telegraph  LRT  that  holds  when  we  make  no  assumption  about 
p  being  equal  to  q.  When  p  ^  q,  Cm  ^  0,  and  there  are  terms 
of  the  form  rlCm  and  r2zC^  in  (18)-(20).  Since  \r\  <  1, 
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rl  —>  0  in  the  limit  as  i  — >  oo.  So  drop  these  terms  to  get: 
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where  C  =  (^)2  is  a  constant.  Define  Ca  = 

and  Ce  =  .  In  order  to  equate  the  coefficients  of  the 

cross-terms  yiyj  between  (24)  and  g  in  (22),  we  require  a  = 
r  =  p  +  q  —  1.  In  g ,  the  ratio  of  the  energy  terms  to  the 
cross-terms  is  For  r  =  a  ~  1  =>  .r-^-  ~  1/2.  The  idea 
is  to  add  the  energy  and  amplitude  statistics  to  g  so  that  all 
three  statistics  are  in  the  same  ratio  as  in  (24).  Let  #hyb  be 
the  “extended”  version  of  g ,  which  we  shall  call  the  hybrid 
filtered  energy/amplitude/energy  detector: 
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We  expect  ghyb  to  have  performance  that  is  similar  to  / 
under  the  conditions  of  large  N,  low  SNR,  and  r  ~  1. 
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We  have  detected  and  manipulated  the  naturally  occurring  statistical  polarization  in  nuclear  spin  en¬ 
sembles  using  magnetic  resonance  force  microscopy.  Using  protocols  previously  developed  for  detecting  single 
electron  spins,  we  have  measured  signals  from  ensembles  of  nuclear  spins  in  a  volume  of  roughly  (150  nm)3 
with  a  sensitivity  of  roughly  2000  net  spins  in  a  2.5  h  averaging  window.  Three  systems  have  been  studied,  19F 
nuclei  in  CaF2,  and  XH  nuclei  (protons)  in  both  polymethylmethacrylate  and  collagen,  a  naturally  occurring 
protein.  By  detecting  the  statistical  polarization,  we  not  only  can  work  with  relatively  small  ensembles,  but  we 
eliminate  any  need  to  wait  a  longitudinal  relaxation  time  Tx  to  polarize  the  spins.  We  have  also  made  use  of  the 
fact  that  the  statistical  polarization,  which  can  be  considered  a  form  of  spin  noise,  has  a  finite  correlation  time. 
A  method  similar  to  one  previously  proposed  by  Carlson  et  al.  [Bull.  Am.  Phys.  Soc.  44,  541  (1999)]  has  been 
used  to  suppress  the  effect  of  the  statistical  uncertainty  and  extract  meaningful  information  from  time-averaged 
measurements.  By  implementing  this  method,  we  have  successfully  made  nutation  and  transverse  spin  relax¬ 
ation  time  measurements  in  CaF2  at  low  temperatures. 

DOI:  10. 1103/PhysRevB. 72. 024413  PACS  number(s):  76.60. -k,  05.40. -a,  07.55.-w 


I.  INTRODUCTION 

Magnetic  resonance  force  microscopy  (MRFM)  combines 
ultrasensitive  force  detection  with  the  principles  of  magnetic 
resonance  to  probe  small  volumes  of  electron  and  nuclear 
spins.1-6  If  the  goal  of  single  nuclear  spin  sensitivity  can  be 
reached,  the  technique  may  ultimately  allow  for  three- 
dimensional  molecular  imaging  with  atomic  resolution.2  It 
might  also  function  as  a  readout  scheme  for  quantum 
computation.7,8  While  the  technique  is  still  far  from  achiev¬ 
ing  these  challenging  goals,  steady  progress  has  been  made, 
as  evidenced  by  the  recent  demonstration  of  detection  of  a 
single  electron  spin  with  a  spatial  resolution  of  25  nm.9 

The  measurement  protocols  previously  developed  for 
small  ensembles  of  electron  spins  did  not  require  a  mean 
spin  polarization,  but  rather  took  advantage  of  the  naturally 
occurring  statistical  polarization.  For  an  ensemble  of  N  spins, 
this  polarization  is  of  order  \/a,  and  can  actually  exceed  the 
thermal  (Boltzmann)  polarization  for  small  N.  One  of  the 
purposes  of  the  current  work  is  to  demonstrate  that  these 
protocols  are  applicable  to  detecting  the  statistical  polariza¬ 
tion  of  nuclear  spins  as  well  as  electron  spins.  This  approach 
is  particularly  valuable  for  nuclear  spins,  since  the  time  re¬ 
quired  to  polarize  a  nuclear  spin  ensemble  can  be  exceed¬ 
ingly  long  at  low  temperatures.  In  this  work,  we  verify  the 
technique  by  studying  ensembles  of  order  A/=10 8  nuclear 
spins  in  a  volume  of  roughly  (150  nm),3  where  we  detect  and 
manipulate  the  statistical  polarization  of  order  \^V=104  net 
spins. 

A  second  purpose  of  this  work  is  to  demonstrate  a  scheme 
to  extract  meaningful  information  in  the  presence  of  the  spin 
noise10-13  associated  with  the  statistical  polarization.  It  has 
recently  been  proposed  that  manipulation  of  such  fluctuating 
statistical  ensembles  should  be  possible  by  making  use  of  the 
finite  correlation  time  of  the  fluctuations.14,15  We  have  imple¬ 
mented  a  version  of  this  scheme  using  MRFM,  and  have 


performed  spin  manipulations  on  statistical  ensembles  of  19F 
and  'll  nuclei.  Despite  the  fact  that  the  magnitude  and  even 
the  sign  of  the  polarization  fluctuates  over  time,  this  ap¬ 
proach  has  allowed  us  to  perform  nutations  and  transverse 
spin  relaxation  measurements. 


II.  EXPERIMENTAL  DETAILS  AND  PROTOCOLS 

Two  species  of  nuclear  spins  were  studied  in  three  differ¬ 
ent  samples.  First,  we  studied  19F  in  a  polished  single  crystal 
of  CaF2  (99.99%  pure).16  We  also  studied  protons  in  two 
different  systems,  the  polymer  polymethylmethacrylate 
(PMMA),  and  the  triple-helix  protein  collagen.  PMMA  has 
previously  been  studied  with  MRFM  using  the  Boltzmann 
polarization.17  The  collagen  was  chosen  to  demonstrate  the 
applicability  of  MRFM  to  the  study  of  biomolecules.  The 
PMMA  sample  was  a  200  nm  thick  film  that  was  spun  onto  a 
silicon  wafer  and  baked  at  175  °C.  It  was  then  metal  coated 
with  75  A  Cr/200  A  Au  to  reduce  charging  effects.  The  col¬ 
lagen  sample  was  prepared  from  a  commercially  available 
solution  of  rat-tail  collagen  in  acetic  acid  (Type  1,  BD  Bio¬ 
sciences).  A  small  drop  of  the  solution  was  placed  onto  a 
silica  substrate  and  allowed  to  air  dry.  The  sample  was  fur¬ 
ther  dried  in  vacuum  at  room  temperature  for  4  days  before 
use,  but  otherwise  no  special  preparations  were  made.  All 
three  systems  have  expected  spin  densities  in  the  range  of 
4-6  X  1022  spins/cm3. 

As  shown  in  Fig.  1  and  described  in  Ref.  12,  our  MRFM 
technique  uses  a  magnetic  tip  mounted  on  the  end  of  a  ver¬ 
tically  oriented  ultrasensitive  cantilever.  The  tip  generates  a 
strong  magnetic  field  gradient  (~1  G/nm).  A  rf  field  Bx  at 
frequency  o)rf  is  used  to  excite  magnetic  resonance  within  a 
thin  resonant  slice  where  the  condition  B0(x,y ,z)  =  (orf/ y  is 
met.  Here  B0  is  due  to  the  tip  field  plus  an  optional  external 
field  BexP  and  y  is  the  gyromagnetic  ratio  (4.2  kHz/G  for 
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FIG.  1 .  Schematic  of  the  MRFM  experiment.  Left-right  statisti¬ 
cal  imbalance  of  the  spin  polarization  within  the  resonant  slice 
gives  rise  to  forces  on  the  magnetic  tip. 

protons  and  4.0  kHz/G  for  19F).  In  the  experiment  on  CaF2, 
for  example,  the  rf  frequency  was  50.8  MHz,  so  that  the 
resonance  condition  was  B0(x,y,z)  =  12  680  G.  The  rf  field 
Bi~  15  G  was  generated  using  a  300  jmm  diameter  coil  that 
was  part  of  an  LC  resonant  circuit.  The  microscope  was  op¬ 
erated  in  a  small  vacuum  chamber  within  the  bore  of  a  su¬ 
perconducting  magnet.  The  temperature  of  the  microscope 
with  the  rf  coil  energized  was  roughly  7  K. 

The  same  single-crystal  silicon  cantilever  was  used  for  all 
three  samples.  It  consisted  of  a  roughly  100  nm  thick,  90  jmm 
long  shaft,  with  a  1  jmn l  X  10  jum  X  4  jmm  silicon  mass  on  the 
end  to  suppress  the  motion  of  the  higher  order  modes.18  The 
cantilever  spring  constant  k  was  roughly  8.6  X  10-5  N/m.  In 
zero  magnetic  field,  the  cantilever  resonant  frequency  fc  was 
3830,  and  the  quality  factor  Q  was  roughly  60  000  at  4.2  K. 
A  SmCo  particle  was  attached  to  the  end  of  the  cantilever 
with  its  magnetic  axis  carefully  oriented  in  a  magnetic  field 
to  within  5°  of  the  cantilever  axis,  and  then  shaped  with  a 
focused  ion  beam  to  submicron  dimensions  to  reduce  its  total 
moment.  This  preparation  was  necessary  to  prevent  exces¬ 
sive  cantilever  bending  in  the  applied  magnetic  field.  The  tip 
of  the  particle  was  further  shaped  to  a  roughly  250  nm  wide 
apex,  resulting  in  field  gradients  in  the  range  of  1  G/nm.  In 
a  field  of  12  000  G,  the  cantilever  frequency  increased  to 
about  4200  Hz  due  to  field-induced  restoring  torque,  and  the 
Q  was  reduced  to  approximately  5000  due  to  various  mag¬ 
netic  loss  mechanisms.19,20  This  undesirable  effect,  while  not 
a  major  problem,  can  be  eliminated  through  the  use  of  an 
alternative  cantilever  orientation.21 

The  spin-manipulation  and  detection  scheme  is  based  on 
the  OSCAR  protocol  used  previously  for  electron  spins,9,12,22 
where  OSCAR  stands  for  oscillating  cantilever-driven  adia¬ 
batic  reversals.  The  cantilever  is  self-oscillated  at  its  natural 
resonant  frequency  using  a  positive  feedback  loop.  As  the 
cantilever  position  oscillates  according  to  xc(t)=xpk  cos (a)ct), 
the  field  B0  at  a  given  sample  location  is  modulated  because 
of  the  field  gradient  G=dB0/dx  from  the  tip.  In  the  language 
of  magnetic  resonance,  the  effective  field  B in  the  rotating 
frame  is  given  by12,23 

Beif=Blx  +  Gxc(t)z.  (1) 

For  Gxpk>Bu  J$eff  will  change  from  the  +£  to  —z.  direction  as 
the  cantilever  oscillates.  As  long  as  the  adiabatic  condition  is 


FIG.  2.  Timing  diagram  for  conventional  iOSCAR  protocol.  As 
the  cantilever  oscillates,  the  z  component  of  spin  is  modulated  in 
response  to  the  cantilever  motion,  except  when  the  rf  field  is  inter¬ 
rupted.  These  interruptions  last  for  precisely  one-half  an  oscillation 
period  and  occur  every  96  cantilever  cycles,  leading  to  a  modula¬ 
tion  in  the  cantilever  frequency. 

met,  the  spins  will  be  locked  (or  antilocked)  to  the  effective 
field,  and  will  thereby  be  synchronously  inverted  with  each 
cantilever  oscillation.23  The  backaction  force  of  the  spins  on 
the  cantilever  gives  rise  to  a  slight  cantilever  frequency  shift 
Sfc,  whose  sign  depends  on  whether  the  spins  are  locked  or 
antilocked  to  the  effective  field. 

Rather  than  measure  a  dc  frequency  shift,  we  impose  a 
modulation  on  the  signal  using  an  “interrupted”  version  of 
the  protocol  referred  to  as  iOSCAR.  As  shown  in  Fig.  2,  we 
turn  off  the  rf  power  for  one-half  of  a  cantilever  cycle  every 
96  cycles,  which  skips  one  adiabatic  reversal,  which,  in  turn, 
results  in  a  change  in  the  sign  of  the  frequency  shift  Sfc.  This 
imposes  a  modulation  on  the  frequency  shift  at  one-half  the 
rf  interrupt  frequency,  giving  the  signal  a  very  distinctive 
signature.  The  frequency  shift  signal  is  analyzed  using  either 
a  Fourier  power  spectrum  or  lock-in  detector  implemented  in 
software.  For  statistical  polarization,  the  mean  lock-in  signal 
is  zero,  so  we  detect  the  variance,  or  signal  “energy.”  A  zero- 
base  line  signal  is  constructed  by  subtracting  the  in-phase 
and  quadrature  variances.  The  optimal  signal-to-noise  ratio 
will  be  obtained  when  the  lock-in  time  constant  is  properly 
matched  to  the  spin-lock  correlation  time  rm.24  We  use  a 
bank  of  filters  to  obtain  the  signal  as  a  function  of  measure¬ 
ment  bandwidth,  which  is  used  to  deduce  rm.9 

In  the  limit  of  Gxpk>Bh  the  cantilever  frequency  shift 
due  to  the  spins  is  given  by25 


Sfc  =  ( 2fchrkxpk)meffG 

(2a) 

=  {2fc/TTkxpk)Fs, 

(2b) 

where  mejj  is  the  component  of  magnetic  moment  along  the 
effective  field,  and  Fs  =  mefjG  represents  the  peak  force  ex¬ 
erted  by  the  spins  on  the  cantilever.  Equation  (2a)  can  be 
used  to  extract  the  effective  number  of  spins  contributing  to 
the  signal.  Alternatively,  Eq.  (2b)  can  be  used  to  convert  the 
signal  to  a  net  force,  which  is  useful  when  comparing  signals 
at  different  cantilever  amplitudes  xph  for  example. 
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FIG.  3.  Power  spectral  density  of  the  frequency  shift  signal  for 
the  CaF2  sample.  An  external  field  of  12  600  G  was  used  to  bring 
the  sample  into  resonance,  resulting  in  the  statistical  spin  signal, 
seen  as  the  peak  at  22  Hz.  The  peak  width  of  1.9  Hz  indicates  a 
correlation  time  of  roughly  170  ms.  The  estimated  tip-sample  sepa¬ 
ration  was  —  150nm.  The  thermal  noise  limit  for  our  cantilever 
with  2  =  5000  is  indicated  by  the  dotted  line.  The  background  noise 
level  is  higher  than  this  level  because  of  sample-induced  excess 
frequency  noise,  which  has  a  1  If- like  spectrum. 

III.  MEASUREMENTS  OF  NUCLEAR  STATISTICAL 
POLARIZATION 

The  power  spectral  density  of  the  iOSCAR  signal  for  the 
CaF2  sample  is  shown  in  Fig.  3.  The  rf  field  was  turned  off 
every  96  cantilever  cycles  (fint=fc/96  —  44  Hz),  resulting  in 
a  signal  peak  at  22  Hz.  The  spectral  width  of  the  peak  is 
about  1.9  Hz,  indicating  a  spin-lock  correlation  time  rm  of 
0.17  s.  This  time  is  closely  related  to  Tlp ,  the  spin-lattice 
relaxation  time  in  the  rotating  frame.12,23  It  can  be  affected 
by  numerous  parameters,  such  as  the  strength  of  the  rf  field, 
magnetic  noise  at  the  Rabi  frequency  due  to  higher  order 
cantilever  modes,26,27  violation  of  the  adiabatic  condition, 
and  the  OSCAR  protocol  itself.12,22  The  latter  two  effects 
will  be  discussed  later  in  the  paper. 

Scanning  the  external  magnetic  field  Bext  moves  the  reso¬ 
nant  slice  relative  to  the  tip,  resulting  in  a  signal  when  the 
slice  penetrates  into  the  sample.  We  plot  the  signal  as  a  func¬ 
tion  of  external  field  in  Fig.  4(a),  where  the  data  represent  the 
signal  power  from  the  lock-in  detector  in  mHz2.  A  zero  base 
line  signal  is  observed  (within  experimental  error),  with  a 
large  peak  present  in  the  range  of  —  12  000- 13  000  G,  where 
the  width  of  the  peak  is  inhomogeneously  broadened  due  to 
the  field  gradient  from  the  tip.  Assuming  an  estimated  field 
gradient  of  1  G/nm,  the  maximum  signal  power  of 
1 100  mHz2  corresponds  to  meff=  1.7  X  10-19  emu,  or  of  order 
10  000  net  nuclear  spins.  The  noise  level  obtained  with  an 
averaging  time  of  2.5  h  was  —25  mHz2,  corresponding  to 
roughly  2000  net  nuclear  spins.  From  Fig.  3,  the  integrated 
base  line  noise  in  the  natural  1.9  Hz  bandwidth  was 
—5000  mHz2,  corresponding  to  roughly  30  000  net  spins. 

Since  we  are  not  measuring  the  Boltzmann  polarization, 
there  is  no  need  to  wait  a  longitudinal  relaxation  time  Tl 
between  measurements.  For  our  CaF2  sample,  the  measured 
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FIG.  4.  MRFM  signal  vs  externally  applied  magnetic  field,  (a) 
19F  in  CaF2.  The  cantilever  oscillation  amplitude  was  16  nm,  and 
the  applied  rf  frequency  was  50.8  MHz,  resulting  in  magnetic  reso¬ 
nance  when  the  total  field  is  equal  to  12  680  G.  The  continuous 
curve  is  calculated  from  a  model  based  on  a  spherical  tip  with  a 
diameter  of  440  nm.  In  order  to  fit  the  model  to  the  data,  we  had  to 
assume  a  shift  of  +250  G  in  the  field  value  at  the  location  of  the 
spins.  This  could  be  due  to  remanent  magnetization  in  the  body  of 
our  microscope,  or  a  slight  (2%)  error  in  our  field  measurement,  (b) 
Protons  in  collagen.  The  lower  rf  frequency,  along  with  the  higher 
gyromagnetic  ratio  for  protons,  result  in  a  shift  in  the  resonant  field 
condition  compared  to  19F. 

Tx  was  over  800  s  at  room  temperature,  and  can  be  hours, 
days  or  longer  in  pure  samples  at  low  temperatures.28,29  This 
consideration  illustrates  the  significant  advantage  that  can 
potentially  be  derived  by  working  with  the  statistical  polar¬ 
ization. 

The  MRFM  signals  from  ]H  nuclei  in  PMMA  and  col¬ 
lagen  were  also  readily  detected.  We  plot  the  iOSCAR  signal 
power  as  a  function  of  Bext  for  collagen  in  Fig.  4(b).  The 
signal  is  seen  in  the  range  of  Bext~  10  000-11  000  G.  As 
expected,  the  field  range  is  lower  than  in  the  case  of  the 
CaF2,  partly  because  of  the  5%  higher  gyromagnetic  ratio  in 
compared  to  19F,  and  partly  because  the  rf  frequency  was 
slightly  lower.  In  addition,  the  exact  shape  and  position  of 
the  peak  will  depend  on  the  distance  between  the  tip  and 
sample,  which  may  have  been  somewhat  less  in  the  collagen 
sample. 

The  rm  correlation  times  extracted  from  the  lock-in  detec¬ 
tor  data  were  disappointingly  short,  with  rm  —  200  ms  for 
CaF2  and  rm—  100  ms  for  collagen.  In  the  case  of  PMMA, 
the  correlation  times  were  less  than  50  ms,  and  it  was  not 
possible  to  make  an  accurate  measurement.  In  comparison, 
values  of  rm  up  to  several  seconds  have  been  measured  for 
E'  centers  (electron  spins).9,12  The  shorter  times  observed 
here  may  be  at  least  partly  due  to  issues  with  these  particular 
materials:  CaF2  requires  a  large  field  Bx  to  overcome  the 
local  fields,30  for  example,  while  in  PMMA,  it  is  known  that 
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FIG.  5.  Signal  power  vs  cantilever  oscillation  amplitude  for 
CaF2.  The  open  circles  are  for  lower  rf  power  (Bl  =  15.5  G)  and  the 
filled  circles  are  for  higher  rf  power  (51  =  18.1  G).  If  the  adiabatic 
reversals  were  equally  effective  under  all  conditions,  the  signal 
power  would  rise  linearly  from  zero  as  a  function  of  oscillation 
amplitude.  In  the  case  shown,  the  signal  is  diminished  at  the  lowest 
amplitudes  by  the  fact  that  the  reversals  are  incomplete.  At  the 
higher  amplitudes,  the  adiabatic  condition  is  no  longer  met,  which 
also  reduces  the  signal. 

tunneling  of  the  methyl  groups  can  occur  over  a  wide  fre¬ 
quency  range  and  can  affect  the  spin  relaxation  behavior.31 
At  the  same  time,  it  is  also  likely  that  the  iOSCAR  protocol 
itself  played  a  role. 

For  the  iOSCAR  protocol  to  work  effectively,  two  condi¬ 
tions  must  be  met.  First,  the  field  B0  must  be  swept  well  off 
resonance;  i.e.,  Gxpk>Bl.  As  seen  by  Eq.  (1),  this  ensures 
that  the  effective  field,  and  thus  the  spins,  are  very  nearly  in 
the  z  direction  when  the  cantilever  reaches  its  extremum. 
Otherwise,  each  time  the  rf  field  is  interrupted,  the  projected 
spin  component  along  the  effective  field  will  be  reduced. 
Second,  the  adiabatic  condition  dB0/ dt<yB\  must  be 
satisfied.23  A  large  value  of  xph  which  drives  the  field  B0 
farther  off  resonance,  helps  to  satisfy  the  first  condition. 
However,  it  hurts  the  second  condition,  since  dB0/  dt^xpk. 

We  observed  that  there  was  a  fairly  restricted  range  of 
cantilever  oscillation  amplitudes  for  which  the  iOSCAR  pro¬ 
tocol  worked  well.  This  is  illustrated  in  Fig.  5,  where  the 
signal  power  vs  oscillation  amplitude  xpk  is  plotted  for  two 
different  values  of  the  rf  field.  For  low  oscillation  ampli¬ 
tudes,  the  signal  increased  as  xpk  was  increased.  This  is  partly 
because  of  the  larger  volume  of  spins  being  swept  through  by 
the  resonant  slice,  but  also  because  the  field  B0  was  being 
swept  farther  off  resonance,  leading  to  more  complete  adia¬ 
batic  inversions.  As  xpk  was  further  increased,  the  signal 
reached  a  peak  and  then  diminished,  presumably  due  to  vio¬ 
lation  of  the  adiabatic  condition.  We  see  that  this  decrease 
occurred  for  adiabaticity  values,  defined  as  (dB0/ dt)  /  yB\, 
greater  than  —8%.  With  a  larger  value  of  Bx  (upper  curve), 
the  signal  was  greater,  and  the  optimum  value  for  xpk  was 
somewhat  larger  as  well. 

The  conflicting  requirements  of  iOSCAR  and  the  adia¬ 
batic  condition  can  be  eased  in  at  least  two  ways.  First,  the 
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FIG.  6.  Pulse  sequences  for  use  with  modified  iOSCAR  that 
allow  for  manipulation  of  the  statistical  polarization,  (a)  Nutations. 
The  rf  field  is  interrupted  for  one  full  cantilever  cycle  every 
96  cycles.  Within  the  interruptions,  the  rf  field  is  pulsed  at  the  zero 
crossings.  Those  spins  that  are  precisely  on  resonance  will  be  begin 
to  precess  at  the  Rabi  frequency  about  the  effective  field,  which  is 
in  the  x  direction.  Thus  the  z  component  will  be  a  sinusoidal  func¬ 
tion  of  the  pulse  width  T.  (b)  Transverse  relaxation.  Two  tt/2  pulses 
are  inserted  a  time  r  apart.  For  short  r,  the  effect  is  identical  to  a 
single  77  pulse.  For  long  r,  the  spins  are  left  in  a  random  direction  in 
the  x-z  plane,  and  the  spins  are  equally  likely  to  be  in  the  +z  and  -z 
directions  when  the  field  is  turned  back  on. 

condition  Gxpk>Bx  can  be  satisfied  by  reducing  Bx  as  the 
cantilever  reaches  its  extremum.  This  type  of  rf  pulse  shap¬ 
ing  was  in  fact  performed  in  these  measurements,  and  in¬ 
creased  the  signal  by  up  to  nearly  a  factor  of  2.  Second, 
decreasing  the  cantilever  resonant  frequency  would  allow  for 
greater  values  of  xpk  without  violating  the  adiabatic  limit. 
There  may  be  lower  practical  limits  on  cantilever  frequency, 
however,  due  to  factors  such  as  environmental  vibrational 
noise. 

IV.  NUCLEAR  SPIN  MANIPULATION  VIA  MODIFIED 
iOSCAR 

The  above  results  verify  that  the  protocols  developed  for 
detecting  statistically  polarized  ensembles  of  electron  spins 
are  also  applicable  to  nuclear  spins.  In  addition  to  detection, 
the  statistical  polarization  can  also  be  manipulated  into  de¬ 
sired  states  through  proper  modification  of  the  protocol.  For 
example,  we  have  previously  demonstrated  manipulation  of 
the  statistical  polarization  of  electron  spin  ensembles  using 
iOSCAR  with  extra  rf  interrupts  in  order  to  the  rectify  the 
spin  fluctuations.32 

In  this  section,  we  describe  additional  manipulation 
schemes  that  essentially  compare  the  net  polarization  before 
and  after  a  series  of  rf  pulses.  Because  of  its  statistical  na¬ 
ture,  the  initial  state  may  be  different  each  time  the  measure¬ 
ment  is  performed.  Nevertheless,  the  pulses  will  affect  the 
net  polarization  in  some  deterministic  way  that  can  be  de¬ 
tected  due  to  the  finite  correlation  time  of  the  polarization. 
The  scheme  is  functionally  very  similar  to  a  previously  pro¬ 
posed  scheme  called  CONQUEST  (correlated  observations 
narrow  quantum  uncertainty,  enhancing  spectroscopic  tran¬ 
sients),  which  is  based  on  a  second  order  correlation  function 
of  the  signal  before  and  after  the  pulse  sequence. 14,15,33 

Our  version  of  the  detection  protocol,  illustrated  in  Fig.  6, 
is  based  on  the  iOSCAR  protocol,  with  one  key  difference. 
The  periodic  interruptions  in  the  rf  field  are  now  for  one  full 
cantilever  cycle,  rather  than  one  half  cycle.  When  the  rf  field 
is  interrupted  for  one  full  cantilever  cycle,  a  spin  is  left  in 
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FIG.  7.  Nutations  of  the  statistical  polarization  in  CaF2.  The 
modified  iOSCAR  signal  is  shown  as  a  function  of  the  pulse  width 
tp.  The  oscillatory  nature  of  the  signal  is  due  to  the  precession  of  the 
spins  about  the  effective  field  at  the  Rabi  frequency.  The  period  of 
the  oscillations  indicates  a  field  strength  ^  =  17  G. 

exactly  the  same  phase  relative  to  the  cantilever  motion  as 
when  it  started.  Therefore  no  modulation  of  the  frequency 
shift  should  occur.  However,  suppose  that  at  every  full-cycle 
rf  interrupt,  a  rf  pulse  is  inserted  precisely  on  resonance,  as 
shown  in  Fig.  6(a).  The  spins  will  nutate  around  the  effective 
field,  which  in  this  case  is  in  the  x  direction.  If  the  pulse 
width  tp  corresponds  to  a  77  pulse,  the  spins  will  end  up 
inverted,  and  the  situation  will  be  equivalent  to  normal 
iOSCAR:  A  modulated  frequency  shift  will  be  observed  at 
half  the  interrupt  frequency,  which  can  be  synchronously  de¬ 
tected  using  a  lockin  amplifier.  If  the  pulse  width  is  2  77,  then 
spins  are  left  noninverted,  and  the  modulation  should  again 
disappear.  Thus  the  signal  should  show  oscillatory  behavior, 
with  a  period  given  by  2ttI  yBx. 

This  protocol  has  allowed  us  to  make  our  detection  of 
nutations  using  a  statistically  polarized  ensemble.  The  results 
for  CaF2  are  shown  in  Fig.  7.  For  the  case  in  which  no 
nutation  pulse  is  inserted  during  the  rf  interrupt  (^=0),  the 
signal  is  nearly  zero,  as  expected.  With  a  pulse  inserted,  the 
signal  shows  clear  oscillatory  behavior  as  a  function  of  pulse 
width,  with  a  period  of  roughly  15  juls.  The  measurement  of 
the  period  allows  for  an  absolute  determination  of  the  rf  field 
strength,  in  this  case  Bl  =  ll  G. 

There  is  a  close  connection  between  this  detection  scheme 
and  the  CONQUEST  scheme,  as  suggested  by  the  following 
argument:  Consider  a  model  for  the  synchronous  lock-in  de¬ 
tection  in  which  the  basic  signal  8fc(t)  is  multiplied  by  a 
square  wave  of  period  IT  and  then  averaged.  The  resulting 
output  signal  then  is  made  up  of  terms  of  the  form 
Sfc(t)-8fc(t+T).  (Here  we  have  assumed  that  the  square 
wave  has  no  dc  component.  We  also  ignore  correlations 
for  times  ^T,  so  that  we  can  break  up  the  time  record 
into  individual  records  of  length  2  T,  each  of  which  is 
equivalent.)  In  this  model,  the  mean  square  signal  is 
simply  the  time  average  of  \8fc(t)~  8fc(t+T)]2,  or 
2( Sfc(t)2)-2(Sfc(t)Sfc(t+T )),  where  (•••)  denotes  the  time 
average.  In  normalized  form,  this  is  simply  1-A,  where  A  is 
the  same  correlation  function  (Iz(t)Iz(t+T))  measured  by 


FIG.  8.  Modified  iOSCAR  signal  versus  pulse  delay  r.  An  ex¬ 
ponential  decay  to  the  baseline  is  observed,  indicating  a  transverse 
relaxation  time  T2  =  4  juls. 

CONQUEST  in  its  most  straightforward  version.33  In  either 
case,  the  technique  relies  upon  the  systematic  effect  on 
the  signal  (such  as  inversion)  caused  by  the  inserted  pulses. 
This  results  in  a  time-averaged  correlation  function 
( 8fc(t)8fc(t+T ))  that  has  a  nonzero  value,  even  if  the  initial 
values  of  8fc(t )  are  completely  random.  Note  that  since  the 
modified  iOSCAR  scheme  imposes  a  modulation  on  the  sig¬ 
nal,  it  has  the  added  benefit  of  reducing  the  influence  of  low 
frequency  noise. 

A  slightly  different  pulse  sequence,  illustrated  in  Fig.  6(b), 
has  been  used  to  measure  the  inhomogeneous  transverse  re¬ 
laxation  time  T2.  Two  it  1 2  pulses  are  inserted  with  variable 
spacing  r  between  them.  The  signal  as  a  function  of  r  for 
CaF2  is  shown  in  Fig.  8.  As  expected,  the  signal  is  a  maxi¬ 
mum  for  small  r.  This  situation  is  equivalent  to  inserting  a 
single  77  pulse,  causing  the  spin  polarization  after  the  pulse 
sequence  to  be  anticorrelated  to  the  signal  before  the  se¬ 
quence.  As  seen  with  the  nutation  pulses,  the  result  is  a  maxi¬ 
mum  in  the  iOSCAR  signal.  For  large  r,  the  polarization 
after  the  pulse  sequence  is  completely  uncorrelated  to  that 
before,  due  to  transverse  relaxation,  and  one  should  get  a 
signal  that  is  one-half  the  maximum  signal.  (Alternatively, 
for  small  r,  the  final  polarization  is  the  negative  of  the  initial 
polarization,  so  the  difference  is  twice  the  initial  polarization. 
For  large  r,  the  final  polarization  has  decayed  to  zero,  so  the 
difference  is  half  the  difference  in  the  short  r  case.)  This 
behavior  is  seen  in  Fig.  8,  with  a  maximum  signal  for  short  r, 
and  an  exponential  decay  to  a  base  line  that  is  roughly  half 
the  initial  value.  From  the  exponential  decay,  we  find  that 
that  the  transverse  relaxation  time  T2  =  4  juls.  This  short  time 
is  a  reflection  of  the  highly  inhomogeneous  field  from  the  tip, 
which  leads  to  a  loss  of  phase  coherence  of  the  spins  in 
physically  different  locations. 

V.  CONCLUSIONS 

In  summary,  we  have  successfully  applied  the  protocols 
developed  for  MRFM  of  electron  spins  to  the  case  of  nuclear 
spins.  By  performing  iOSCAR  to  detect  the  statistical  polar¬ 
ization  in  ensembles  of  nuclear  spins,  we  have  demonstrated 
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detection  with  sensitivity  on  the  order  of  2000  net  nuclear 
spins  in  our  2.5  h  averaging  window.  The  statistical  detection 
avoids  the  problem  of  the  long  Tx  times  that  exists  for  many 
nuclear  spin  systems  at  low  temperatures.  We  have  applied 
the  technique  to  the  system  CaF2,  and  have  demonstrated  our 
use  of  MRFM  on  a  biomolecule.  We  have  also  demonstrated 
the  ability  to  manipulate  the  naturally  occurring  statistical 
polarization  in  order  to  make  nonequilibrium  measurements, 
such  as  relaxation  times  and  nutation  measurements.  Given  a 
roughly  1000X  improvement  in  signal-to-noise  ratio,  which 
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remains  a  considerable  challenge,  these  techniques  should  in 
principle  be  extendible  to  individual  nuclear  spins. 

ACKNOWLEDGMENTS 

We  thank  M.  Sherwood,  C.  M.  Jefferson,  J.  Logan,  and  C. 
Rettner  for  discussions  and  technical  assistance.  We  grate¬ 
fully  acknowledge  support  from  the  DARPA  Three- 
Dimensional  Atomic-Scale  Imaging  program  administered 
through  the  US  Army  Research  Office. 


1J.  A.  Sidles,  Phys.  Rev.  Lett.  68,  1124  (1992). 

2J.  A.  Sidles,  J.  L.  Garbini,  K.  J.  Bruland,  D.  Rugar,  O.  Zueger,  S. 

Hoen,  and  C.  S.  Yannoni,  Rev.  Mod.  Phys.  67,  249  (1995). 

3D.  Rugar,  O.  Zuger,  S.  Hoen,  C.  S.  Yannoni,  H.  M.  Vieth,  and  R. 

D.  Kendrick,  Science  264,  1560  (1994). 

4Z.  Zhang,  P.  C.  Hammel,  and  P.  E.  Wigen,  Appl.  Phys.  Lett.  68, 
2005  (1996). 

5K.  R.  Thurber,  L.  E.  Harrell,  R.  Fainchtein,  and  D.  D.  Smith, 
Appl.  Phys.  Lett.  80,  1794  (2002). 

6S.  R.  Garner,  S.  Kuehn,  J.  M.  Dawlaty,  N.  E.  Jenkins,  and  J.  A. 

Marohn,  Appl.  Phys.  Lett.  84,  5091  (2004). 

7D.  P.  DiVincenzo,  Phys.  Rev.  A  51,  1015  (1995). 

8G.  P.  Berman,  G.  D.  Doolen,  P.  C.  Hammel,  and  V.  I.  Tsifrinov- 
ich,  Phys.  Rev.  B  61,  14694  (2000). 

9D.  Rugar,  R.  Budakian,  H.  J.  Mamin,  and  B.  W.  Chui,  Nature 
(London)  430,  329  (2004). 

10 F.  Bloch,  Phys.  Rev.  70,  460  (1946). 

11 T.  Sleator,  E.  L.  Hahn,  C.  Hilbert,  and  J.  Clarke,  Phys.  Rev.  Lett. 
55,  1742  (1985). 

12 H.  J.  Mamin,  R.  Budakian,  B.  W.  Chui,  and  D.  Rugar,  Phys.  Rev. 
Lett.  91,  207604  (2003). 

13  S.  A.  Crooker,  D.  G.  Rickel,  A.  V.  Balatsky,  and  D.  L.  Smith, 
Nature  (London)  431,  49  (2004). 

14 P.  L.  Carlson,  L.  Madsen,  G.  M.  Leskowitz,  and  D.  Weitekamp, 
Bull.  Am.  Phys.  Soc.  44,  541  (1999). 

15  P.  L.  Carlson,  L.  A.  Madsen,  G.  M.  Leskowitz,  and  D.  P. 

Weitekamp,  U.S.  Patent  Nos.  6,078,082  and  6,081,119  (2000). 
16 Corning  (formerly  Optovac),  North  Brookfield  MA. 

17  J.  H.  Choi,  U.  M.  Mirsaidov,  C.  W.  Miller,  Y.  J.  Lee,  S.  Guchhait, 

M.  D.  Chabot,  W.  Lu,  and  J.  T.  Markert,  Proc.  SPIE  5389,  399 
(2004). 

18  B.  W.  Chui,  Y.  Hishinuma,  R.  Budakian,  H.  J.  Mamin,  T.  W. 


Kenny,  and  D.  Rugar,  Technical  Digest  of  the  12th  International 
Conference  on  Solid-State  Sensors  and  Actuators  (Transduc- 
ers’03 ),  (IEEE,  Boston,  MA,  2003),  p.  1120. 

19 B.  C.  Stipe,  H.  J.  Mamin,  T.  D.  Stowe,  T.  W.  Kenny,  and  D. 
Rugar,  Phys.  Rev.  Lett.  86,  2874  (2001). 

20  T.  N.  Ng.,  N.  E.  Jenkins,  and  J.  A.  Marohn  (private  communica¬ 

tion). 

21  J.  A.  Marohn,  R.  Fainchtein,  and  D.  D.  Smith,  Appl.  Phys.  Lett. 

73,  3778  (1998). 

22 B.  C.  Stipe,  H.  J.  Mamin,  C.  S.  Yannoni,  T.  D.  Stowe,  T.  W. 

Kenny,  and  D.  Rugar,  Phys.  Rev.  Lett.  87,  277602  (2001). 

23  C.  P.  Slichter,  Principles  of  Magnetic  Resonance  (Springer, 
Heidelberg,  1996). 

24 M.  Ting,  A.  O.  Hero,  D.  Rugar,  C.-Y.  Yip,  J.  A.  Fessler,  preprint 
at  http://xxx.lanl.gov/abs/quant-ph/0312139  (2003). 

25  G.  P.  Berman,  D.  I.  Kamenev,  and  V.  I.  Tsifrinovich,  Phys.  Rev.  A 

66,  023405  (2002). 

26  D.  Mozyrsky,  I.  Martin,  D.  Pelekhov,  and  P.  C.  Hammel,  Appl. 

Phys.  Lett.  82,  1278  (2003). 

27  G.  P  Berman,  V.  N.  Gorshkov,  D.  Rugar,  and  V.  I.  Tsifrinovich, 

Phys.  Rev.  B  68,  094402  (2003). 

28 B.  T.  Beal  and  J.  Hatton,  Phys.  Lett.  15,  210  (1965). 

29 P.  L.  Kuhns,  P.  C.  Hammel,  O.  Gonen,  and  J.  S.  Waugh,  Phys. 
Rev.  B  35,  4591  (1987). 

30 K.  Wago,  O.  Zuger,  R.  D.  Kendrick,  C.  S.  Yannoni,  and  D.  Rugar, 
J.  Vac.  Sci.  Technol.  B  14,  1197  (1996). 

31  A.  J.  Horsewill,  Prog.  Nucl.  Magn.  Reson.  Spectrosc.  35,  359 
(1999). 

32 R.  Budakian,  H.  J.  Mamin,  B.  W.  Chui,  and  D.  Rugar  (unpub¬ 
lished). 

33  L.  A.  Madsen,  Ph.D.  thesis,  Caltech,  2002. 


024413-6 


